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In a photosynthetic reaction 
centre an electron is stepwise 
transferred between chlorophyll 
and other molecules 
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Marcus Joins the Nobel Ranks 


Professor Rudolph A. Marcus receiving the 1992 Nobel Prize in chemistry from King Carl Gustaf of Sweden for his pioneering 
research which has led to a better understanding of what os when atoms join to make new molecules. He is the Arthur 
A. Noyes Professor of Chemistry at the California Institute of Technology and for 40 years has been a principal investigator for 
the Office of Naval Research. Under the sponsorship of ONR, Marcus published between 1956 and 1964 his theory of elec- 
tron-transfer in a series of now classic papers in the scientific literature. The concepts and equations, which Marcus devel- 
oped, he has since expanded and refined into what is known today as the Marcus Theory. 


Currently with the support of ONR, Marcus and his research associates are delving into new aspects of electron transfer. A 
current focus, says Marcus, is on reactions in which electron transfer occurs over distances of a billionth of a meter or so -- 
quite long in the molecular world — in cytochrome C and other protiens important in photosynthesis and cellular respiration 
Such studies should reveal how energy flows through these large molecules and pherhaps point toward possible spinofts, 
such as improved devices for tapping solar energy. (Photo is reproduced through the courtesy of Tobbe Gustavsson, Report- 
agebild, Stockholm, Sweden) 
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In photosynthesis plants use light to form energy-rich compounds from water and carbon dioxide from the air. The 
decisive reactions occur in a reaction center in the cell. The incoming light is caught in an antenna system and finds 
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Introduction 


Peter Schmidt, Guest Editor 
Chemistry and Materials Divisions 
Office of Naval Research 


Being invited to participate in the activities during the 
week of the award of the Nobel Prizes is both an honor and a 
thrill that is hard to put exactly in words. Anticipation and 
excitement are words that capture, in part, feelings I had before 
the week of events. The week itself had a dream-like quality 
that, in retrospect, makes it difficult to believe I actually was 
in Stockholm. The thrill of being there was greatly enhanced 
by the fact that someone | have known for more than twenty 
years, Professor R. A. Marcus, was to receive the Nobel Prize 
in Chemistry. It is therefore an equally great honor for me to 
act as a guest editor of this edition of Naval Research Reviews 
that honors Rudy’s accomplishment by reprinting his Nobel 
address as well as publishing articles that reflect the influence 
of his contribution to electron transfer theory, the work for 
which he received the Prize. It is a proud moment for the Navy 
to acknowledge its contribution to Marcus’s work through the 
support provided by ONR to the original research on which 
the award was based. 

In the following few paragraphs I want to do three things. 
I want first to try to give you an idea of the singular importance 
of Marcus’s research to the general understanding we now 
have of chemical processes. I then want to indicate briefly how 
the following articles in this issue of Naval Research Reviews 
build on and reflect the strength of Marcus’s original contri- 
bution. And, finally, I want to end with a personal note, a short 
chronology, as I experienced it, of the events that took place 
during the week of the award of the Nobel Prizes. 

For more than thirty years now, scientists in a number of 
disciplines that have connections to chemistry — for example, 
chemistry itself, biochemistry, solid state chemistry, and ma- 
terials science — have been using Marcus’s name to refer to the 
only successful theory of electron transfer. Marcus theory has 
come simply to mean the fundamental theory of oxidation and 
reduction processes. Let me take a minute to explain what this 
means. Readers familiar with chemistry will recall that very 
early in the study of the subject, one is introduced to chemical 
reactions that involve the gain or loss of charge in individual 
reactants and products. We have known for some time now 
that the charged particle involved is primarily the electron; 
although in a few cases the proton figures in the charge 
transfer. The electron can be exchanged either between chemi- 
cal species in solution or between species in solution at the 
boundary with a metal, an electrochemical charge transfer. In 
many Cases, it is possible to harness the potential for electron 
transfer stored in chemical species to produce electricity: thus, 
the battery. It may seem hard to believe, but before 1956 there 
was no complete understanding at the molecular level, even at 
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Dr. Peter Schmidt (left ) of the Office of Naval Research and long 
time friend of Dr. Marcus congratulates him at the Nobel Cere- 
monies. Beside Schmidt is his son Nicholas. 





the simplest level of approximation, of how these chemical 
processes take place. 

Marcus’s great contribution built on his understanding of 
chemical reactions in general (he had been involved with O. 
K. Rice in the development of unimolecular reaction rate 
theory in what is now referred to as the RRKM theory — the 
last M being Marcus) and on insight into the way in which 
solvent can play a role in activating the electron transfer. He 
details all of this in his Nobel address and there is no need for 
me to repeat it here. What I want to point out, however, is the 
fact that Marcus’s theory, and his approach to the electron 
transfer theory in particular, was far more general that the title 
electron transfer theory would imply. Imitation — whether 
intentional or not — is surely the finest form of flattery. And, 
Marcus’s fundamental theoretical approach was imitated. 


In 1964 Marcus published several papers in which he 
generalized reaction rate theory and the electron transfer the- 
ory in particular. A crucial point in this generalization was the 
conceptual recognition that in a reacting chemical system, 
which can only be described by a countless number (at least 
Avagadro’s number) of coordinates and momenta, the dynam- 
ics of the system along only one coordinate is important. One 
can imagine a surface in a many-dimensional space to which 
this one important coordinate is perpendicular. The particular 
surface represents the condition for which the initial and final 
state energies of the reacting system are the same. In more 
concise terms, at this special surface, energy is conserved for 
a reaction that proceeds from a reactant to a product state. But, 
any point on this surface is not necessarily the optimal point 
for reaction. Marcus was perhaps the first to recognize that the 
optimal point is found through a saddle point; there are well- 
established methods to find this point. It is interesting to note 
that this method yields a functional form generally referred to 





as “Gaussian” for the dependence of the reaction on various 
energy quantities. The appearance of this functional form led 
Marcus to predict unusual behavior for chemiluminescent 
electron transfer reactions [R. A. Marcus, J. Chem. Phys., 43, 
2654 (1965)]}. The verification of the Gaussian dependence of 
the rate constant on free energy has been verified. 

Several years later, a very nice paper by Glyde appeared 
in the Reviews of Modern Physics [Rev. Mod. Phys., 39, 373 
(1967)] concerning a conundrum of the time — the temperature 
dependence of the self-diffusion coefficient for fluoride in an 
FCC lattice. The analysis is virtually the same as Marcus 
originally provided. And, although Marcus’s work was appar- 
ently not known to Glyde, the omission is understandable 
[Glyde was working in a separate area] and must not be held 
against him. Glyde provides a convincing independent verifi- 
cation of the solidity of Marcus’s pioneered approach. 

There are other areas of investigation of phenomena that 
overlap the Marcus theory. One notable area is that of the 
mobility of the small polaron. This particular work, given a 
beautiful treatment by Holstein in 1959 [T. Holstein, Ann. 
Phys. (NY), 8, 325 (1959)], derives expressions for the mobil- 
ity of a local excess charge in weakly polar lattice that are very 
much the same as the expressions derived by Marcus. It should 
be noted, however, that Marcus used a dielectric response 
method developed by Pekar to handle the effect of a polar 
environment on a transferring electron. It is natural, therefore, 
to expect that theory developing in the areas of polar solids 
and polar liquids might share some similarities. The similari- 
ties are indeed remarkable. 

An additional facet of Marcus’s research that is known 
perhaps to fewer people than it should be has to do with the 
effect of polar solvents on spectral properties. The Marcus 
theory was cast into a different form in 1959-60 by the late V. 
G. Levich and R. R. Dogonadze [V. G. Levich and R. R. 
Dogonadze, Doklady Akad, Nauk, SSSR. 124, 123 (1959); 
Coll. Czech. Chem. Comm., 26, 293 (1961)]. Beginning with 
a Born-Oppenheimer representation of the electron trans- 
fer/solvent system, they were able to use time-dependent 
perturbation theory to derive an expression for the electron 
transfer rate constant. Much of the expression they derived was 
the same as that found by Marcus, and they noted the similari- 
ties. The frequency factor was in fact the same as derived by 
Marcus using the Landau-Zener theory for what is called the 
"non-adiabatic limit”; this is a limit in which there is very weak 
interaction between the reactant and product states of a reac- 
tion. A number of subsequent papers by a variety of authors 
during the next fifteen years built on various aspects of the 
Marcus and Levich-Dogonadze approaches to electron trans- 
fer. But, the interesting aspect of this work was the fact that 
the quantum mechanical treatment to obtain spectral line 
shapes (giving the shape of an absorption band as a function 
of frequency) was formally very similar to the quantum me- 
chanical treatment of the electron transfer reaction. Indeed, the 


similarity between the use of the polar solvent modes in the 
electron transfer reaction and their role in determining spectral 
properties was not lost originally on Marcus. A paper [R. A. 
Marcus, J. Chem. Phys., 24, 979, (1956)] that accompanies his 
Nobel Prize-winning paper of 1956 deals with general polar 
solvent effects that, in the work of others, eventually grew into 
the first account of polar solvent effects on spectra [E. G. 
McRay, J. Phys. Chem., 61, 562 (1957)]. The connection 
between spectral shape functions and chemical reactivity is 
now well established and is much more than a theoretical 
chemistry cottage industry; it is the foundation of many recent 
advances in chemical reaction rate theory. 

Thus, the influence of Marcus’s theoretical contributions 
is great. The theory has survived the test of time. There is firm 
experimental corroboration now that did not exist twenty years 
ago. We are indeed now in possession of a much greater and 
more sophisticated understanding of electron transfer reac- 
tions and chemical reactions in general. 


I turn now to introduce the papers that follow Marcus’s 
Nobel address. The research described by these authors builds 
on or corroborates Marcus theory. We have selected a variety 
to reflect some of the exciting recent developments in electron 
transfer. Although not represented fully here, I want to note 
that one of the most exciting new areas of investigation for the 
electron transfer reaction is in the biological area. Electron 
transfer, it seems obvious to say, has to be an important part of 
many processes in living organisms. Uncovering these proc- 
esses alone has to be a difficult experimental challenge. Un- 
ravelling the electron transfer pathways is being aided directly 
by the unifying conceptual approach provided by Marcus 
theory. That aspects of the theory need to be modified to apply 
to the complexity of a biological system is a change in degree, 
not substance. 


In recent years, drawing on the wealth of formalism 
offered by correlation function theory, it has been possible to 
incorporate Feynman path integral methods into Marcus the- 
ory. At first sight, the proof that Feynman’s methods work may 
not seem to add more than the effort is worth, but this is an 
illusion. Greg Voth provides us with an account that shows 
nicely how the added sophistication of the Feynman path 
integral approach yields much greater insight than would have 
been possible plodding along using methods that many (in- 
cluding me) have used for years. Moreover, the path integral 
approach may indeed bring closer to fruition the hope (long 
held by me) that a computational theory of electron transfer, 
similar to the currently available quantum chemical programs, 
might be available to all. 


Mark Ratner provides a wealth of experience and knowl- 
edge of chemical reactivity. His interest in the basic theory of 
chemical reactivity is long standing. Mark has had a continu- 
ing interest in electron transfer and some his interest focuses 
on the use density functional formalism. In addition, over the 
years he has tackled a number of physically interesting and 
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important issues such as the mobility of charged ions in 
random (often polymeric) networks. Much of that work has 
enhanced our current understanding of the operation of a 
number of new high energy density polymer film battery 
systems. Some of this work has been carried cut with ONR 
participation. 

Professors Weaver and Hupp make contact with experi- 
mental reality. I have known Mike Weaver for a number of 
years now and I am always impressed by his enthusiastic and 
constant probing the limits of applicability of all aspects of the 
electron transfer theory. His, and Joe Hupp’s, relentless inves- 
tigations have provided us with a valuable body of valuable 
information concerning the role of the polar solvent in oxida- 
tion-reduction reactions. They occupy, therefore, a command- 
ing vantage from which to survey the match of theory to 
experiment. 


Finally, Michael Zerner provides us with an account of 
his recent successful work to improve the incorporation of 
polar solvent effects into quantum chemical calculations. One 
of the great difficulties encountered in the past few years in 
accounting for polar solvent effects particularly in calculations 
of biochemical and biophysical importance has been the fact 
that all too often simple dielectric models do not work. In a 
review I wrote about twenty years ago [P. P. Schmidt, Special- 
ist Periodical Reports (The Chemical Society, London), Elec- 
trochemistry, 5, 21-131 (1975)}, I concluded that the polar 
solvent model embodied in the dielectric continuum approxi- 
mation would serve for some time to come. This prediction 
has proved accurate. The seminal contributions to this ap- 
proximation are, one, Marcus’s original recognition that the 
polar continuum could be incorporated properly into theories 
of chemical reactivity and spectroscopy in solution, and, now, 
two, Zerner’s methods to incorporate dielectric continuum 
manipulations properly into quantum chemical calculations. 
There is a thread that begins with Marcus’s work and continues 
to present-day applications. 


To close this introduction, I thought it would be of interest 
to some readers to have a brief account of the activities that 
take place during the week of the Nobel Ceremonies. 

The prizes in physics, chemistry, physiology or medicine, 
and literature and a prize in economics in the name of Nobel 
are awarded in Stockholm on the anniversary of the death of 
Alfred Nobel, December 10. The Peace Prize is awarded in 
Oslo, Norway at the same time. In the time leading up to the 
actual awards, each of the prize winners gives a Nobel address, 
of which Marcus’s is reprinted here. 


Because of the way I scheduled travel, I unfortunately 
missed being able to attend Professor Walcott’s address on 
literature on Monday of the week of the awards. I arrived with 
my family Monday morning in time to get a little rest before 
the reception for the chemistry, physics, and economics prize- 
winners at the Academy of Sciences that evening. This recep- 
tion is one of the first of many that the prize-winners attend 
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and it was an opportunity for each of us to meet the recipients 
and give them our congratulations. For me, it was pure pleas- 
ure to be able to see Rudy receiving the recognition for his 
work that I have admired so long. It was also a fitting recog- 
nition for a fine human being. 

The main events of the next day were the lectures given 
by Professors Marcus and Charpak in chemistry and physics 
and by Professors Fischer and Krebs for physiology and 
medicine. The lectures in physics and chemistry were held at 
the Academy of Sciences and those in physiology at the 
Karolinska Institute. Because of the distance between these 
two institutions, it proved not possible to attend all lectures 
and also attend the Nobelkonsert (concert) at Concert Hall 
Tuesday evening. 

Although it would seem that with only a few lectures 
during each day there would be plenty of time to go sightsee- 
ing, we found this difficult. The bus ride from the Grand Hotel 
to the Academy of Science was sufficiently long that trips back 
and forth — to get lunch or to prepare for the next event — took 
more time than anticipated. Of course, there were ample 
opportunities to converse with other attendees, among whom 
were Harrison Shull, Provost of the Naval Postgraduate 
School, and his wife. 


The main address Wednesday was in Economics by Pro- 
fessor Becker. This was also one of the days for some sight- 
seeing. The award-winners, their families and dignitaries were 
feted at various embassies. We used the time to tour the old 
city of Stockholm. 


Thursday was the jewel of the week because of the award 
ceremony and banquet afterward. It can only be said that these 
two events are magic. The Concert Hall where the awards were 
presented was filled with men and women in formal attire. 
White tie and tails for men and formal gowns for women were 
required. We were given seats with the diplomatic corps with 
a fine view of the stage. This is an event for which those 
entitled to wear them display all their medals. Green, blue or 
crimson ribbons supporting prominent distinctions abounded. 
It was truly a glittering occasion. 

The stage was occupied by the Royal Family, members of 
the Academy of Science and the Karolinska Institute, and the 
winners of the Nobel Prizes. The symphony orchestra, under 
Valdimir Ashkenazy, provided music. Each prize winner was 
introduced to King Carl X VI Gustav and his accomplishments 
noted in Swedish. Following this introduction, the King 
awarded the prize-winner with the gold Nobel medal and 
diploma. Each diploma was individually crafted by an artist to 
capture the essence of the work for which the prize was given. 

If the awarding of the prizes was glittering, the Award 
Banquet that followed was spectacular. It was a splendid 
celebration of human excellence and accomplishment at- 
tended by approximately 1,500 people, some of whom had 
received the Nobel Prize in earlier years. It is almost impossi- 
ble for the common mortal in this century to imagine state 





dinners if one is not part of the diplomatic service or attached 
at a high level to a government office. The Nobel Award 
banquet has its own service that is used only for that occasion. 
The City Hall, in which the banquet is held, literally sparkles. 
Nothing is left to chance. The dinner is a lavish affair with each 
course introduced by music or music and dance. In a touch that 
I found appealing, Swedish university students serve as ushers 
and guides. They were easily identified by their characteristic 
student caps and sashes with emblems that indicated the par- 
ticular “Student Nations” (fraternities or sororities) to which 
they belonged. 

Each of us was assigned a place at the banquet. Our seat 
assignments were listed in a program that contained all of our 
names, including the names of the Royal Family. My wife and 
I found ourselves with friends and the friends and acquain- 
tances of Professor Marcus. Even though the food was delight- 
ful, and memorable, the conversation with those at our table 
was equally enjoyable. And, even though the meal, with its 
many courses, took hours to complete, the time passed all too 
quickly. At the end of the banquet, we were ushered to a hall 
upstairs in which the Ball took place. There were two dance 
bands that entertained well into the night. The students, we 
understand, generally celebrate until dawn. Being somewhat 
past that age, and relying at that hour, on bus transportation 
back to the Grand Hotel, most of the rest of us reluctantly 
ended the evening with a trip home around midnight. 

For the prize winners, the ceremonies did not end with the 
banquet. They then moved on throughout Sweden to present 
their Nobel Addresses at several universities. For us, it was a 
time to wonder if it all really happened, visit one more sight, 
and to see friends and colleagues before boarding the plane for 
the trip home. 
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Electron Transfer: Simple 
Model, Successful Theory 


Peter Schmidt 
Chemistry and Materials Divisions 
Office of Naval Research 


Although many people have a basic knowledge of chem- 
istry, the somewhat sanitized-sounding title “Electron Transfer 
Reaction” masks a well-known process. It is the familiar 
oxidation-reduction reaction. While there are many examples 
of homogeneous oxidation-reduction reactions that take place 
in solution, electrochemical reactions that take place at solid 
surfaces also belong. Oxidation and reduction involve the loss 
and gain of charge, generally electronic charge, in chemical 
species that undergo reaction. An interesting feature of some 
of these reactions is the fact that charge can transfer without 
the formation or destruction of distinct chemical bonds. This 
fact made it possible for Marcus to construct a theory for the 
reaction that has since been verified. My purpose is to discuss 
the model of these reactions in a simple costume, in much the 
same way as Marcus originally did. 

There has been an increasing effort during the past decade 
or so to extend the electron transfer theory; yet there has been 
no modification of the foundation of the theory. It is notewor- 
thy that with the use of the clearly understandabie arguments 
based in electrostatics and statistical mechanics, Marcus was 
able to provide a clear understanding of how these reactions 
occur’. His original equations are accurate in the adiabatic 
limit for the model of the solvent as a polar continuum. There 
has been no subsequent contribution to the body of theory to 
contravene this discovery. Indeed, when the theory is devel- 
oped beginning with the quantum mechanical expression for 
the transition probability”, the same expression emerges that 
Marcus derived making additional use of the Landau-Zener’ 
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theory to describe weak overlap, or non-adiabatic, reactions. 
It is useful, therefore, to examine the model of the electron 
transfer and its interpretation that leads to the Marcus theory. 
A system in which simultaneous oxidation and reduction 
occurs is the following. Two charged, spherical species, A™* 
and B™, dissolved in a polar solvent come into proximity and 
exchange charge according to the following reaction: 


A™ + B™ — A&* + BOT, 


The available electronic charge on one species, A, migrates to 
the second species, B. In the process, energy contained in the 
whole system must localize in the region of the individual 
molecular reacting species to provide the driving force to 
activate the transfer. The transition state is higher in energy 
than either the initial or final state. Thus the energy needed to 
activate the reaction passes back to the surroundings after the 
system moves through the transition to the final state. As long 
as the system is free of external sources of additional energy 
— unconnected to a battery, for example — there is no net gain 
or loss of energy over time. But, energy within the system can 
redistribute itself, and is therefore available to activate chemi- 
cal reactions. If, for example, one of the reactants, B, is 
replaced by a metal surface, M, the reaction with an individual 
state of the metal surface can be written in essentially the same 
form: 


A™ +M—>A®™)* 4 M* 


Although this short article reviews the basic notions of the Marcus theory, it is not intended as any kind of exhaustive review. As a result, I do not provide a 


large list of references. References to some of the earlier reviews are given. The articles to follow in this issue of the Naval Research Reviews will contain 
additional references to current work. 





A™ + M** > A®)* + Me 


but now, as indicated, there are many of these reactions to 
consider, k of them. The index & also labels a level within a 
band of energy levels at the metal surface. An electrochemical 
electron transfer is therefore a sum of all the possible surface 
reactions. Because the energy of a metal surface state can be 
changed with the imposition of an external potential, the 
expressions for the electrochemical electron transfer rate con- 
stant will differ from the homogeneous case by terms that 
reflect the effect of the external potential. 

As is customary with the thermodynamic analyses of 
chemical systems, we imagine the system composed of two 
main parts: (1) the individual reactants, the A and B-species, 
and (2) the solvent. The model is further refined by considering 
the fact that, because of the great difference between the mass 
of the electron and the masses of the other atomic and molecu- 
lar species in the system, the motion of the electron to a first 
approximation can be referred to the remainder of the heavier 
species held at rest. In molecular quantum mechanics this 
assumption is referred to as the Born-Oppenheimer adiabatic 
approximation. It is essential here*”. 

With tis subdivision of the system into local reactants 
and the larger surroundings, it is now possible to account for 
a reaction as follows’. Given a solution in which oxidation and 
reduction takes place, the migrating electronic charge marks 
this change and moves in response to the fluctuations of all the 
other atomic and molecular species. Even in a system that is 
in a state of thermodynamic equilibrium, with overall total 
energy a constant quantity, local energy within the system is 
continually redistributed. It is therefore possiole that sufficient 
energy can collect at the point of reaction for an instant in time 
to activate the process. 

It is simpler first to consider the case in which two species 
A and B exchange charge and both species are chemically 
different: for example, the reaction 


Fe** + Ce* — Fe** + Ce**. 


The (thermodynamic) energy of the system described by 
the reactants on the left of the equation is clearly different from 
the products on the right. When an electron migrates, the 
transfer takes place at a point in the many-dimensional phase 
space where the energy of the initial matches the energy of the 
final state. There are many places where this condition is 
satisfied but are not transition states. One of the states is the 
transition state and eventually will be identified*’. 

Marcus has shown that the free energy of charging two 
ions, labelled ] and 2, in solution is'*? 


aga(t._1)(441 , 462 , Serer 
© op €, 2a, 2a R 


The dielectric constant t. is associated with the optical modes 
and is generally equal to the square of the index of refraction. 
The second dielectric constant, €,, is the static constant for the 
solvent — for water it is about 78. The charge difference Ae is 
e* - e where one has taken the system from a charge state e to 
a new State e* and back. 

The optimum free energy of activation, AG", is deter- 
mined through the requirement that (1) the energies of the 
initial and final states be equal at the transition state, and (2), 
that the point at which the energies are equal be the lowest 
possible transition state energy. This condition is ensured 
through the requirement that the variation of the free energy 
of formation for the overall reaction be zero: 5AG, = 0. 
Together with the requirement for the conservation of charge, 
one can take the variation of AG to find the free energy of 
formation of the transition state’: 


(Ag + AGy )? 
Ao 


l 1 1 1 1 

=Ae|——— |i—+——-— 

ro Cop Ep 2a, 2a, R 
and now Ae is the quantity of charge transferred in the reaction. 
AG' is the famous Marcus expression for the free energy of 


activation of an electron transfer reaction. The expression for 
the rate constant is then! 


AG™ 


k=Zexp(-AG'/RT ) 


and Z is the “collision” number. The unusual feature, for the 
time, was the quadratic dependence of the free energy of 
activation on the square of the free energy of reaction. It is just 
this dependence that has been demonstrated to hold for reac- 
tions in the extremely endo- and exothermal regions, the 
so-called abnormal reactions. 

It seems evident that the simple model of a simple electron 
transfer reactions is indeed too simple. The modes of the 
dielectric continuum alone do not account for all of the process 
of activation. Nevertheless, consideration of classical polari- 
zation does account for much of the activation even though 
more sophistication seems warranted. Marcus, in his early 
electron transfer work, improved the original theory by ac- 
counting for discrete vibrations of solvent and ligand in the 
immediate vicinity of the reactive species''. This adds a term 
for the inner sphere contribution, in effect a Stokes energy. 
There are in addition many systems within which electronic 
charge transfers but which are not well represented as dielec- 
tric systems. The vast collection of biological electron transfer 
systems fall largely into this class'*. In these systems, charge 
migrates over great distances. It is generally believed that 
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superexchange mechanisms are involved. The migratory elec- 
tron makes use of intervening molecular species as bridges 
between donor and acceptor. The polarization of the system 
still figures in the calculations, but must be handled carefully. 

One of the more interesting developments over the past 
two decades has been, I believe, the realization that the process 
of electron transfer resembles to a significant degree radiation- 
less transitions in molecules suspended in condensed phases". 
The understanding that one could cast the rate constant into a 
dependence on a time-dependent correlation function opened 
the analysis of, first, the electron transfer reaction, and then 
the remainder of reactions to useful quantum mechanical 
treatment. It is also interesting to note that many of the indi- 
viduals involved in the development of radiationless transition 
theory branched into the development of the electron transfer 
reaction as a natural outlet'*. This process continues. One of 
the main difficulties in both the radiationless transition theory 
in spectroscopy and electron transfer theory has to do with the 
problem of bond break that accompanies charge transfer. Bond 
break in general is a process that begs to go beyond the 
harmonic limit of analysis. The harmonic limit is generally the 
savior of the electron transfer and radiationless transition 
theory. Some attempts have been made to try to incorporate 
adequate account of anharmonic effects. This has not been 
done to my satisfaction, but the continuing effort is well worth 
the price. What continues to be noteworthy, however, is the 
simple fact that the harmonic limit of analysis of the contribu- 
tion of molecular vibrational modes to the activation of any 
process seems to give us the majority of the energy we expect. 
Were this not the case, Marcus’s analysis never would have 
been possible. With the polaron type of analysis of the polar 
mode representation — as was originally done by Levich and 
Dogonadze” — however, one immediately sees that it is the 
representation of energy states of the polarization that admit 
the harmonic representation for two reasons: (1) the energy of 
the polarization field depends on | P e P dV, and (2) boson 
quantization. This is of course implicit in Marcus’s original 
analysis’ because it is a basic part of (linear) electrostatics. 
Non-linear polarization effects might play some role, but one 
can anticipate that they would be very small. 

Where are the uncharted regions of the electron transfer 
theory now? I think there are a few. The important develop- 
ment for the near term, I believe, will continue to be with the 
construction of an adequate model and analysis of the contri- 
bution of discrete molecular vibrations to the activation proc- 
ess. For, while these modes may not dominate the activation 
in comparison to the polar continuum, they do nevertheless 
manifest subtle effects that may be very important in our 
efforts to interpret important biological and materials science 
applications. Finally, there has been, and continues to be, a 
desire among many to move well past the continuum repre- 
sentation of the solvent. Efforts have been made’ to formalize 
a theoretical dependence on various better-defined, more rig- 
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orous models of polar solutions, e.g., the mean spherical 
approximation. There need to be reasonable computational 
algorithms that allow one to use effectively Monte Carlo or 
molecular dynamics routines to find rate parameters. Some of 
this effort is discussed in this issue of the Reviews. 

The Marcus theory is elegant in its design, its execution, 
and its apparent simplicity of expression. The evolution of the 
theory will add texture to a strong foundation. Chemistry is 
now richer for the greater understanding of nature that Marcus 
theory has helped us to attain. 
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Looking Back 


My first encounters with McGill University came when I 
was still in a baby carriage. My mother used to wheel me about 
the campus when we lived in that neighborhood and, as she 
recounted years later, she would tell me that I would go to 
McGill. There was some precedent for my going there, since 
two of my father’s brothers received their M.D.’s at McGill. 

I have always loved going to school. Since neither of my 
parents had a higher education, my academic “idols” were 
these two paternal uncles and one of their uncles, my great-un- 
cle, Henrik Steen (ne Markus). My admiration for him, living 
in faraway Sweden, was not because of a dr.telol. (which he 
received from the University of Uppsala in 1915) nor because 
of the many books he wrote — I knew nothing of that — but 
rather because he was reputed to speak 13 languages. I learned 
decades later that the number was only 9! Growing up, mostly in 
Montreal, I was an only child of loving parents. I admired my 
father’s athletic prowess — he excelled in several sports — and my 
mother’s expressive singing and piano playing. 

My interest in the sciences started with mathematics in 
the very beginning, and later with chemistry in early high 
school and the proverbial home chemistry set. My education 


at Baron Byng High School was excellent, with dedicated 
masters (boys and girls were separate). I spent the next years 
at McGill University, for both undergraduate and, as was the 
custom of the time, graduate study. Our graduate supervisor, 
Cari A. Winkler, specialized in rates of chemical reactions. He 
himself had received his Ph.D. as a student of Cyril Hinshel- 
wood at Oxford. Hinshelwood was later the recipient of the 
Nobel Prize for his work on chemical kinetics. Winkler 
brought to his laboratory an enthusiastic joyousness in re- 
search and was much loved by his students. 

During my McGill years, I took a number of math courses, 
more than other students in chemistry. Upon receiving a Ph.D. 
from McGill University in 1946, I joined the new post-doctoral 
program at the National Research Council of Canada in 
Ottawa. This program at NRC later became famous but at the 
time it was still in its infancy and our titles were Junior 
Research Officers. The photochemistry group was headed by 
E.W.R. Steacie, an international figure in the study of free- 
radical reactions and a major force in the development of the 
basic research program at NRC. I benefitted from the quality 
of his research on gas phase reaction rates. Like my research 
on chemical reaction rates in solution at McGill (kinetics of 
nitration), it was experimental in nature. There were no theo- 
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retical chemists in Canada at the time, and as students I don’t 
think we ever considered how or where theories were con- 
ceived. 

About 1948 a fellow post-doctoral at NRC, Walter Trost, 
and I formed a two-man seminar to study theoretical papers 
related to our experimental work. This adventure led me to 
explore the possibility of going on a second post-doctoral, but 
in theoretical work, which seemed like a radical step at the 
time. I had a tendency to break the glass vacuum apparatus, 
due to a still present impetuous haste, with time-consuming 
consequences. Nevertheless, the realization that breaking a 
pencil point would have far less disastrous consequences 
played little or no role, I believe, in this decision to explore 
theory! 

I applied in 1948 to six well-known theoreticians in the 
U.S. for a postdoctoral research fellowship. The possibility 
that one of them might take on an untested applicant, an 
applicant hardly qualified for theoretical research, was prob- 
ably too much to hope for. Oscar K. Rice at the University of 
North Carolina alone responded favorably, subject to the suc- 
cess of an application he would make to the Office of Naval 
Research for this purpose. It was, and in February 1949 I took 
the train south, heading for the University of North Carolina 
in Chapel Hill. I was impressed on arrival there by the red clay, 
the sandy walks, and the graciousness of the people. 

After that, I never looked back. Being exposed to theory, 
stimulated by a basic love of concepts and mathematics, was 
a marvelous experience. During the first three months I read 
everything I could lay my hands on regarding reaction rate 
theory, including Marcelin’s classic 1915 theory which came 
within one small step of the Transition State Theory of 1935. 
I read numerous theoretical papers in German, a primary 
language for the “chemical dynamics” field in the 1920s and 
1930s, attended my first formal course in quantum mechanics, 
given by Nathan Rosen in the Physics Department, and was 
guided by Oscar in a two-man weekly seminar in which I 
described a paper I had read and he pointed out assumptions 
in it that I had overlooked. My life as a working theorist began 
three months after this preliminary study and background 
reading, when Oscar gently nudged me toward working on a 
particular problem. 

Fortunately for me, Oscar’s gamble paid off. Some three 
months later, I had formulated a particular case of what was 
later entitled by B. Seymour Rabinovitch, RRKM theory 
(“Rice-Ramsperger-Kassel-Marcus’”). In it, I blended statisti- 
cal ideas from the RRK theory of the 1920s with those of the 
transition state theory of the mid-1930s. The work was pub- 
lished in 1951. In 1952 I wrote the generalization of it for other 
reactions. In addition, six months after arrival in Chapel Hill, 
I was also blessed by marriage to Laura Hearne, an attractive 
graduate student in sociology at UNC. She is here with me at 
this ceremony. Our three sons, Alan, Kenneth and Raymond, 
and two daughters-in-law are also present today. 
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In 1951, I attempted to secure a faculty position. This 
effort met with little success (35 letters did not yield 35 no’s, 
since not everyone replied!). Very fortunately, that spring I met 
Dean Raymond Kirk of the Polytechnic Institute of Brooklyn 
at an American Chemical Society meeting in Cleveland, which 
I was attending primarily to seek a faculty position. This 
meeting with Dean Kirk, so vital for my subsequent career, 
was arranged by Seymour Yolles, a graduate student at UNC 
in a course I taught during Rice’s illness. Seymour had been a 
student at Brooklyn Poly and learned, upon accidentally en- 
countering Dr. Kirk, that Kirk was seeking new faculty. After 
a subsequent interview at Brooklyn Poly, I was hired, and life 
as a fully independent researcher began. 

I undertook an experimental research program on both gas 
phase and solution reaction rates, wrote the 1952 RRKM 
papers, and wondered what to do next in theoretical research. 
I felt at the time that it was pointless to continue with RRKM 
since few experimental data were available. Some of our 
experiments were intended to produce more. 

After some minor pieces of theoretical study that I worked 
on, a student in my statistical mechanics class brought to my 
attention a problem in polyelectrolytes. Reading everything I 
could about electrostatics, I wrote two papers on that topic in 
1954/55. This electrostatics background made me fully ready 
in 1955 to treat a problem I had just read about on electron 
transfers. I comment on this next period on electron transfer 
research in my Nobel Lecture. About 1960, it became clear 
that it was best for me to bring the experimental part of my 
research program to a close — there was too much to do on the 
theoretical aspects — and I began the process of winding down 
the experiments. I spent a year and a half during 1960-61 at 
the Courant Mathematical Institute at New York University, 
auditing many courses which were, in part, beyond me, but 
which were, nevertheless, highly instructive. 

In 1964, I joined the faculty of the University of Illinois 
in Urbana-Champaign and I never undertook any further ex- 
periments there. At Illinois, my interests in electron transfer 
continued, together with interests in other aspects of reaction 
dynamics, including designing “natural collision coordi- 
nates”, learning about action-angle variables, introducing the 
latter into molecular collisions, reaction dynamics, and later 
into semiclassical theories of collisions and of bound states, 
and spending much of my free time in the astronomy library 
learning more about classical mechanics, celestial mechanics, 
quasiperiodic motion, and chaos. I spent the academic year of 
1975-76 in Europe, first as Visiting Professor at the University 
of Oxford and later as a Humboldt Awardee at the Technical 
University of Munich, where I was first exposed to the prob- 
lem of electron transfer in photosynthesis. 

In 1978, I accepted an offer from the California Institute 
of Technology to come there as the Arthur Amos Noyes 
Professor of Chemistry. My semiclassical interlude of 1970-80 
was intellectually a very stimulating one, but it involved for 





me less interaction with experiments than had my earlier work 
on unimolecular reaction rates or on electron transfers. Ac- 
cording, prompted by the extensive experimental work of my 
colleagues at Caltech in these fields of unimolecular reactions, 
intramolecular dynamics and of electron transfer processes, as 
well as by the rapidly growing experimental work in both 
broad areas world-wide, I turned once again to those particular 
topics and to the many new types of studies that were being 
made. Their scope and challenge continues to grow to this day 
in both fields. Life would be indeed easier if the experimen- 
talists would only pause for a little while! 

There was a time when I had wondered about how much 
time and energy had been lost doing experiments during most 
of my stay at Brooklyn Poly — experiments on gas phase 
reactions, flash photolysis, isotopic exchange electron trans- 
fer, bipolar electrolytes, nitration, and photoelectrochemistry, 
among other — and during all of my stay at NRC and at McGill. 
In retrospect, I realized that this experimental background 
heavily flavored my attitude and interests in theoretical re- 
search. In the latter I drew, in most but not all cases, upon 
experimental! findings or puzzles for theoretical problems to 
study. The growth of experiments in these fields has served as 
a continually rejuvenating influence. This interaction of ex- 
periment and theory, each stimulating the other, has been and 
continues to be one of the joys of my experience. 

Honors received for the theoretical work include the 
Irving Langmuir and the Peter Debye Awards of the American 
Chemical Society (1978, 1988), the Willard Gibbs, Theodore 
William Richards, and Pauling medals, and the Remsen and 





Figure 1. 


Examples of topics in the electron transfer field (Marcus and 
Siddarth, ref. 2) 
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Edgar Fahs Smith Awards, from various sections of the ACS, 
(1988, 1990, 1991, 1991, 1991), the Robinson and the Centen- 
ary Medals of the Faraday Division of the Royal Society of 
Chemistry (1982, 1988), Columbia University’s Chandler 
Medal (1983) and Ohio State’s William Lloyd Evans Award 
(1990), a Professorial Fellowship at University College, Ox- 
ford (1975 to 1976) and a Visiting Professorship in Theoretical 
Chemistry at Oxford during that period, the Wolf Prize in 
Chemistry (1985), the National Medal of Science (1989), the 
Hirschfelder Prize in Theoretical Chemistry (1993), election 
to the National Academy of Sciences (1970), the American 
Academy of Arts and Sciences (1973), the American Philo- 
sophical Society (1990), honorary membership in the Royal 
Society of Chemistry (1991), and foreign membership in the 
Royal Society (London) (1987) and in the Royal Society of 
Canada (1993). Honorary degrees were conferred by the Uni- 
versity of Chicago and by Goteborg, Polytechnic, McGill, and 
Queen’s Universities and by the University of New Brunswick 
(1983, 1986, 1987, 1988, 1993, 1993). A commemorative issue 
of the Journal of Physical Chemistry was published in 1986. 


Electron Transfer Experiments 
Since The Late 1940s 


Since the late 1940s, the field of electron transfer proc- 
esses has grown enormously, both in chemistry and biology. 
The development of the field, experimentally and theoreti- 
cally, as well as its relation to the study of other kinds of 
chemical reactions, represents to us an intriguing history, one 
in which many threads have been brought together. In this 
lecture, some history, recent trends, and my own involvement 
in this research are described. 

The early experiments in the electron transfer field were 
on “isotopic exchange reactions” (self-exchange reactions) 
and, later, “cross reactions.” These experiments reflected two 
principal influences. One of these was the availability after the 
Second World War of many radioactive isotopes, which per- 
mitied the study of a large number of isotopic exchange 
electron transfer reactions, such as 


Fe* + Fe>* — Fe** + Fe?” (1) 
and 


Ce** + Ce 4+ C™+Ce™, (2) 
in aqueous solution, where the asterisk denotes a radioactive 
isotope. 

There is a two-fold simplicity in typical self-exchange 
electron transfer reactions (so-called since other methods be- 
sides isotopic exchange were later used to study some of 
them): (1) the reaction products are identical with the reac- 
tants, thus eliminating one factor which usually influences the 
rate of a chemical reaction in a major way, namely the relative 
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thermodynamic stability of the reactants and products; and (2) 
no chemical bonds are broken or formed in simple electron 
transfer reactions. Indeed, these self-exchange reactions rep- 
resent, for these combined reasons, the simplest class of reac- 
tions in chemistry. Observations stemming directly from this 
simplicity were to have major consequences, not only for the 
electron transfer field but also, to a lesser extent, for the study 
of other kinds of chemical reactions as well (cf Shaik et al, ref. 
2). 

Asecond factor in the growth of the electron transfer field 
was the introduction of new instrumentation, which permitted 
the study of the rates of rapid chemical reactions. Electron 
transfers are frequently rather fast, compared with many reac- 
tions which undergo, instead, a breaking of chemical bonds 
and a forming of new ones. Accordingly, the study of a large 
body of fast electron transfer reactions became accessible with 
the introduction of this instrumentation. One example of the 
latter was the stopped-flow apparatus, pioneered for inorganic 
electron transfer reactions by N. Sutin. It permitted the study 
of bimolecular reactions in solution in the millisecond time 
scale (a fast time scale at the time). Such studies led to the 
investigation of what has been termed electron transfer “cross 
reactions,” i.e., electron transfer reactions between two differ- 
ent redox systems, as in 


Fe** + Ce** — Fe** + Ce**, (3) 


which supplemented the earlier studies of the self-exchange 
electron transfer reactions. A comparative study of these two 





Figure 2. 


Potential energy contours for reaction (4), AB + C + A + BC, 
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types of reaction, self-exchange and cross-reactions, stimu- 
lated by theory, was also later to have major consequences for 
the field and, indeed, for other areas. 

Again, in the field of electrochemistry, the new post-war 
instrumentation in chemical laboratories led to methods which 
permitted the study of fast electron transfer reactions at metal 
electrodes. Prior to the late 1940s only relatively slow electro- 
chemical reactions, such as the discharge of an H°O* ion at an 
electrode to form H?, had been investigated extensively. They 
involved the breaking of chemical bonds and the forming of 
new ones. 

Numerous electron transfer studies have now also been 
made in other areas, some depicted in Figure 1. Some of these 
investigations were made possible by a newer technology, 
lasers particularly, and now include studies in the picosecond 
and subpicosecond time regimes. Just recently, (non-laser) 
nanometer-sized electrodes have been introduced to study 
electrochemical processes that are still faster than those hith- 
erto investigated. Still other recent investigations, important 
for testing aspects of the electron transfer theory at electrodes, 
involve the new use of an intervening ordered adsorbed mono- 
layer of long chain organic compounds on the electrode to 
facilitate the study of various effects, such as varying the 
metal-solution potential difference on the electrochemical 
electron transfer rate. 

In some studies of electron transfer reactions in solution 
there has also been a skillful blending of these measurements 
of chemical reaction rates with various organic or inorganic 
synthetic methods, as well as with site-directed mutagenesis, 
to obtain still further hitherto unavailable information. The use 
of chemically modified proteins to study the distance depend- 
ence of electron transfer, notably by Gray and coworkers, has 
opened a whole new field of activity. 

The interaction of theory and experiment in these many 
electron transfer fields has been particularly extensive and 
exciting, and each has stimulated the other. The present lecture 
addresses the underlying theory and this interaction. 


The Early Experience 


My own involvement in the electron transfer field began 
in a rather circuitous way. In an accompanying biographical 
note I have commented on my earlier background, which was 
in experimental measurements of reaction rates as a chemistry 
graduate student at McGill University (1943-46) and as a 
post-doctoral associate at the National Research Council of 
Canada (1946-49). A subsequent post-doctoral study at the 
University of North Carolina (1949-51) on the theory of 
reaction rates resulted in what is now known in the literature 
as RRKM theory (Rice, Ramsperger, Kassel, Marcus). 

This unimolecular reaction field reflects another long and 
extensive interaction between theory and experiment. RRKM 
theory enjoys widespread use and is now usually referred to 





in the literature only by its acronym (or by the texts written 
about it, ref. 4), instead of by citation of the original articles. 

After the theoretical post-doctoral I joined the faculty of 
the Polytechnic Institute of Brooklyn in 1951 and wondered 
what theoretical research to do next after writing the RRKM 
papers (1951-52). I remember vividly how a friend of mine, a 
colleague at Brooklyn Poly, Frank Collins, came down to my 
office every day with a new idea on the liquid state transport 
theory which he was developing, while I, for theoretical re- 
search, had none. Perhaps this gap in not doing anything 
immediately in the field of theory was, in retrospect, fortunate: 
In not continuing with the study of the theory of unimolecular 
reactions, for which there were too few legitimate experimen- 
tal data at the time to make the subject one of continued 
interest, I was open for investigating quite different problems 
in other areas. I did, however, begin a program of experimental 
Studies in gas phase reactions, prompted by my earlier studies 
at NRC and by the RRKM work. 

In the biographical note I have also recalled how a student 
in my Statistical mechanics class in this period (Abe Koltliar) 
asked me about a particular problem in polyelectrolytes. It led 
to my writing two papers on the subject (1954-55), one of 
which required a considerable expansion in my background in 
electrostatics, so as to analyze different methods for calculat- 
ing the free energy of these systems: In polyelectrolyte mole- 
cules, it may be recalled, the ionic charges along the organic 
or inorganic molecular backbone interact with each other and 





Figure 3. 


Typical nuclear configurations for reactants, products, and 
surrounding solvent molecules in reaction (1). The longer 
M-OHz2 bond length in the +2 state is indicated schematically 
by the larger ionic radius. (Sutin, ref. 2) 
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with the solvent. In the process I read the relevant parts of the 
texts that were readily available to me on electrostatics (Cal- 
tech’s Mason and Weaver’s was later to be particularly help- 
ful!). When shortly thereafter I encountered some papers on 
electron transfer, a field entirely new to me, I was reasonably 
well prepared for treating the problems which lay ahead. 


Developing An Electron 
Transfer Theory 


Introduction 


My first contact with electron transfers came in 1955 as a 
result of chancing upon a 1952 symposium issue on the subject 
in the Journal of Physical Chemistry. An article by Bill Libby 
caught my eye a use of the Franck-Condon principle to 
explain some experimental results, namely, why some isotopic 
exchange reactions which involve electron transfer between 
pairs of small cations in aqueous solution, such as reaction (1), 
are relatively slow, whereas electron transfers involving larger 
ions, such as Fe(CN),* — Fe(CN),* and MnO,- -MnO,”, are 
relatively fast. 

Libby explained this observation in terms of the Franck- 
Condon principle, as discussed below. The principle was used 
extensively in the field of spectroscopy for interpreting spectra 
for the excitation of the molecular electronic-vibrational quan- 
tum states. An application of that principle to chemical reac- 
tion rates was novel and caught my attention. In that paper 
Libby gave a “back-of-the-envelope” calculation of the result- 
ing solvation energy barrier which slowed the reaction. How- 
ever, I felt instinctively that even though the idea — that 
somehow the Franck-Condon principle was involved — 
seemed strikingly right, the calculation itself was incorrect. 
The next month of study of the problem was, for me, an 
especially busy one. To place the topic in some perspective | 
first digress and describe the type of theory that was used for 
other types of chemical reaction rates at the time and continues 
to be useful today. 


Reaction Rate Theory 


Chemical reactions are often described in terms of the 
motion of the atoms of the reactants on a potential energy 
surface. This potential energy surface is really the electronic 
energy of the entire system, plotted versus the positions of all 
the atoms. A very common example is the transfer of an atom 
or a group B from AB to form BC 


AB+C—>A+BC. (4) 
An example of reaction (4) is the transfer of an H, such as 
in IH + Br + 1 + HBr, or the transfer of a CH, group from one 
aromatic sulfonate to another. To aid in visualizing the motion 


of the atoms in this reaction, this potential energy function is 
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Figure 4. 


Profile of potential energy surfaces for reactants plus environ- 
ment, R, and for products plus environment, P. Solid curves 
schematic. Dashed curves: schematic but slightly more real- 
istic. The typical splitting at the intersection of U; and Up is not 
shown in the Figure (Marcus and Siddarth, ref. 2) 
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frequently plotted as constant energy contours in a space 
whose axes are chosen to be two important relative coordi- 
nates such as, in reaction (4), a scaled AB bond length and a 


scaled distance from the center of mass of AB toC, as in Figure 
> 


A point representing this reacting system begins its trajec- 
tory in the lower right region of the figure in a valley in this 
plot of contours, the “valley of the reactants.” When the system 
has enough energy, appropriately distributed between the vari- 
ous motions, it can cross the “mountain pass” (saddle-point 
region) separating the initial valley from the products’ valley 
in the upper left, and so form the reaction products. There is a 
line in the figure, XY, analogous to the “continental divide” in 
the Rocky Mountains in the U.S., which separates systems 
which could spontaneously flow into the reactants’ valley from 
those which could flow into the products’ one. In chemists’ 
terminology this line represents the “transition state” of the 
reaction. 

In transition state theory a quasi-equilibrium between the 
transition state and the reactant is frequently postulated, and 
the reaction rate is then calculated using equilibrium statistical 
mechanics. A fundamental dynamical basis, which replaces 
this apparently ad hoc but common assumption of transition 
state theory and which is perhaps not as well known in the 
chemical literature as it deserves to be, was given many years 
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ago by the physicist and one-time chemical engineer, Eugene 
Wigner (1938). He used a classical mechanical description of 
the reacting system in the many-dimensional space (of coor- 
dinates and momenta). Wigner pointed out that the quasi-equi- 
librium would follow as a dynamical consequence, if each 
trajectory of a moving point representing the reacting system 
in this many-dimensional space did not recross the transition 
state (and if the distribution of the reactants in the reactants’ 
region were a Boltzmann one). In recent times, the examina- 
tion of this recrossing has been a common one in classical 
mechanical trajectory studies of chemical reactions. Usually, 
recrossings are relatively minor, except in nonadiabatic reac- 
tions, where they are readily treated (cf discussion, later). 

In practice, transition state theory is generalized, so as to 
include as many coordinates as are needed to describe the 
reacting system. Further, when the system can “tunnel” quan- 
tum mechanically through the potential energy barrier (the 
“pass”) separating the two valleys, as for example frequently 
happens at low energies in H-transfer reactions, the method of 
treating the passage across the transition state region needs, 
and has received, refinement. (The principal problem encoun- 
tered here has been the lack of “dynamical separability” of the 
various motions in the transition state region.) 





Figure 5. 


Free energy of reactants plus environment vs. the reaction 
coordinate Q (R curve), and free energy of products plus 
environment vs. reaction coordinate Q (P curve). The three 
vertical lines on the abscissa denote, from left to right, the 
value for the reactants, for the transition state, and for the 
products. (Marcus and Siddarth, ref. 2) 
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Electron Transfer Theory. Formulation 


In contrast to the above picture, we have already noted 
that in simple electron transfer reactions no chemical bonds 
are broken or formed and so a somewhat different picture of 
the reaction is needed for the electron transfer reaction. 

In his 1952 symposium paper, Libby noted that when an 
electron is transferred from one reacting ion or molecule to 
another, the two new molecules or ions formed are in the 
wrong environment of the solvent molecules, since the nuclei 
do not have time to move during the rapid electron jump: in 
reaction (1) a Fe”* ion would be formed in some configuration 
of the many nearby dipolar solvent molecules that was appro- 
priate to the original Fe** ion. Analogous remarks apply to the 
newly formed Fe** ion in the reaction. On the other hand, in 
reactions of “complex ions,” such as those in the Fe(CN),° 
—-Fe(CN),* and MnO,- —Mn0O,” self-exchange reactions, the 
two reactants are larger, and so the change of electric field in 
the vicinity of each ion, upon electron transfer, would be 
smaller. The original solvent environment would therefore be 
less foreign to the newly formed charges, and so the energy 
barrier to reaction would be less. In this way Libby explained 
the faster self-exchange electron transfer rate for these com- 
plex ions. Further confirmation was noted in the ensuing 


discussion in the symposium: the self-exchange Co(NH,),"° 


-Co(NH,),”* reaction is very slow, and it was pointed out that 


there was a large difference in the equilibrium CoN bond 
lengths in the 3+ and the 2+ ions, and so each ion would be 
formed in a very “foreign” configuration of the vibrational 
coordinates, even though the ions are “complex ions.” 

After studying Libby’s paper and the symposium discus- 
sion, I realized that what troubled me in this picture for 
reactions occurring in the dark was that energy was not con- 
served: the ions would be formed in the wrong high-energy 
environment, but the only way such a non-energy conserving 
event could happen would be by the absorption of light (a 
“vertical transition”), and not in the dark. Libby had percep- 
tively introduced the Franck-Condon principle to chemical 
reactions, but something was missing. 

In the present discussion, as well as in Libby’s treatment, 
it was supposed that the electronic interaction of the reactants 
which causes the electron transfer is relatively weak. That 
view is still the one that seems appropriate today for most of 
these reactions. In this case of weak-electronic interaction, the 
question becomes: how does the reacting system behave in the 
dark so as to satisfy both the Franck-Condon principle and 
energy conservation? I realized that fluctuations had to occur 
in the various nuclear coordinates, such as in the orientation 
coordinates of the individual solvent molecules and indeed in 
any other coordinates whose most probable distribution for the 
products differs from that of the reactants. With such fluctua- 
tions, values of the coordinates could be reached which satisfy 
both the Franck-Condon and energy conservation conditions 
and so permit the electron transfer to occur in the dark. 





Figure 6. 


Plot of the free energy G versus the reaction coordinate q, for 
reactants’ (R) and products’ (P), for three different values of 


AG” the cases ! to Ill indicated in Rigure & (Marcus and 
Siddarth, ref. 2) 
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For a reaction such as reaction (1), an example of an initial 
and final configuration of the solvent molecules is depicted in 
Figure 3. Fluctuations from the original equilibrium ensemble 
of configurations were ultimately needed, prior to the electron 
transfer, and were followed by a relaxation to the equilibrium 
ensemble for the products, after electron transfer. 

The theory then proceeded as follows. The potential en- 
ergy U, of the entire system, reactants plus solvent, is a 
function of the many hundreds of relevant coordinates of the 
system, coordinates which include, among others, the position 
and orientation of the individual solvent molecules (and hence 
of their dipole moments, for example), and the vibrational 
coordinates of the reactants, particularly those in any inner 
coordination shell of the reacting ions. (E.g., the inner coordi- 
nation shell of an ion such as Fe”* or Fe** in water is known 
from EXAFS experiments to contain six water molecules.) No 
longer were there just the two or so important coordinates that 
were dominant in reaction (4). 

Similarly, after the electron transfer, the reacting mole- 
cules have the ionic charges appropriate to the reaction prod- 
ucts, and so the relevant potential energy function U, is that 
for the products plus solvent. These two potential energy 
surfaces will intersect if the electronic coupling which leads 
to electron transfer is neglected. For a system with N coordi- 
nates this intersection occurs on an (N-1) dimensional surface, 
which then constitutes in our approximation the transition state 
of the reaction. The neglected electronic coupling causes a 
well-known splitting of the two surfaces in the vicinity of their 
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intersection. A schematic profile of the two potential energy 
surfaces in the N-dimensional space is given in Figure 4. (The 
splitting is not shown.) 

Due to the effect of the previously neglected electronic 
coupling and the coupling between the electronic motion and 
to the nuclear motion near the intersection surface S, an 
electron transfer can occur at S. In classical terms, the transfer 
at S occurs at fixed positions and momenta of the atoms, and 
so the Franck-Condon principle is satisfied. Since U_ equals 
U_ at S, energy is also conserved. The details of the electron 
transfer depend on the extent of electronic coupling and how 
rapidly the point representing the system in this N- dimen- 
sional space crosses § . (It has been treated, for example, using 
aS an approximation the well-known one-dimensional Lan- 
dau-Zener expression for the transition probability at the near- 
intersection of two potential energy curves.) 

When the splitting caused by the electronic coupling 
between the electron donor and acceptor is large enough at the 
intersection, a system crossing S from the lower surface on the 
reactants’ side of S continues onto the lower surface on the 
products’ side, and so an electron transfer in the dark has then 
occurred. When the coupling is, instead, very weak, (“nonadi- 
abatic reactions”) the probability of successfully reaching the 
lower surface on the products’ side is small and can be calcu- 
lated using quantum mechanical perturbation theory, for ex- 
ample, using Fermi’s “Golden Rule,” an improvement over the 
1-dimensional Landau-Zener treatment. 

Thus, there is some difference and some similarity with a 
more conventional type of reaction such as reaction (4), whose 
potential energy contour plots were depicted in Figure 2. In 
both cases, fluctuations of coordinates are needed to reach the 
transition state, but since so many coordinates can now play a 
significant role in the electron transfer reaction, because of the 
major and relatively abrupt change in charge distribution on 
passing through the transition state region, a rather different 
approach from the conventional one was needed to formulate 
the details of the theory. 


Electron transfer theory. Treatment 


In the initial paper (1956) I formulated the above picture 
of the mechanism of electron transfer and, to make the calcu- 
lation of the reaction rate tractable, treated the solvent as a 
diclectric continuum. In the transition state the position-de- 
pendent diclectric polarization P (0) of the solvent, due to the 
orientation and vibrations of the solvent molecules, was not 
the one in equilibrium with the reactants’ or the products’ ionic 
charges. It represented instead, some macroscopic fluctuation 
from them. The electronic polarization for the solvent mole- 
cules, on the other hand, can rapidly respond to any such 
fluctuations and so is that which is dictated by the reactants’ 
charges and by the instantaneous P(r). 

With these ideas as a basis, what was then needed was a 
method of calculating the electrostatic free energy G of this 
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system with its still unknown polarization function P(r). | 
obtained this free energy G by finding a reversible path for 
reaching this state of the system. Upon then minimizing G, 
subject to the constraint imposed by the Franck-Condon prin- 
ciple (reflected in the electron transfer occurring at the inter- 
section of the two potential energy surfaces), I was able to find 
the unknown P(r) and, hence, to find the G for the transition 
state. That G was then introduced into transition state theory 
and the reaction rate calculated. 

In this research | also read and was influenced by a lovely 
paper by Platzmann and Franck (1952) on the optical absorp- 
tion spectra of halide ions in water and later by work of 
physicists such as Pekar and Frohlich (1954) on the closcly 
related topic of polaron theory. As best as I can recall now, my 
first expressions for G during this month of intense activity 
seemed rather clumsy, but then with some rearrangement a 
simple expression emerged that had the right “feel” to it and 
that I was also able to obtain by a somewhat independent 
argument. The expression also reduced reassuringly to the 
usual one, when the constraint of arbitrary P(r) was removed. 
Obtaining the result for the mechanism and rate of electron 
transfer was indeed one of the most thrilling moments of my 
scientific life. 

The expression for the rate constant k of the reaction is 
given by 


—AG * 


k=Aexp = 
k,l 
where AG’ , in tum, is given by 


AG *= 


aA, Ae 
4 


(1+ (5b) 


Si 


The A in Eq. (Sa) is a term depending on the nature of the 
electron transfer reaction (e.g., bimolecular or intramolecu- 
lar), AG® is the standard free energy of reaction (and equals 
zero for a self-exchange reaction), A is a “reorganization 
term,” composed of solvational (A) and vibrational (A,) com- 
ponents. 


haa) +A, (6) 


In a two-sphere model of the reactants, 4, was expressed 
in terms of the two ionic radii a, and a, (including in the radius 
any inner coordination shell), the center-to-center separation 
distance R of the reactants, the optical (D_ ) and static (D,) 
dielectric constants of the solvent, and the charge transferred 
Ae from one reactant to the other: 





For a bimolecular reaction, work terms, principally elec- 
trostatic, are involved in bringing the reactants together and in 
separating the reaction products, but are omitted from Eq. (5) 
for notational brevity. The expression for the vibrational term 
A, is given by 


1 
hi=5 LK G-OQ (8) 


J 


where Q‘ and Q? are equilibrium values for the jth normal 
mode coordinate Q and k., is a reduced force constant 
2k *k F/(k; +k Pk * being the force constant for the reactants 
and k? being that for the products. (I introduced a “symmetri- 
zation” approximation for the vibrational part of the potential 
energy surface, to obtain this simple form of Eqs. (5) to (8), 
and tested it numerically.) 

In 1957 I published the results of a calculation of the A. 
arising from a stretching vibration in the innermost coordina- 
tion shell of each reactant, (the equation used for A, was given 
in the 1960 paper). An early paper on the purely vibrational 
contribution using chemical bond length coordinates and ne- 
glecting bond-bond correlation had already been published for 
self-exchange reactions by George and Griffiths in 1956. 

I also extended the theory to treat electron transfers at 
electrodes, and distributed it as an Office of Naval Research 
Report in 1957, the equations being published later in a journal 
paper in 1959. I had little prior knowledge of the subject, and 
my work on electrochemical electron transfers was facilitated 
considerably by reading a beautiful and logically written sur- 
vey article of Roger Parsons on the equilibrium electrostatic 
properties of electrified metal-solution interfaces. 

In the 1957 and 1965 work I showed that the electro- 
chemical rate constant was again given by Egs. (5)-(8), but 
with A now having a value appropriate to the different “ge- 
ometry” of the encounter of the participants in the reaction. 
The 1/2a, in Eq. (7) was now absent (there is only one reacting 
ion) and R now denotes twice the distance from the center of 
the reactant’s charge to the electrode (it equals the ion-image 
distance). A term en replaced the AG? in Eq.(Sb), where e is 
the charge transferred between the ion and the electrode, and 
7) is the activation overpotential, namely the metal-solution 
potential difference, relative to the value it would have if the 
rate constants for the forward and reverse reactions were equal. 
These rate constants are equal when the minima of the two G 
curves in Figure 5 have the same height. 

When leni<A, most electrons go into or out of quantum 
States in the metal that are near the Fermi level. However, 
because of the continuum of states in the metal, the inverted 
effect was now predicted to be absent for this process, i.c., the 
counterpart of Eq. (5) is applicable only in the region lenl<A: 
In the case of an intrinsically highly exothermic electron 
transfer reaction at an electrode, the electron can remove the 
immediate “exothermicity” by (if entering) going into a high 





Figure 7. 


Plot of Ln kr vs AG® Points | and Ill are in the normal and 
inverted regions, respectively, while point Il, where In k; is a 


maximum, occurs at -AG° = | (Marcus and Siddarth, ref. 2) 
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unoccupied quantum state of the metal, or (if leaving) depart- 
ing from a low occupied quantum state, each far removed from 
the Fermi level. (The inverted region effect should, however, 
occur for the electron transfer when the electrode is a narrow 
band semiconductor.) 

After these initial electron transfer studies, which were 
based on a dielectric continuum approximation for the solvent 
outside the first coordination shell of each reactant, I intro- 





Figure 8. 


Inverted region effect in chemical electron transfer reactions 
(Miller, et al, ref. 3) 
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duced a purely molecular treatment of the reacting system. 
Using statistical mechanics, the solvent was treated as a col- 
lection of dipoles in the 1960 paper, and later in 1965 a general 
charge distribution was used for the solvent molecules and for 
the reactants. At the same time I found a way in this 1960 paper 
of introducing rigorously a global reaction coordinate in this 
many-dimensional (N) coordinate space of the reacting sys- 
tem. The globally defined coordinate so introduced was 
equivalent to using U, - U,, the potential energy difference 
between the products plus solvent (U_) and the reactants plus 
solvent (U ) (cf A. Warshel, 1987). It was, thereby, a coordinate 
defined everywhere in this N-dimensional space. 

The free energy G_ of a system containing the solvent and 
the reactants, and that of the corresponding system for the 
products, G,, could now be defined along this globally defined 
reaction coordinate. (In contrast, in reactions such as that 
depicted by Figure 2, it is customary, instead, to define a 
reaction coordinate locally, namely, in the vicinity of a path 
leading from the valley of the reactants through the saddle 
point region and into the valley of the products.) 

The potential energies U_and U_ in the many-dimensional 
coordinate space are simple functions of the vibrational coor- 
dinates but are complicated functions of the hundreds of 
relevant solvent coordinates: there are many local minima 
corresponding to locally stable arrangements of the solvent 
molecules. However, I introduced a “linear response approxi- 





Figure 9. 


A favored formation of an electronically excited state of the 
products (Marcus and Siddarth, ref. 2) 
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mation,” in which any hypothetical change in charge of the 
reactants produces a proportional change in the dielectric 
polarization of the solvent. (Recently, I utilized a central limit 
theorem to understand this approximation better beyond sim- 
ple perturbation theory, and plan to submit the results for 
publication shortly.) With this linear approximation the free 
energies G_ and G, became simple quadratic functions of the 
reaction coordinate. 

Such an approach had major consequences. This picture 
permitted a depiction of the reaction in terms of parabolic free 
energy plots in simple and readily visualized terms, as in 
Figure 5. With them the trends predicted from the equations 
were readily understood. It was also important to use the free 
energy curves, instead of oversimplified potential energy pro- 
files, because of the large entropy changes which occur in 
many electron transfer cross-reactions, due to changes in 
strong ion-polar solvent interactions, (The free energy plot is 
legitimately a one-coordinate plot while the potential energy 
plot is at most a profile of the complicated U_ and U, in 
N-dimensional space. 

With the new statistical mechanical treatment of 1960 and 
1965 one could also see how certain relations between rate 
constants initially derivable from the dielectric continuum- 
based equations in the 1956 paper could also be valid more 
generally. The relations were based, in part, on Equations (4), 
(5) and (initially via (7) and (8)) on the approximate relation 


hio=5 Any +A) (9) 


where A> is the A for the cross-reaction and the Ay and Lu are 
those of the self-exchange reactions. 


Predictions 


In the 1960 paper I had listed a number of theoretical 
predictions resulting from these equations, in part to stimulate 
discussion with experimentalists in the field at a Faraday 
Society meeting on oxidation-reduction reactions, where this 
paper was to be presented. At the time | certainly did not 
anticipate the subsequent involvement of the many experi- 
mentalists in testing these predictions. Among the latter was 
one which became one of the most widely tested aspects of the 
theory, namely, the “cross-relation.” This expression, which 
follows from Egs. (5) and (9), relates the rate constant k., of 
a cross-reaction to the two self-exchange rate constants, k, , 
and k,,, and to the equilibrium constant K, , of the reaction. 


kya = (ky Koz Ki2f,2)” (10) 


where f,, isa known function of k, ,,k,, and K, , and is usually 
close to unity. 

Another prediction in the 1960 paper concerned what I 
termed there the inverted region: In a series of related reac- 
tions, similar in 1 but differing in AG°, a plot of the activation 





free energy AG" vs. AG° is seen from Eq.(5) to first decrease 
as AG° is varied from 0 to some negative value, vanish at AG° 
= i, and then increase when AG? is made still more negative. 
This initial decrease of AG’ with increasingly negative AG° is 
the expected trend in chemical reactions and is similar to the 
usual trend in “Bronsted plots” of acid or base catalyzed 
reactions and in “Tafel plots” of electrochemical reactions. I 
termed that region of AG° the “normal” region. However, the 
prediction for the region where AG® A, the “inverted region,” 
was the unexpected behavior, or at least unexpected until the 
present theory was introduced. 

This inverted region is also easily visualized using Figures 
6 and 7: Successively making AG° more negative, by lowering 
the products’ G curve vertically relative to the reactant curve, 
decreases the free energy barrier AG" (given by the intersection 
of the reactants’ and products’ curves): that barrier is seen in 
Figure 6 to vanish at some AG? and then to increase again. 

Other predictions dealt with the relation between the 
electrochemical and the corresponding self-exchange electron 
transfer rates, the numerical estimate of the reaction rate 
constant k and, in the case of non-specific solvent effects, the 
dependence of the reaction rate on solvent dielectric properties. 
The testing of some of the predictions was delayed by an extended 
sabbatical in 1960-61, which I spent auditing courses and attend- 
ing seminars at the nearby Courant Mathematical Institute. 


Comparisons of Experiment and Theory 


Around 1962 during one of my visits to Brookhaven 
National Laboratory, I showed Norman Sutin the 1960 predic- 
tions. Norman had either measured via his stopped-flow appa- 
ratus or otherwise knew rate constants and equilibrium 
constants which permitted the cross-relation Eq. (10) to be 
tested. There were about six such sets of data which he had 
available. I remember vividly the growing sense of excitement 
we both felt as, one by one, the observed k,,"s more or less 
agreed with the predictions of the relation. I later collected the 
results of this and of various other tests of the 1960 predictions 
and published them in 1963. Perhaps by showing that the 
previously published expressions were not mere abstract for- 
mulae, but rather had concrete applications, this 1963 paper, 
and many tests by Sutin and others, appear to have stimulated 
numerous subsequent tests of the cross-relation and of the 
other predictions. A few examples of the cross-relation test are 
given in Table I. 

The encouraging success of the experimental tests given 
in the 1963 paper suggested that the theory itself was more 
general than the approximations (e.g., solvent dipoles, un- 
changed force constants) used in 1960 and stimulated me to give 
a more general formulation (1965). The latter paper also contains 
a unified treatment of electron transfers in solution and at metal 
electrodes, and served, thereby, to generalize my earlier (1957) 
treatment of the electrochemical electron transfers. 





Table 1. 
Comparison of Calculated and Experimental k12 Values 





ky. (M~'sec™") 


Reaction Observed Calculated 





IrCl,?~ + WI(CN),*~ 
IrCl,?~ + FeiCN),*~ 
IrCi,?~ + MoiCN),*~ 
MoiCN),’~ + W(CN),*~ 
MoiCN),’~ +FeiCN),*~ 3.0 x 10° 
FeiCN),’~ +WICN),* 4.3x10* 
Ce’ + WiCN),*~ > 10° 

Ce'* + Fei(CN),*~ 1.9 x 10° 8 x10° 
Ce'¥ + Mo(CN),*~ 1.4 10’ 1.3 10’ 
L-Co{( — )PDTA}?~ + Feibipy),’* 81x10 2 >10° 
L-Fe{((—)PDTA}*- +Coi\EDTA)™ 1.3 x 10! 1.3 10! 
L-Fe[( — )PDTA]*~ +Cotox),’ 2.2x« 10 1.0x 10° 
CriEDTA)*~ +Fe(EDTA)~ > 10° 10° 
CriEDTA)*~ +Co(/EDTA)~ =3 x10° 4 x10’ 
Fe(EDTA)*?~ +MniCyDTA)~ =~4 x10’ 6 x10 
Co\EDTA)?~ +MniCyDTA) 9 x107' 2.1 
Fe(PDTA)?~ +CoiCyDTA)~ 1.210! 1.8x 10! 
Cotterpy),’* +Co(bipy),’* 6.410 3.210 
Coiterpy),’* +Co(phen),’* 2.810 Lix1e 
Coiterpy),’* +Co(bipyH,O),’* 68x10 6.4x 10° 
Coiterpy),?* +Co(phen\(H,0),’* 14x10 6.4 10° 
Coiterpy),’* +Co(H,0),’* 7.4 10° 2 x10" 
Fe(phen),’* +MnO,~ 6 x10 4 x10 
Fe(CN),*~- +MnO,~ 1.3 10* $ x10 
V(H,O),?* +Ru(NH;),’* 15x10 42x10 
Ru(en),’* + Fe(H,O),’* 8.4.x 10° 4.2x10° 
Ru(NH;,),?* +Fe(H,O),’* 34X10 7.5x 10° 
Fe(H,O),’* +Mn(H,0),’* 1.5 10° 3 x10* 


6.1 x 10” 
3.8x 10° 
1.9 10° 
5.0 x 10° 


6.110’ 
7 x10° 
9 x10 
48x 10° 
2.9 10° 
6.3 x 10° 
4 x10 








*Bennett, Ref. 3. 





The best experimental evidence for the inverted region 
was provided in 1984 by Miller, Calcaterra and Closs, almost 
25 years after it was predicted. This successful experimental 
test, which was later obtained for other electron transfer reac- 
tions in other laboratories, is reproduced in Figure 8. Possible 
reasons for not observing it in the earlier tests are several-fold 
and have been discussed elsewhere. 

Previously, indirect evidence for the inverted region had 
been obtained by observing that electron transfer reactions 
with a very negative AG” may result in chemiluminescence: 
when the Gand G. curves intersect at a high AG” because of 
the inverted region effect, there may be an electron transfer to 
a more easily accessible G, curve, one in which one of the 
products is electronically excited and which intersects the G_ 
curve in the normal region ata low AG", as in Figure 9. Indeed, 
experimentally in some reactions 100% formation of an elec- 
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tronically excited state of a reaction product has been observed 
by Bard and coworkers, and results in chemiluminescence. 

Another consequence of Eq. (5) is the linear dependence 
of k,T In k on with a slope of 1/2, when |AG°/AI is small, and 
a similar behavior at electrodes, with AG° replaced by en, the 
product of the charge transferred and the activation overpoten- 
tial. Extensive verification of both these results has been 
obtained. More recently, the curvature of plots of In k vs. en, 
expected from these equations, has been demonstrated in 
several experiments. The very recent use of ordered organic 
molecular monolayers on electrodes, either to slow down the 
electron transfer rate or to bind a redox-active agent to the 
electrode, but in either case to avoid or minimize diffusion 
control of the fast electron transfer processes, has considerably 
facilitated this study of the curvature in the In k vs. en) plot. 

Comparison of experiment and theory has also included 
that of the absolute reaction rates of the self-exchange reac- 
tions, the effect on the rate of varying the solvent, an effect 
sometimes complicated by ion pairing in the low dielectric 
constant media involved, and studies of the related problem of 
charge transfer spectra, such as 


DA + hv — D*A™ (11) 


Here, the frequency of the spectral absorption maximum is 
given by 


WVimax = A + AG®. 





Figure 10. 


Comparison of isotopic exchange electron transfer rates in solu- 
tion, covering 20 orders of magnitude, with rates of correspond- 
ing electron transfers at metal electrodes. (Cannon, ref. 2) 
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Comparisons with Eq. (12), using Eq.(7) for A, have included 
those of the effects of separation distance and of the solvent 
dielectric constant. 

Comparisons have also been made of the self-exchange 
reaction rates in solution with the rates of the corresponding 
electron transfer reactions at electrodes. An example of the 
latter is the plot given in Figure 10, where the self-exchange 
rates are seen to vary by some twenty orders of magnitude. 
The discrepancy at high k’s is currently the subject of some 
reinvestigation of the fast electrode reaction rates, using the 
new nanotechnology. Most recently, a new type of interfacial 
electron transfer rate has also been measured, electron transfer 
at liquid-liquid interfaces. In treating the latter, I extended the 
“cross relation” to this two-phase system. It is clear that much 
is to be learned from this new area of investigation. (The study 
of the transfer of ions across such an interface, on the other 
hand, goes back to the time of Nernst and of Planck, around 
the turn of the century.) 


Other Applications and Extensions 


As noted in Figure 1, one aspect of the electron transfer 
field has been its continued and, indeed, ever-expanding 
growth in so many directions. One of these is in the biological 
field, where there are now detailed experimental and theoreti- 
cal studies in photosynthetic and other protein systems. The 
three-dimensional structure of a photosynthetic reaction cen- 
ter, the first membrane protein to be so characterized, was 
obtained by Deisenhofer, Michel and Huber, who received the 
Nobel Prize in Chemistry in 1988 for this work. A bacterial 
photosynthetic system is depicted in Figure 11, where the 
protein framework holding fast the constituents in this reaction 
center is not shown. 

In the photosynthetic system there is a transfer of elec- 
tronic excitation from “antenna” chlorophylls (not shown in 
Figure 11) to a special pair BChl,. The latter then transfers an 
electron to a pheophytin BPh within a very short time (~3 
picoseconds) and from it to a quinone Q, in 200 psec and 
thence to the other quinone Q,- (Other chemical reactions then 
occur with these separated charges at each side of the mem- 
brane, bridged by this photosynthetic reaction center.) 

To avoid wasting the excitation energy of the BChl,* 
unduly it is necessary that the of this first electron transfer to 
BPh be small. (It is only about 0.25 eV out of an overall 
excitation energy of BChl,* of 1.38 eV.) In order that this 
electron transfer also be successful in competing with two 
wasteful processes, the fluorescence and the radiationless 
transition of BChL,* , it is also necessary that for that first 
electron transfer step be small and hence, by Eq. (5), that the 
i be small. The size of the reactants is large, and the immediate 
protein environment is largely nonpolar, so leading to a small 
i (cf Eq.(7)). Nature appears, indeed, to have constructed a 
system with this desirable property. 





Furthermore, to avoid another form of wasting the energy, 
it is also important that an unwanted back electron transfer 
reaction from the BPh to the BChL,* not compete successfully 
with a second forward electron transfer step from BPh toQ,. 
That is, it is necessary that the back transfer, a “hole-electron 
recombination” step, be slow, even though it is a very highly 
exothermic process (~1.leV). It has been suggested that the 
small A (~0.25eV) and the resulting inverted region effect play 
a significant role in providing this essential condition for the 
effectiveness of the photosynthetic reaction center. 

There is now a widespread interest in synthesizing sys- 
tems which can mimic the behavior of nature’s photosynthetic 
systems, and so offer other routes for the hamessing of solar 
energy. The current understanding of how nature works has 
served to provide some guidelines. In this context, as well as 
that of electron transfer in other proteins, there are also rele- 
vant experiments in long range electron transfer. Originally the 
Studies were of electron transfer in rigid glasses and were duc 
to Miller and coworkers. More recently the studies have in- 
volved a donor and receptor held together by synthetically 
made rigid molecular bridges. The effect of varying the bridge 
length has been studied in the various systems. A theoretical 
estimate of the distance dependence of electron transfers in a 
photosynthetic system was first made by Hopfield, who used 
a square barricr model and an approximate molecular estimate 
of the barrier height. 

Recently, in their studies of long range electron transfer 
in chemically modified proteins, Gray and coworkers have 
studied systematically the distance or site dependence of the 
electronic factor, by attaching an appropriate electron donor 
or acceptor to a desired site. For cach such site the reactant 
chosen should be such that -AG® = A, i.e., which has a k at the 
maximum of the In k vs. -AG®° curve (cf Egs. (4) -(5)). The 
value of k then no longer depends on a AG”. Since AG" is 
distance-dependent (cf Eq.(7)), it is particularly desirable to 
make AG* = 0, so that the relative k’s at the various sites now 
reflect only the electronic factor. Dutton and coworkers have 
treated data similarly for a number of reactions by using, where 
possible, the k at the maximum of each In k vs. AG® curve. Of 
particular interest in such studies is whether there is a simple 
exponential decrease of the electronic factor on the separation 
distance between donor and acceptor, or whether there are 
deviations from this monotonic behavior, due to local struc- 
tural factors. 

In a different development, the mechanism of various 
organic reactions has been explored by several investigators, 
notably by Eberson (ref. 2), in the light of current electron 
transfer theory. Other organic reactions have been explored by 
Shaik and Pross, in their analysis of a possible electron transfer 
mechanism vs. a conventional mechanism, and by Shaik et al 
(ref.2). 

Theorctical calculations of the donor-acceptor electronic 
interactions, initially by McConnell and by Larsson, and later 


by others, our group among them, have been used to treat 
long-range electron transfer. The methods have recently been 
adapted to large protein systems. In our studies with Siddarth 
we used an “artificial intelligence” searching technique to 
limit the number of amino acids used in the latter type of study. 

Another area of much current activity in electron transfers 
is that of solvent dynamics, following the pioneering treatment 
for general reactions by Kramers (1940). Important later de- 
velopments for electron transfer were made by many contribu- 
tors. Solvent dynamics affects the electron transfer reaction 
rate when the solvent is sufficiently sluggish. As we showed 
recently with Sumi and Nadler, the solvent dynamics effect 
can also be modified significantly, when there are vibrational 
(A,) contribution to A. 

Computational studies, such as the insightful one of David 
Chandler and coworkers on the Fe** + Fe’ self-exchange 
reaction, have also been employed recently. Using computer 
simulations they obtained a verification of the parabolic G 
curves, even for surprisingly high values of the fluctuation in 
G. They also extended their studies to dynamical and quantum 
mechanical effects of the nuclear motion. Studies of the quan- 
tum mechanical effects on the nuclear motion on electron 
transfer reactions were initiated in 1959 by Levich and Dogo- 
nadze, who assumed a harmonic oscillator model for the polar 
solvent medium and employed perturbation theory. Their 
method was related to that used for other problems by Huang 
and Rhys (1951) and Kubo and Toyozawa (1954). 

There were important subsequent developments by vari- 
ous authors on these quantum effects, including the first dis- 





Figure 11. 


Redox-active species involved in the initial charge separation 
for a photosynthetic bacterium (cf Deisenhofer el al, ref. 3; and 
Yeates et al, ref 3), with labels added, to conform to the 
present text; they include a missing Qg 
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cussion of quantum effects for the vibrations of the reactants 
by Sutin in 1962 and the important work of Jortner and 
coworkers in 1974-75, who combined a Levich and Dogo- 
nadze type approach to treat the high frequency vibrations of 
the reactants with the classical expression which I described 
earlier for the polar medium. These quantum effects have 
implications for the temperature dependence of k, among other 
effects. Proceeding in a different (classical) direction Savéant 
recently showed how to extend Eq.(5) to reactions which 
involved the rupture of a chemical bond by electron transfer 
and which he had previously studied experimentally: M(e) + 
RX — M+R-+X, where R is an alkyl group, X a halide and 
M a metal electrode. 

A particularly important early development was that by 
Taube in the 1950s, who received the Nobel Prize for his work 
in 1983. Taube introduced the idea of different mechanisms 
for electron transfer outer sphere and inner sphere electron 
transfers, which he had investigated experimentally. His ex- 
perimental work on charge transfer spectra of strongly inter- 
acting systems (“Creutz-Taube” ion, 1959, 1973) and of 
weakly interacting ones has been similarly influential. Also 
notable has been Hush’s theoretical work on charge transfer 
spectra, both of intensities and absorption maxima (1967), 
which supplemented his earlier theoretical study of electron 
transfer rates (1961). 

There has been a “spin-off” of the original electron trans- 
fer theory to other types of chemical reactions as well. In 
particular, the vs relation and the cross-relation have been 
extended to these other reactions, such as the transfer of atoms, 
protons, or methyl groups. (Even an analog of Eqs. (5) and (9), 
but for binding energics instead of energy barriers has becn 
introduced to relate the stability of isolated proton-bound 
dimers AHB* to those of AHA* and BHB*!) 

Since the wansfcr of these nuclei involves strong elec- 
tronic interactions, it is not well represented by intersecting 
parabolic free energy curves, and so a different theoretical 
approach was needed. For this purpose I adapted (1968) a 
“bond-cnergy-bond-order” model of H. Johnston, in order to 
treat the problem for a reaction of the type given by Eq.(4). 
The resulting simple expression for is similar to Eq.(5), when 
|AG"/A is not large ( 1/2), but differs from it in not having any 
inverted region. It has the same property as that given by 
Eq.(9), and has resulted in a cross-relation analogous to Eq. 
(10). The cross-relation has been tested experimentally for the 
transicr of methyl groups by E. Lewis, and the vs relation has 
been used or tested for other transfers by Albery and by 
Kreevoy and their coworkers, among others. 

It is naturally gratifying to see one’s theories used. A 
recent article, which showed the considerable growth in the 
use of papers such as the 1956 and 1964 articles (ref. 5), points 
up the impressive and continued vitality of the ficld itself. The 
remarks above on many arcas of electron transfer and on the 
spin-off of such work on the study of other types of reactions 
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represent a necessarily brief picture of these broad-based 
investigations. 
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Introduction 


During the past decade, the uses of computer simulation 
and modeling techniques in scientific research have grown 
enormously. The remarkable power of the computer, when 
harnessed through intelligent and efficient simulation method- 
Ologies, has led in many instances to a modification of the 
scientific method. In particular, modern scientific research 
often involves a dynamic interplay between theory, experi- 
ment, and computer simulation. This “triad”, when properly 
balanced, can offer a considerable advantage over the more 
traditional interplay between theory and experiment alone. 

The exciting field of electron transfer (ET) has also begun 
to be influenced by computer simulation. Most of these simu- 
lations have focused on classical calculations of the famous 
Marcus free energy curves’ for realistic models of several ET 
systems.” While these realistic classical calculations have 
confirmed the validity of the linear response assumptions 
employed by Marcus in his original theory,' a first principles 
computational treatment of electron transfer requires a method 
which can directly and simultaneously address both the influ- 


ence of medium/vibronic effects and the quantum mechanical 
character of the electron being transferred between donor and 
acceptor. (A similar statement also holds true for the study of 
other quantum rate phenomena such as condensed phase pro- 
ton and hydride transfer reactions, hydrogen diffusion in and 
on materials, etc.) For kinetic processes involving heavier 
particles (e.g., most chemical reactions), such a computational 
method has existed for many years in the form of the well 
known classical transition state theory (TST).’ Only recently, 
however, has the formulation of a quantum mechanical TST 
(QTST) been proposed*® which provides a consistent and 
computationally powerful generalization of the classical TST. 
This quantum theory, which is based on Feynman path inte- 
grals,’ will be briefly described below and representative 
applications to computer simulation discussed. 


Classical Rate Theory 


To begin, it is useful to first review the structure of 
classical TST which provides a basis for understanding the 
attributes of a more general quantum theory.*® The TST ap- 
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proximation to the classical rate constant k_, is often expressed 
as 


pst kpT 


d Wp. exp (-BF7,), (1) 


where Q,.1 is the classical partition function of the total system 
in the reactant state, and the excess classical free energy at the 
dividing surface, F”,, is given by 


Ft, =—kgT In{ Q2) (m/2nh’ 8) (2) 


where m is the effective mass of the reaction coordinate 
motion and B equals (k,T)'. In Eq. (2), Q*, is the classical 


phase space density at the position of the dividing surface 
along the reaction coordinate q, given by 


> l 


2.,.= 
" (2nh ) 


w Jae, Jap, fda fax 8 (a" - @) exp (-BH), 3) 


In the above equations, q* is the transition state value of the 
reaction coordinate, N is the number of degrees-of-freedom of 
the entire system, and H is the classical Hamiltonian. The 
coordinates x and their conjugate momenta p, are for all 
degrees-of-freedom other than the reaction coordinate. The 
key factor in Eq. (1) is the dimensionless classical excess free 


energy factor exp(—BF-, ) which, in relation to the reactant 


State partition function, leads to the activation factor in TST. 
The latter factor dominates the temperature dependence of the 
rate constant. The pre-exponential factor in Ea. (1) is the 
mathematical result of the TST assumption? (i.e., all trajecto- 
ries which reach the transition state configuration from the 
reactants will continue onward to the products). Marcus was 
able to employ classical TST in his original formulation of the 
theory of adiabatic electron transfer by choosing the collective 
solvation variable as the reaction coordinate.’ Not surpris- 
ingly, the theoretical basis for classical computer simulations 
of ET reactions” has been the classical TST perspective (i.c., 
these simulations focus on the classical solvent activation free 
energy for ET). 

The most important feature of the path integral formula- 
tion of quantum TST is that it specifies both a quantum 
activation factor and a TST prefactor,*® both of which are the 
analogs to the classical TST quantities in Eqs. (1)-(3). As in 
classical TST, the activation factor is the most important 
feature of the quantum theory since it dominates the tempera- 
ture dependence of the rate constant and can be determined 
without any dynamical information (i.e., by using statistical 
mechanics). The quantum activation factor also captures any 
tunneling and zero-point energy effects which can give rise to 
substantial deviations from the classical TST rate constant. 
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Feynman Path Integrals 


The approach to a quantum mechanical TST described in 
the following section*® draws heavily on Feynman’s path 
integral formulation of quantum statistical mechanics.’ The 
path integral formalism is, in fact, the essential basis for the 
conceptual and computational power of the quantum TST. The 
Feynman path integral for, e.g., the partition function of the 
system is derived by first expressing the trace of the 
Boltzmann operator in the coordinate representation, i.c., as 


Q = Tre ™ = [dg <qie* Ig. (4) 


It should be noted here that for simplicity Eq. (4) has been 
written for a single degree-of-freedom, but an extension to 
multiple dimensions and/or particles is straightforward. The 
Boltzmann operator can now be written as a product of P 
components such that 


BH = (eel (5) 


where e equals § / P. The partition function is then rewritten 
exactly as 


Q= Jaq <qie Me eg. (6) 


It is now possible to insert P complete sets of coordinate space 
states, each given formally by 


Jag | g><q'=1 (7) 


between the individual pieces of the Boltzmann operator. The 
partition function is then given by 


Q=f day dg... f dap<q, lem 


| qo >< qn le 1g, >... <qple™ Iq, >, 


where a numbering scheme from | to P has been adopted to 
keep track of the different integrations over complete sets of 
coordinate states. For large values of P, the expressions for the 
matrix elements <q.le“"ig, > are known’ since the Boltzmann 
operator exp(-€H) = exp(-H/k, T,,), with T, =T xP, is equiva- 
lent to one for a very high temperature system. Upon substi- 
tution of these expressions for each matrix element in Eq. (8), 
the partition function is given by the classical-like expression 


Q = lim (m/2nh 26) faa, = dap exp[-BV 4 (4;,---Gp))(9) 
P —se0 


where the effective potential is given by 





P 
mP a 
Vi Qe =D laa Guy) +pY a) (10) 
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Remarkably, the seemingly daunting quantum expression in 
Eq. (8) is now reduced to an effective classical partition 
function similar to a “polymer” consisting of P quasiparticles 
located at the coordinates {q,}.* The quasiparticles interact 
with each other through nearest-neighbor harmonic potentials, 
and each quasiparticle feels the potential energy function via 
the term V(q,) / P. In addition, the quasiparticle coordinates 
are cyclic so that qi equals q, ,,. Figure 1 schematically depicts 
the discretized path integral representation of an electron 
solvated in a heavy classical fluid such as argon. 

The path integral “isomorphism” between a single quan- 
tum mechanical particle and a polymer of effective classical 
quasiparticles is the centerpiece of modern computational path 
integral techniques.*” If one chooses a large enough value of 
P in Eq. (8) to insure convergence, the evaluation of quantum 
equilibrium quantities is simply mapped onto an isomorphic 
classical calculation. Such calculations can be readily per- 
formed on a computer using Monte Carlo or Molecular Dy- 
namics techniques. Many studies have now appeared in the 
literature in which path integral Monte Carlo has been em- 
ployed to simulate a quantum particle (or particles) in very 
complex environments involving hundreds of atoms and/or 
molecules.’ 

The level of abstraction in the path integral theory in- 
creases when one takes the formally exact P — @ limit in Eq. 
(8). In that case, the path integral for the partition function 
becomes a true functional integral,’ i.e., 





Figure 1. 


A schematic diagram of the isomorphic path integral “poly- 
mer” for an electron in a condensed phase environment of 
heavy classical particles. Equilibrium properties are calcu- 
lated by performing an appropriate Boltzmann-weighted sum 
over all possible configurations of the path integral polymer 
and bath particles 














Q=J ... [Dg (x) expl-S tq mV/r }, (11) 


where q(t) is the path of the quantum particle in imaginary 
time and S[{q(t),x(t)] is the imaginary time action functional’ 


S [9(t)] = fra 3 ary +V Iq col} (12) 


The above exact path integral representation is necessary when 
one develops new formalism or analytic theory, but the dis- 
crete “polymer” picture embodied in Eq. (9) and Fig.1 is all 
this is required for one to understand and visualize the path 
integral perspective.® 


Quantum Rate Theory 


In general terms, the calculation of the quantum activated 
rate constant can be accomplished with a general TST-like 
theory if such a theory provides expressions for both the 
quantum activation free energy and the prefactor as in the 
classical theory [cf. Eqs. (1)-(3)]. The latter factor must be the 
QTST equivalent to the classical pre-exponential factor in Eq. 
(1). Such a theory has been under development during the past 
few years.** 

The progress towards a path integral quantum TST is 
partly due to the importance of the path “centroid” variable 
G, in path integration.*® This variable is defined along the q 
direction, for example, by the expression 


eS B 
=— dt q(t) (13) 
i= J aca 


or, in the discretized path integral notation,* by 


P 

- 1 

1o=p ¥4; (14) 
i=] 


The centroid variable is seen from Eq. (13) to be the imaginary 
time average of the quantum path q(t) or, in the discretized 
path integral expression in Eq. (14), to be the center-of-mass 
of the discretized “polymer” of classical quasiparticles having 
coordinates {q,}. 

The path integral quantum TST rate constant is given by*® 


kpT 
OTST iO, exp (-BF2) (15) 


where the excess path integral centroid free energy at the 
dividing surface is defined as 


F* =—kgT In { Q2/(m/2nh *B)*} (16) 
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The reduced centroid density Q” at the dividing surface along 
the reaction coordinate q is the central quantity in determining 
the value of the rate constant. The latter quantity is given by**® 


Q" =f... J Da (t) Dx (x) 30g" -G,) exp {-Sad.x(@) Vm }.17) 


The above formula may be readily applied to a variety of 
problems because, in the spirit of classical transition state 
theory,’ no explicit quantum dynamical information is re- 
quired to estimate the rate constant. Additionally, the activa- 
tion free energy in Eq. (15) can be calculated with discretized 
path integral Monte Carlo or Molecular Dynamics techniques” 
for highly complex systems. 


Representative Applications 


In this section, two applications of the path integral QTST 
are described. The interested reader is referred to papers 
containing other applications path integral QTST, including 
homogeneous electron transfer theory'®"' and simulation,'*! 
proton transfer theory’* and simulation,'*'> hydrogen diffu- 
sion metal surfaces, '° molecular diffusion on metals,’ and the 
theory of condensed phase effects in quantum activated dy- 
namics.4>"!!+!8 


Proton Transfer Reactions in Polar 
Solvents 

Proton transfer reactions in condensed phases are of con- 
siderable importance in chemistry. Since proton transfer in- 


volves a redistribution of solute electronic charge density, a 
substantial contribution to the activation free energy may 





Figure 2. 


A graphical depiction of the discretized Feynman path integral 
representation of a quantum particle being transferred be- 
tween two heavy classical-like particles. The reacting system 
is immersed in a simple classical fluid. The centroid of the 
quantum electron paths is constrained to be at the transition 
State in order to calculate the relevant quantum activation free 
energy. 
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come from solvent reorganization effects just as in electron 
transfer’ (i.c., solvent fluctuations are necessary to create a 
degeneracy between the two proton binding states so that the 
proton can transfer). In some ways, however, proton transfer 
is even more complicated than ET, since intramolecular vibra- 
tions may also play a central in modulating the proton transfer 
process.'*-!9 At the very least, proton transfer reactions in polar 
solvents provide fascinating systems through which to study 
the interplay between solvent activation dynamics, intra- 
molecular mode coupling, and quantum tunneling effects. Not 
surprisingly, they have attracted a considerable degree of 
experimental and theoretical attention in recent years. 

In this section, the application of path integral QTST to 
proton transfer reactions described in Refs. 13 and 14 is 
reviewed. The benefit of this approach is that it includes both 
the solvent and intramolecular contributions to the proton 
transfer activation free energy while treating any number of 
particles in the system quantum mechanically (including, of 
course, the proton). The method also bridges the adiabatic and 
nonadiabatic limits of the activation factor and is applicable 
to situations where the proton itself becomes thermally acti- 
vated. In this implementation of path integral QTST, the actual 
proton transfer coordinate is treated as the reaction coordinate 
(e.g., a proton asymmetric stretch mode),'* rather than the 
collective solvent polarization coordinate as in a Marcus-like 
approach.'? This formulation of the problem allows the quan- 
tum tunneling of the proton to be explicitly treated on the same 
footing as the solvent activation free energetics. It also facili- 
tates detailed examination of specific intramolecular interac- 
tions between the proton transfer coordinate and other 
molecular modes in the charge transfer complex.'* 

According to the path integral QTST equation, the central 
quantity in any quantum activated rate process in the quantum 
activation free energy. This quantity can be cast into a form 
particularly well suited for computer simulation,’* i.e., 


AF* =—kgT In [P(q,9')] (18) 


where Pc(qr A q*) is the quantum probability to move the 
reaction coordinate centroid variable from the reactant con- 
figuration to the transition state. This probability is readily 
calculated14 by path integral Monte Carlo techniques9 com- 
bined with umbrella sampling.20 In the latter computational 
technique, a number of “windows” are set up which confine 
the path centroid variable of the reaction coordinate to differ- 
ent regions.14 These windows connect in a piece wise fashion 
the possible centroid positions in going from the reactant state 
to the transition state. A series of Monte Carlo calculations are 
then performed, one for each window, and the centroid prob- 
ability distribution in each window is determined. These indi- 
vidual window distributions are then smoothly joined to 
calculate the overall probability function Pc(gr A q*) in Eq. 
(18). Alternatively, the mean force on the path integral centroid 
reacuion coordinate could be calculated in the simulation and 





then integrated from q_ to q’ to give the activation free energy. 
In Fig. 2, a schematic picture depicts the path integral QTST 
representation of a quantum particle (e.g., an electron or 
proton) at the transition state between the donor and acceptor 
species in a classical fluid. The centroid of the path integral 
“polymer” for the quantum particle (cf. Eq. (9)] is constrained 
to be halfway between the donor and acceptor (i.e., at its 
transition state). The function P(q, — q°) in Eq. (18) is the 
probability for the path integral of the quantum particle to 
achieve the transition state configuration relative to the reac- 
tant State. 

By using path integral QTST approach, one can directly 
study the solvent contribution to a proton or an electron 
transfer reaction without making any classical assumptions. 
As an example, consider the model proton transfer reaction 


A*OS 57°05 .., Af A... HOS... A**]' 2A... HO -AO 


where “A” denotes some heavy atom. In Ref. 14, the path 
integral QTST was used to calculate the total quantum activa- 
tion free energy for the model three-body proton transfer 
complex solvated in the SPC/E model”! for water. Changes in 
the total activation free energy, and hence the reaction prob- 
ability, due to the proton tunneling were also a focus of the 
study. For more details of the model or the path integral Monte 
Carlo calculations, the reader is referred to Ref. 14. 

For the system studied in Ref. 14, the total quantum 
activation free energy curves as a function of the proton 
asymmetric stretch coordinate are depicted in Fig. 3. Shown 
are the activation curves for the complex in isolation and for 
the complex in the water solvent. The effect of the solvent in 
the total activation process is obvious, contributing ~5 
kcal/mole to the overall activation free energy.'* In fact, the 
solvent activation is apparently the rate determining step for 
this proton transfer system just as in many ET reactions. The 
activation free energy curves were also calculated in the clas- 
sical limit of the proton transfer reaction. By comparison with 
the fully quantum result, the effect of the proton tunneling on 
the total activation free energy was found to be large (~5 
kcal/mole), increasing the total reaction probability over the 
classical case by roughly a factor of 2400. A similar path 
integral calculation has been performed’? for the Fe?+/Fe*+ 
ET reaction in water using a simplified, but accurate, tight- 
binding model for the electron. Interestingly, the water libra- 
tional modes were found to exhibit a large degree of tunneling, 
leading to an unexpected quantum enhancement of the ET rate 
by a factor of 40 over the classical Marcus theory prediction. 


Heterogeneous Electron Transfer at 
Electrodes 


While path integral QTST allows one to directly simulate 
many outer-sphere homogeneous electron transfer reactions, 





Figure 3. 


Quantum activation free energy curves for the model proton 
transfer described in the text and Ref. 14. The free energy is 
calculated by virtue of Eq. (18). The solid line is for the solute 
immersed in the water solvent, while the dashed line is for the 
solute in isolation 








the simulation of heterogencous ET at electrodes represents an 
even greater challenge due to the many-electron nature of the 
problem, as well as the increased complexity of the molecular 
interactions in such systems. As a first step, electron transfer 
between an ion in solution and a metal electrode can be 
described??? using a version of a Hamiltonian referred to in 
the literature as the Anderson-Newns Hamiltonian.2*”’ This 
Hamiltonian can be written as 

H=H,.,+H (19) 


1 el 


where H_, is the Hamiltonian for the electrolyte solution 
(including solution-metal interactions) and H.,, is the “electron 
transfer” part of the Hamiltonian, given by 


4 » + : wa % 
H,=(€,+4E)n,+ > (en, +! c'c.+V. cc ) 


kk akak kaka) 
k y 


Here, €, is the vacuum energy level of the ion electron orbital 
which is involved in the electron transfer (orbital “a”) and an 
infinitely dilute solution is assumed for simplicity, though not 
necessary. This acceptor orbital has an occupancy number n_ 
equal to either 0 or 1, depending on whether or not it is empty, 
and AE is the instantaneous energy shift of the ion electron 
State | a> due to the interaction with the fluctuating solvent 
(i.e., the Marcus solvation variable' or “coordinate”’). The sum 
over k in Eg. (20) represents the electronic states in the 
semi-infinite metal electrode. The latter term clearly includes 
contributions arising from both the occupancy of the metal 
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Figure 4. 


Adiabatic free energy curves for the Fe**/Fe** heterogeneous 
electron transfer reaction between an ion and a Pt electrode. 
The parabolic fits to the Marcus diabatic curves are shown by 
the dashed lines 
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electron states as well as a piece which is a “transfer” interac- 
tion leading to the charge transfer between the metal states and 
the ion orbital (i.c., c’, and c, are electron creation and annihi- 
lation operators, respectively, and Vi, are the matrix elements 
which couple the relevant states). 

The quantum mechanical Hamiltonian of Eq. (19) is not 
yet in a useful form for the purpose of simulating heterogene- 
ous electron transfer. Following the work of Grimley,” and 
Muscat and Newns,”° as well as invoking the physically rea- 
sonable assumption that the solvent and electronic degrees of 
freedom are separable due to different timescales, one must 
first concentrate on the electronic Hamiltonian H, in Eq. (20). 
After a series of mathematical manipulations involving the 
resolvent operator for this Hamiltonian, the ground state elec- 
tronic energy E,(AE) can be analytically obtained as a function 
of four quantities: (1) the Fermi level of the metal electrode 
€,., (2) the vacuum energy level €, of the ion, (3) the Marcus 
solvation variable AE, and (4) the vacuum broadening of the 
ion orbital la> due to the interaction with the metal electronic 
states. The latter parameter A is given by the expression 


A(e)=n > 1 V4 17 Ke —-e,) (21) 
k 


and is evaluated at € = €,.. Because one can relate the Marcus 
solvauon variable AE to the nuclear configurations of the 
solvent molecules, one can use the overall adiabatic ground 
state Hamiltonian function H, = H_, + E,(AE) to describe the 
forces on the solvent molecules.22 
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At this point, simulations of heterogeneous electron trans- 
fer at metal electrodes can now be designed.”* One needs to 
first construct a realistic microscopic potential which de- 
scribes the interactions between the solvent molecules, the 
solvent and the electrode, the donor/acceptor and the solvent, 
and the donor/acceptor and the electrode.’ Once completed, 
umbrella sampling techniques can then be applied to calculate 
the adiabatic Marcus free energy curve as a function of AE and 
thus find the classical solvent activation free energy barrier to 
the charge transfer reaction involving the ionic orbital la>. The 
first such curves to ever be calculated” are depicted in Fig. 4 
for the Fe?*/Fe** redox reaction near a Pt(111) electrode. 
Shown by the dashed lines are the parabolic Marcus-like 
curves”? which best fit the diabatic product and reactant states. 
In Fig. 5, a “snapshot” is shown from the classical simula- 
tions”* which produced Fig. 4. The next step in this research 
will be to employ the path integral QTST along with a quan- 
tized version of the water solvent to search for any quantum 
solvent effects such as those found in the homogeneous 
Fe”*/Fe** ET simulations.'° 


Outlook for the Future 


The path integral QTST formulation*® reviewed in the 
present article represents one way to introduce the consider- 
able advantages of computer simulation into the study of 
electron transfer reactions and other quantum rate phenomena. 
Such simulations have helped to confirm many of the theoreti- 





Figure 5. 


A snapshot from the simulation of the Fe**/Fe** heterogene- 
ous electron transfer reaction between an ion and a Pt elec- 
trode. The single iron ion is the dark sphere, surrounded by 
water molecules, and is at a distance of 5 A from the Pt 
electrode surface 
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cal ideas developed by Marcus 40 years ago' and may indeed 
help to uncover some unexpected effects in the field of electron 
transfer. One real challenge for the future will be to simulate, 
on a very large scale, increasingly complex systems such as 
those found in heterogeneous ET reactions.”*”? If these efforts 
are to be truly influential, however, new “first principles” 
simulation methods will need to be developed which, for 
example, unify path integral QTST with novel methods” to 
calculate the electronic structure of materials. Given the cur- 
rent rapid development of computer simulation techniques for 
homogencous electron transfer and other quantum rate proc- 
esses, as well as the advances in the fundamental theory 
underlying those techniques, the outlook seems bright for the 


role of computer simulation in field of heterogeneous electron 
transfer. 
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Introduction 


The kinetics of electron-transfer (ET) processes is a sub- 
ject graced by a remarkable diversity as well as importance in 
chemistry and biology. Of the manifold types of chemical ET 
reactions, those occurring in solution and at electrode-solution 
interfaces have been the subject of a notably concerted and 
fruitful research effort. This activity was in large part spawned, 
and subsequently spurred, by the development of the theoreti- 
cal treatments by Rudy Marcus and others, now broadly 
known as “Marcus theory.”'© The happy circumstances that 
enabled the formulation of such quantitative theoretical rela- 
tionships for solution-phase ET kinetics having true predictive 
and interpretative value have acted as a considerable impetus 
to the development of an increasingly diverse and impressive 
range of experimental activities. These now encompass ther- 
mal and photoinduced charge transfers for inorganic, organic, 
and biological systems in homogeneous solution, ET at met- 
al-solution interfaces, and ET and photo-ET at semiconduc- 
tor-solution interfaces. Indeed, the ability of the Marcus 
theoretical framework to interrelate, interpret, and predict ET 
rates and other dynamical phenomena across such diverse 
classes of reactions constitutes a major achievement in con- 
temporary physical chemistry, and was undoubtedly responsi- 
ble in part for the award of the 1992 Nobel Prize in Chemistry 
to Rudy Marcus. 


Along with numerous other experimentalists, our re- 
search into electron transfer has been influenced centrally, 
persistently, and encouragingly not only by Marcus theory, but 
also by Rudy Marcus himself. In this article we outline some 
pertinent issues and experimental inquiries into solution-phase 
ET phenomena, illustrated in part by findings from our labo- 
ratories. Some emphasis is placed on solvent effects, for which 
the interaction between experiment and theory has been par- 
ticularly lively and beneficial. For the benefit of the more 
general reader, we preface this discussion with a short sum- 
mary of some relevant conceptual and theoretical material. 
The article also in this issue by Ratner contains further detailed 
information along these lines. 


Some Conceptual Background 


The central challenge, met so ably by Marcus theory, is 
to rationalize, interpret, and predict the seemingly dazzling 
array of rates observed for solution-phase ET reactions, corre- 
sponding to reaction half lives ranging from picoseconds to 
many months or years! A generalized expression for bimolecu- 
lar or electrochemical rate constants based on activated-com- 
plex theory is’ 


k = Z exp (-AG’/RT) (1) 
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where Z is the prefactor and AG” is the overall free energy of 
activation. The original Marcus formulations focussed atten- 
tion on the latter component, understandably so since one 
expects that the very wide range of observed rate constants can 
often be attributed primarily to differences in the activation 
energetics. The framework of Marcus theory provides an 
extremely useful way in which the various factors which 
influence AG" can be separated, identified, and understood. 
Particularly useful is the distinction between “intrinsic” 
and “extrinsic” (or thermodynamic) components, denoting 
contributions to AG” arising from factors present in the ab- 
sence, and additional presence, of the free-energy driving 
force, AG°. With the presumption of parabolic free-energy 
surfaces, the following well-known relation is obtained:° 


AG’ = (A + AG°)?/4A (2) 


In eq. 2, A is the so-called “reorganization energy”, equal to 
four times the “intrinsic barrier”, AG’. the component of AG" 
which remains in the absence of a driving force (see Figure 1). 
Kinetic-based tests of Eq. 2 have attracted much attention, 
formerly concerning the predicted quadratic log k - AG® 
dependence in the “normal” driving-force region (where -AG° 
< A)’, and latterly the diminishing rates anticipated with in- 
creasing driving force in the “inverted” region (where -AG° < 
A) (see below)*. While Eq. 2 has limitations under some 
conditions, the reasonable concordance often observed with 
experiment, especially in the “normal” free-energy region, is 
a notable feature of ET phenomena. 

Another central aspect of Marcus theory concerns the 
nature of the intrinsic barrier, AG"... The Marcus treatment of 
the outer-shell (solvent) contribution, AG", to the intrinsic 
barrier, based on a nonequilibrium dielectric-continuum 
model, has proven remarkably resilient; representative experi- 
menial tests based on the energetics of optical ET transitions 
are discussed below. 

Among the more significant consequences of the theoreti- 
cal treatments emphasized originally by Marcus, are the inter- 
relationships predicted between rate constants for different 
outer-sphere reactions featuring common redox couples.*° 
Probably the best known is a connection between the rate 
constants for homogencous-phase cross reactions and for the 
pair of constituent self-exchange processes. Also significant, 
however, is the relation predicted between the rate constants 
for corresponding self-exchange (K) and electrochemical 
exchange (k* ) reactions: 


Ke /AS Ss (ke,/A")* (3) 


where A° and A° are the relevant prefactors (vide infra). In eq. 
3, the extent of the predicted inequality depends on the extent 
of image stabilization of the electrochemical transition state. 
Another prediction that has proven useful in mapping reactiv- 
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Figure 1. 


Schematic representation of energy relationships for activated 
electron transfer. 














Reaction coordinate 





ity trends in related homogeneous-phase and electrochemical 
processes is* 


(K\/K}) 2 =(Ki/K)p (4) 


referring to electrochemical rate constants for a pair of redox 
couples at a fixed electrode potential, E, in comparison with 
homogeneous-phase rate constants involving the same pair of 
couples reacting with a common reagent (oxidant or reduc- 
tant), R. 

It is instructive to consider further the appropriate form of 
the prefactor terms appearing, for example, in Egs. 1 and 3. 
Presumably for simplicity as well as convenience, such terms 
for outer-sphere ET reactions were initially formulated on the 
basis of collision-theory expressions. More recently, it was 
recognized by Sutin,? Marcus,'° and ourselves"! that an alter- 
native “encounter preequilibrium” treatment can provide a 
more useful description of the preexponential factors in both 
outer- and inner-sphere ET processes. 

A general expression for the overall preexponential factor 
on this basis is" 

A=K,,1 Vp, (5) 
where K_ is an equilibrium constant (statistical factor) for 
forming the precursor state (i.e. with appropriate geometry for 
ET but prior to nuclear activation). The remaining pair of terms 
in Eq. 5 refers to the ET activation step itself: the nuclear 
frequency factor v_ (s') describes the net dynamics of ap- 
proaching the transition state (associated with nuclear reor- 





ganization), and the electronic transmission coefficient «_, 
denotes the fractional probability with which electron tunnel- 
ing (so to consummate the reaction) takes place once the 
transition state is formed. The magnitude of K, for inner- 
sphere reactions is clearly governed by the bonding estab- 
lished between the redox centers. For outer-sphere processes, 
simple statistical expressions can be deduced, depending on 
the reaction geometry. Thus for homogeneous-phase proc- 
esses between spherical reactants in the absence of electro- 
static work terms: 


K,= 4nNrSr (6) 


where N is Avogadro’s number, r is the reactant separation in 
the transition state, and Sr is a “reaction-zone thickness” 
denoting the range of internuclear separations over which 
electron transfer primarily takes place. A related expression 
applies to outer-sphere electrochemical reactions, only now 
K, = &r since the “coreactant” is a planar interface. The values 
of dr are anticipated typically to be around 1A, enabling at 
least approximate estimates of the statistical factor K, to be 
obtained. 

Another useful outcome of this “encounter preequili- 
brium” treatment is that it enables the reactivities of related 
homogeneous-phase and electrochemical processes to be 
compared in an especially direct fashion.'* Specifically, the 
rate constant for a homogeneous reaction involving a pair of 
redox couples, k,,", is usefully compared with the electro- 
chemical rate constant involving one of the reactants, k,°, 
measured st an electrode potential equal to the formal potential 
for the coreactant redox couple. In view of Eq. 6, the homo- 
geneous-phase and electrochemical reactivities can be evalu- 
ated on a common basis by multiplying the latter by 4xNr’. 


4xNr’k, = k,,°, having the same units as k,," and equalling the 
elecuochemical reactivity in the (hypothetical) circumstance 
where the electrode “coreactant” offers the same reaction 
geometry as the homogeneous coreactant. Comparison between 
the values of k,,° and k,," can provide straightforward insight 
into the differing environmental factors that influence electro- 
chemical and homogeneous-phase reactivities. This procedure is 
most useful if the latter rate constants are corrected for the 
inner-shell (vibrational) activation associated with the homoge- 
neous-phase coreactant. For the resulting k_" values, any differ- 
ences in comparison to k,, © values will reflect only differences 
in the outer-shell (solvent) reorganization and/or electronic over- 
lap factors between the two reaction environments. 

An extensive comparison secently undertaken for inor- 
ganic reactions shows that k,.° typically exceeds k" __ , the 
former quantities commonly’ being as much as 107-10" fold 
larger than the latter.'? Some illustrative rate data of this type 
are given in Table I. These substantial rate differences, illus- 
trating the often markedly more facile nature of electrochemi- 
cal outer-sphere processes, most likely reflect the greater 
efficiency of electron tunneling (i.e. larger x_,) and possibly 
smaller outer-shell reorganization energy prevalent at metal- 
solution interfaces. 

The physical factors that control the preexponential fac- 
tor, and especially the nuclear reaction dynamics as reflected 
in v_ (Eq. 5), have been the subject of considerable recent 
interest, particularly with regard to the role of solvation dy- 
namics.'* In spite of the apparent form of Eq. 5, the nuclear 
dynamics will only influence the ET rate significantly when 
the donor-acceptor electronic coupling is sufficiently strong 
so that K,, — 1 (so-called “adiabatic pathway”). For weaker 
coupling (such that K_, << 1) the reaction dynamics will be 
determined instead by electron tunneling. (In other words, 


This yields an “equivalent second-order” rate constant, «Vv (= v,) will be proportional to (H a where H., is the 





Table 1. Comparison Of “second-order Equivalent” Electroreduction Rates For Selected Reactions At Mercury Surface 
With Corresponding Homogeneus-phase Rates 





Ke Kp? 
Solvent cm = mM” s' 
H20 4X10° 13x10 
H20 2.5 85x10 
H20 5 -2x 10° 
H20 0.06 2x 10° 
DMF 2.0 9x10’ 
TMU 0.35 15x10’ 


Oxidant 

Co (NH3) &* 
Ru (NH3)a" Ru (NHs)&* 
Fe (OHz 2 Ru (NH3)&" 
Or Ru (NH3)&" 
Cp2Co* Cp2Co 
Cp2Co* Cp2Co 

See ref. 12 for details. 


Reductant 
Ru (NH3)&* 





a Work-corrected electrochemical rate constant for reduction of oxidant measured at formal potential of 
reductant redox couple. 


b. “Equivalent second-order” rate constant for electroduction, derived from k*; (see text and ref. 12). 
c. Second-order rate constant for homogeneous reaction. 
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electronic coupling matrix element.) Nonetheless, there are 
many experimental ET systems (including outer-sphere sys- 
tems), where reaction adiabaticity is achieved (or ap- 
proached). Consequenily, the dynamics of activated nuclear 
barrier crossing, as denoted by v_, should often hold sway over 
the preexponential factor. 

Generally speaking, the overall nuclear frequency factor 
v,, will be determined by the dynamics of both the reactant 
vibrational (inner-shell) and solvation (outer-shell) compo- 
nents of the activation barrier. In the transition-state theory 
(TST) limit, v_ is usually predicted to be governed chiefly by 
the fastest dynamics (usually reactant vibrations). The role of 
solvation dynamics, therefore, was originally thought to be 
unimportant. Recently, however, it has become apparent that 
significant or even substantial deviations from TST may occur 
as a consequence of so-called “solvent friction”, whereby the 
dynamics of collective solvent motion (restricted rotations, 
etc.) required to surmount the ET transition state fall below the 
dynamics for the rotation of individual solvent dipoles (solvent 
inertial TST limit). Unlike the TST case, the frictional solvent 
motion is predicted to exert an important or even dominant 
influence upon v_ even in the face of considerably more rapid 
inner-shell dynamics. '+16 Interestingly, the frictional dynam- 
ics can be sensitive to the nature of the solvating medium, so 
that distinctly different solvent-dependent ET kinetics are 
anticipated depending on whether the activation energetics, or 
additionally the solvent dynamics, influences the reaction 
rate.'* This issue is considered further below. 


Solvent Barriers 


One of the earliest features of the theoretical treatments now 
collectively known as “Marcus theory” was a unique description 





Figure 2. 


Plot of Eop vs 1/d for the mixed-valence dimers 


[(bpy)eCiRu(L)RuCl(bpy)eP’*. where L = pyrazine, 4,4’ 
bipyridine, trans-1,2-bis(4-pyridyl)ethylene (data from ref. 17) 











of condensed-phase ET activation barriers as nonequilibrium 
solvent polarization barriers.'° A schematic representation 
of the barrier effect is shown below. A more quantitative 


representation followed from the recognition by Marcus that 
a transferring charge (electron) would essentially instantane- 
ously perturb the bound electrons associated with the sur- 
rounding medium (i.e. polar solvent), but would only slowly 


2 
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molecule 





Table ll. Spectroscopic, Structural and Reorganizational Parameters for Electron Transfer from Fe(CN)s* 


to Colloidal 





Relative 
__Intensity 
20.0 


> 


- A 
0.95 
6.61 0.33 
5.44 0.27 
0.27 a8 
1.00 0.59 
0.33 0.24 
1.12 0.89 
0.90 0.82 
0.56 0.69 
0.27 0.43 


Mode 
2118 cm-1 
2072 
2058 

720 

598 

540 

516 

484 

418 

364 





\Aal Ai 


Assignment 





0.048 A 
0.014 
0.026 


9 


0.026* 
0.039 


** 
*** 


*** 


0.059 


1000 cm’! Vc-n bridge 
340 Vc-N radial 
280 Vc-N terminal 
40 ? 
180 
60 
230 
200 
140 
80 


VFe-C 
VFe-c bridge 
VTi-O 
VTi-O 
VTi-O 
VTi-N 


* Taken from (or taken as) the crystallographically determined value for Fe(CN)¢** 
**Value not determined, since the measured normal coordinate displacement (A) may entail more than one type 
of bond length displacement (i.e., local-mode approximation may not be appropriate). 
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(in comparison to the timescale for an electronic transition) 
perturb the solvent nuclei. The relatively slow nuclear re- 
sponse was expected to create an intrinsic barrier to ET that 
could be surmounted by random solvent polarization fluctua- 
tions. The barrier was theoretically characterized, therefore, 
by both high frequency (€,) and low frequency (€,) dielectric 
components (corresponding to fast (electronic) and slow (nu- 
clear) polarization responses). Further consideration of the 
electrostatic consequences c: finite molecular reactant size 
(radius, r) and finite reactant separation distance or image 
distance (d) led to the following barrier expressions: 


AG" = (e7/4){(1/2r,) + (1/2rp) - 1/d)} (1/eg,, - 1/e,)(Ba) 
AG" = (¢/8)(I/t - 1/d)(1€/.,, - We.) (8b) 
where Eq. 8a is applicable to homogeneous solution-phase 
reactions, Eq. 8b to interfacial electrochemical processes, and 


e in both equations is the unit electronic charge. As suggested 
by Eq. 8, the polarization analysis offers some remarkably 





Figure 3. 


Intervalence transfer band energies vs 1kop - 1ks. Key to 
solvents: (1) nitrobenzene, (2) dimethyl sulfoxide, (3) 1-methyl- 
2-pyrrolidone, (4) hexamethylphosphoramide, (5) dimethy- 
lacetamide, (6) dimethylformamide, (7) formamide, (8) 
propylene carbonate, (9) acetone, (10) acetonitrile, and (11) 
deuterium oxide 
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simple predictions as to how barriers and, therefore rates, for 
electron transfer should depend on macroscopic solvent prop- 
erties, as well as molecular reactant size and geometry. Nev- 
ertheless, compelling experimental evidence in support of the 
solvent barrier theory was lacking for a number of years, 
despite several seemingly straightforward kinetics investiga- 
tions (such as rate studies as a function of the solvent, etc). In 
retrospect, many of these studies were bound to fail in this 
context because of the existence of other important solvent- 
variable kinetic factors, most notably: 1) electrostatic work 
terms (including interfacial adsorption and solution-phase ion 
pairing), and 2) dynamics effects (see below). In addition, at 
least some of the systems chosen for evaluation at electro- 
chemical interfaces (for example, reduction of protons to 
dihydrogen) were inappropriate because they involved exten- 
sive bond creation or annihilation. 

Given these circumstances, a key experimental devel- 
opment was the design and synthesis (largely by the 
Meyer'”:'* and Taube'*” groups) of molecular donor-ac- 
ceptor (D-A) systems featuring well-resolved charge- 
transfer absorption bands.2'?* When these systems are 
symmetrical (i.e. when the donor and acceptor sites are 
chemically identical as in Eq. 9), the absorption energy €,) 


Av 
(bpy), CIRu" — bridge — Ru Ci(bpy),** > 


(9) 


(nm QO) we GOO + OD) 


(bpy),CIRu'” — Ru"Cl(bpy),** 


for the lowest charge-transfer transition can be identified 
directly with the (vertical) Marcus reorganization energy, 
1.34 As already noted (Fig. 1), the reorganizational parame- 
ter, in turn, can easily be related to the expected thermal ET 
activation energy, AG”. For parabolic energy surfaces (i.e. 
those derived from displacement of harmonic oscillators) the 
energy relationship is simply: 


AG" = V4 (= E,,,/4) (10) 


Shortly after the appearance of the first valence-localized D-A 
studies,'*?! systematic variations in bridge length were exam- 
ined,'’ leading to a verification of the predicted dependence 
(Eg. 8a) of the solvent reorganization energy on D-A separa- 
tion distance d (Fig. 2). 
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A central aspect of the optical ET strategy as exempli- 
fied by the system illustrated in Fig. 3 (Eq. 11) is the 
evaluation of the reorganization energy (or 


bv 
(NH 5)sRa" XO)O} Ra” (NH 3) ——_—_» 
(NH 3)sRa™ NC ))><C hs a7 (NE 3) 5” (11) 


optical ET barrier) as a function of the solvent dielectric 
properties. The dependence found in this example,”>* and in 
several others, clearly is consistent with the qualitative predic- 
tions of the polarization theory (Eq. 8). We — along with the 
Brookhaven group”’ — have noted, however, that quantitative 
agreement is lacking: the extent of variation of A with solvent is 
(in this instance) markedly less than anticipated from theory. A 
Clue to the origin of this disparate behavior is provided by the ET 
distance-dependent study in Fig. 2. The nominal charge-transfer 
distance in reaction 11 (i.e. the ruthenium-to-ruthenium distance) 
is 11.3 A. The detailed chemistry, however, is such that the 
pertinent dx metal orbitals mix (selectively) with bridging ligand 
orbitals of the same symmetry (both 1 and n°) and become 
significantly polarized along the charge-transfer axis. The effec- 
tive ET distance, as recently determined by electronic Stark-effect 
spectroscopy,” is closer to 6A. A shorter charge-transfer distance 
should lead to correspondingly less polarization of solvent and 
smaller reorganization energies (again, note Fig. 2). With the 
experimentally determined distance modification, we find that 
agreement with the dielectric-continuum theory now becomes 
essentially quantitative. 

Related work on the solvent barrier/ET reactivity problem 
— in this case by the Brookhaven group” — has revealed the 
spectacular consequences of simple variations in molecular 
size. This group showed, for example, that systematically 
increasing the average molecular radii of pairs of simple 
transition-metal complexes (from 7A to 14 A) could in- 
crease bimolecular (homogeneous) electron-exchange rates 
(Eq. 12) by over five orders of magnitude!”? The observed 


M(L),°** + M(L),”* M(L),?* + M(L),** (12) 


reactivity variations are almost undoubtedly solvational in 
origin, as the redox pairs were carefully chosen so as to 
minimize reactant vibrational barriers (see below) as well as 
nonadiabaticity («,,) effects. Interpretation of the experiments 
in terms of nonequilibrium solvent-polarization phenomena 
provides a quantitative explanation for the reactivity pattern 
once the secondary influence of reactant size on separation 
distance (d; Eq. 2) is also taken into account. From a primitive 
physical perspective, the size dependence arises, then, from 
two sources: 1) a diminution of the effective charge density on 
the surfaces of the reacting molecular spheres (and therefore 
a decrease in surrounding solvent polarization) with increasing 
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sphere size, and 2) a diminution of partially compensating 
image-charge interactions as the center-to-center distance (d) 
for each reacting pair is likewise increased. 

Surprisingly, only recently have the corresponding effects 
of reactant size on electrochemical reactivity been explored. 
As Figure 4 illustrates, rates for 


M(L),** + e (electrode) <> M(L),”* (13) 


interfacial ET (Eq. 13) have been observed to increase expo- 
nentially with decreasing effective reactant radius.*° This ob- 
servation is broadly consistent, of course, with a Marcus-type 
solvent barrier effect. Two additional points, however, are 
worth noting. First, the magnitude of the kinetic size effect 
(spanning roughly two orders of magnitude in ET rate, and 
corresponding to ca. 3 kcal mo!" in barrier height) is consid- 
erably less than found in the study of homogeneous bimolecu- 
lar reactivity noted above. The most important difference 
chemically between the two studies is the involvement of only 
a single molecular reactant in each of the electrochemical 
processes, versus two reactants in the solution-phase proc- 
esses. At a crude intuitive level one would expect, therefore, 
roughly half the extent of solvent repolarization in the former 
as the latter. Indeed, this notion — expressed more rigorously 
by Marcus” — underlies the so-called homogeneous/hetero- 


Figure 4. 





Log of ks for the indicated ET reaction (eq. 13; normalized to ks for 
Ru(NH3)63+/2+) versus the inverse radius of the reactant. The elec- 
trode material is low-defect-density HOPG; the electrolyte is 1 M aq 
KCl. Line drawn is a best-fit line for all points except point 6 (open 
circle). Key to data points: (1) Ru(NH3)63+/2+,(2) Ru(NH3)5(pyri- 
dine)3+/2+, (3) (NH3)5Ru(pyrazine)Ru(NH3)55+/4+ (4) 
(NH3)5Ru(pyrazine)Ru(NH3)56 +/5+ ,(5)t-(pyridine)(NH3)4 Ru(pyrazi 

ne)Ru(NH3)4(pyrazine)5+/4+,(6)t-(pyridine) (NH3)4Ru(pyrazine) 
Ru(NH3)4(pyrazine)6+/5+, and (7) Fe(phenanthroline)33+/2+ 
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geneous cross relationship (see Eq. 3) that permits solution- 
phase reactivities to be predicted from electrochemical kinet- 
ics data and vice versa. The second point is that two of the 
reactants in the electrochemical study are actually small di- 
ameter reactants whose effective sizes have been increased by 
covalently linking them to an additional reactant and then 
delocalizing the associated charges over the resulting “super 
molecule”. We suggest that this could prove to be a more 
generally effective strategy for accelerating outer-sphere elec- 
tron-transfer reactivity. In any case, it provides a simple illus- 
tration of the experimental exploitation of basic barrier 
concepts to manipulate reactivity. 

As perhaps suggested by the preceding discussion, one of 
the active areas of current research is the synthetic manipula- 
tion of specific chemical systems to provide detailed control 
of solvent-based barrier effects. Several reports from Curtis 
and co-workers have illustrated how control can be exerted via 
a Clever combination of mixed solvation and selective reac- 
tan/solvent hydrogen bond formation.*'*” Related studies at 
Northwestern have illustrated how local solvation can be 
manipulated to: a) trigger kinetically observable intramolecu- 
lar ET events,” b) create barriers and induce valence localiza- 
tion (or partial localization) in otherwise delocalized, 
multi-site systems,” and c) significantly modulate superex- 
change-based electronic coupling effects.*° 

Another important focus of contemporary research is 
solvent “molecularity”. The nonequilibrium solvent polariza- 
tion theory of Marcus achieves its simplicity and broad appli- 
cability, in part, by representing the solvent as a structureless, 
polarizable medium, i.e. “molecularity” is intentionally ne- 
glected. Nevertheless, molecules — including solvent mole- 
cules — are obviously of intrinsic interest to chemists! Recent 
attempts to gain a more molecular-level understanding of 
solvent reorganization effects have proceeded along at least 
four different lines. One involves the introduction of a mean 
spherical approximation (MSA) for individual solvent mole- 
cules, within the context of a broader electrostatic repre- 
sentation of medium repolarization effects. An interesting 
MSA-based model due to Wolynes treated the nonequilibrium 
solvation problem in terms of a frequency-dependent Gurney 
cosphere engendered by the solvent molecular size.** This 
approach, along with other MSA treatments, usually yields 
reorganization energies that are somewhat (say 15-25%) 
smaller than obtained from the conventional continuum esti- 
mates. A numerical comparison with solvent-dependent E 
values for biferrocene species indicates a rough concordance 
in several aprotic solvents, although not in hydrogen-bound 
media.*’ A more recent MSA treatment by Blum and Fawcett 
for equilibrium solvation (i.e. the zero-frequency part of the 
solvent reorganization energy) emphasizes anticipated differ- 
ences between the solvation of cations and anions.** This 
model has also been applied to an analysis of the solvent-de- 
pendent electrode kinetics of Cp,Co*” versus Cp,Co™ (Cp= 


cyclopentadieny!).** Interestingly, the MSA model can ac- 
count for the observed more sluggish kinetics of the latter 
couple, from the predicted higher solvent activation barrier, 
although other factors (inner-shell barriers, specific double- 
layer effects, ion-pairing, solvation dynamics, etc.) perhaps 
also play a role. The evaluation of E_ data for differently 
charged mixed-valence systems in this context would be 
worthwhile. 

An alternative, perhaps complementary, semi-emperical 
approach to estimating AG” is the so-called “nonlocal” elec- 
tronic treatment due to Kornyshev and Ulstrup,“° which con- 
siders the effects of spatial correlation on the 
wavevector-dependent solvent dielectric properties. One 
worthwhile facet of this model is that it can provide a viable 
rationalization of some deviations from the continuum predic- 
tions observed for metallocenes in hydrogen-bound me- 
dia.?’*! 

Asecond, fundamentally different, approach to molecular 
solvent reorganization makes use of molecular-dynamics 
(MD) simulations and carefully chosen potentials for solvent 
interactions. Some encouraging results based on this strategy 
have emerged from Chandler and coworkers” (see also the 
article in this issue by G. Voth). This group has successfully 
simulated the redox activation behavior of the well-known 
hexaaquo iron(III/II) system in water as solvent, and has 
derived important new insights regarding the role of nuclear 
tunneling.*” Most notably, they found a remarkable concor- 
dance with the Marcus prediction of parabolic potential-en- 
ergy surfaces for solvent reorganization, even though complex 
short-range interactions abound in this highly charged/hydro- 
gen-bonded system. On the other hand, model MD simulations 
for both dipole creation/annihilation® and ET“ in methanol 
show clear evidence for non-linear polarization effects, that is, 
significant deviations from parabolic free-energy surfaces. 
The latter simulations also yield higher activation barriers for 
anion-neutral than for cation-neutral reactant pairs (vide su- 
pra) attributable to short-range solvation, although the effect 
is small for reactant radii approaching those prevalent in ET 
phenomena.“ Overall, while such MD simulations can be 
complex and time consuming, they are beginning to add a new 
molecular-level dimension to our understanding of the ener- 
getics of nonequilibrium solvation. The associated dynamical 
aspects are discussed briefly below. 

A third approach is largely experimental and focuses on 
the controlled assembly of orientationally constrained, local 
solvent environments (for example, large “encapsulating” 
crown ethers, semi-rigid hemicarcerands, etc.). Such assem- 
blies may be amenable to detailed in situ (multidimensional 
NMR) and/or ex situ (crystallographic) structural charac- 
terization as a function of encapsulated reactant or product 
oxidation state. Initial studies have shown that the overall 
solvent reorganizational energetics for such systems can read- 
ily be monitored by optical ET methods (cf. Eq. 12).*° Finally, 
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Figure 5. 


Plot of logarithm of work-corrected standard electrochemical 
rate constants kex° at mercury-aqueous interface versus cal- 
culated inner-shell activation barrier AGis . Data taken from 
J.T. Hupp and M.J. Weaver, Inorg. Chem, 22, 2557 (1983); T. 
Gennett and M.J. Weaver, Anal. Chem., 56, 1444 (1984) 
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a fourth approach — pioneered by Hendrickson and co-workers 
— involves solid-state electron transfer. Single crystals of 
various molecular mixed-valency species are prepared with 
zero, one, Or two solvent molecules of co-crystallization. 
Coupling of isolated solvent motions to ET can then be exam- 
ined via a combination of X-ray and variable-temperature 
Mossbaur spectroscopy techniques. 

The presence of such a myriad of molecular solvent 
reorganization treatments notwithstanding, it is worth reiter- 
ating that the Marcus dielectric-continuum approach remains 
a generally useful, and in most cases, a semiquantitatively (or 
even quantitatively) reliable means of estimating total ET 
solvent reorganization energies. That this should be the case 
even though the relevant physics contains a common element 
with the Born solvation model — well known to be much less 
reliable — lies partly in the usual dominance (at least in polar 
media) of the optical frequency component (e_ *') with respect 
to the Born term (¢, ') in Eq. 8. Indeed, quite marked deviations 
between the equilibrium solvation energies and the Born pre- 
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dictions can occur without influencing substantially the non- 
equilibrium solvation energy ( A or AG"). We discussed this 
point some time ago in an analysis that harnesses experimental 
solvation energies for redox couples (i.e. differential solvation 
energies) in place of the Born term in Eq. 8.*’ 

Other solvent related issues of current interest and impor- 
tance include: 1) “solvent” reorganization in protein environ- 
ments (where again, additional molecular-level theory 
development would be invaluable), 2) solvent reorganization 
at nonmetal/solution interfaces — for example, liquid/liquid 
interfaces, semiconductor/solution interfaces, and mem- 
brane/solution interfaces, and 3) solvent reorganization at 
mesoscopic and nanoscopic particle(or electrode)/solution in- 
te1faces. Some insightful theoretical work on liquid/liquid and 
semiconductor/liquid interfaces has appeared recently from 
Marcus and from Smith and Koval.*? Experiments capable 
of testing these theories would be most useful. Finally, in the 
mesoscopic chemistry field, a number of exciting experimen- 
tal developments relating to interfacial ET and energy conver- 
sion have emerged in the last three years, in particular.°°* 
Fundamental studies of barrier phenomena here could provide 
a rational basis for further advances and could help to explain 
some of the more provocative existing results. 


Vibrational Barriers 


A second major component of most ET reaction barriers 
is vibrational reorganization. This component arises because 
of the oxidation-state dependence of the normal coordinates 
or internal bond lengths. Interconversion of oxidation states 
(electron transfer) therefore requires the displacement of co- 
ordinates (bond compression or bond stretching) and is accom- 
plished, in the Marcus classical limit, by vibrational 
activation.*>*? The combination of net coordinate displace- 
ment and transient vibrational excitation leads to an energy 
barrier shaped much like the one in Figure 1. In this picture, 
the actual charge transfer occurs at the top of the barrier (the 
transition state) where the best compromise, in terms of bond 
lengths, has been achieved between the redox reactant and 
product. According to Marcus, the vibrational activation bar- 
rier (AG* .) or reorganization energy (A,..) for an electron 
exchange reaction (homogeneous or electrochemical) can be 
estimated simply from known bond length changes (Aa) and 
associated force constants (f): 


Nr = WE Aa (14 
J 


In Eq. 14, b is the number of equivalent bonds displaced, and 
the summation (j) is over all bonds displaced. 

The importance of vibrational reorganizational effects 
can hardly be overstated. For example, Brunschwig et al, in 
a study of about a dozen transition-metal redox couples, found 
that bimolecular self-exchange rate variations over some /5 





orders of magnitude could be explained almost entirely on the 
basis of classical vibrational activation effects. We showed in 
a related study that similar effects can also control reactivity 
for cross reactions and at electrochemical interfaces.** Figure 
5, for example, shows a plot of standard double-layer corrected 
electrochemical rate constants for four transition-metal com- 
plexes of similar size versus independently estimated vibra- 
tional activation energies. As in the Brunschwig study, a 
reasonable correlation exists. 

In both the examples above, vibrational barriers were 
derived from X-ray measurements (EXAFS or crystal- 
lography) of metal-ligand bond lengths in oxidized and re- 
duced states, with the bond length differences then 
incorporated into Eq. 14. This strategy is useful when: a) 
normal coordinate displacements can be approximated by 
local coordinate (i.e. bond length) displacements, b) the num- 
ber of coordinates or types of bonds displaced is small, c) both 
redox states are chemically stable and long lived, and d) the 
available crystalline or solution X-ray measurement environ- 
ment is sufficiently similar to the redox reaction environment 
to yield kinetically meaningful structural data. In many cases, 
however, one or more of these conditions is not fulfilled, and 
alternative routes to structural data are required. Viable routes 
include: 1) Franck-Condon analysis of structured emission 
(fluorescence or phosphorescence) from photo excited elec- 
tronic states,°® 2) empirical correlations (e.g. Badger’s rules, 
etc.) of redox-induced vibrational frequency shifts with nor- 
mal or local coordinate displacements,”’ and 3) time-depend- 
ent analysis of resonance Raman scattering intensities.~* 

Over the past few years the Northwestern group’*™, and 
the Myers group at Rochester,” have made extensive use of 
the third strategy and have found it to be extremely powerful. 
To implement the strategy, it is necessary to identify electronic 
absorption bands corresponding to the thermal electron trans- 
fer reaction of interest. This is reasonably straightforward for 
homogeneous solution-phase reactions. We find, however,°! 
that the Raman method can also be applied to interfacial 
reactions such as reaction 15: 


Fe(CN),* Fe&(CN),* 


Fe(CN),” Fe(CN), 


» + 
FelC™), Fe(CN),* 


»n * 
Fe(CN), Fe(CN),* 


Visible region excitation of an intense ferrocyanide to titanium 
dioxide (particle or electrode) charge transfer band leads, in 
this case, to enhanced Raman scattering from nearly a dozen 
vibrational modes. Observation of enhancement implies vi- 
bronic coupling of these modes to the transition of interest, i.e. 
Eq. 15 or the reverse. Time-dependent theory then permits one 
to connect the observed scattering intensities quantitatively to 
normal coordinate displacements (A) — which, in tum, may be 


transformed into bond length changes. Table II summarizes 
the results of such an experiment. A number of points are worth 
noting. First, the symmetry of the molecular system is clearly 
reduced by surface binding, resulting in activation of several 
vibrational modes which are completely undetectable in solu- 
tion or single crystal environments(i.c. X-ray measurement 
environmenis). Second, the electrode (i.e. Ti0,) lattice itself 
is activated by interfacial ET. Third, the overall vibrational 
barrier can be partitioned into high- and low-frequency com- 
ponents. This then permits important tunneling corrections to 
the classical Marcus vibrational analysis to be applied in a 
mode-specific fashion. Finally it is worth noting that this 
particular system would be very difficult to examine by con- 
ventional structural methods because the redox form on the 
right hand side of Eq. 15 returns to the initial form in less than 
400 ns.“ The Raman method overcomes this problem by 
interacting with the short-lived state only in a resonance 
fashion; significant real lifetimes, therefore, are not required. 

It is appropriate to conclude this section by outlining some 
of the current issues and problems in this area of redox 
chemistry. One of these concerns bridged systems where cou- 
pling of bridge-localized vibrations to intramolecular ET is of 
appreciable theoretical, as well as experimental, interest. Also 
currently being investigated — primarily by resonance Raman 
in the extended near infrared — is the question of how specific 
vibrational modes conspire to cause electronic localization 
versus delocalization in more strongly electronically coupled 
systems. Conceptually related to the Raman studies, and 
especially the time-dependent scattering studies, are new theo- 
retical studies involving wave packet propagation methods. 
These studies, being carried out primarily by Nitzan and 
Ratner™ (see accompanying article), may offer significant 
insights into the real-time progress of ultrafast ET reactions. 
The increasing experimental availability of ultrafast tech- 
niques (i.c. laser-based techniques which can probe reactions 
in the 15 to 50 femtosecond time regime) should soon permit 
some very exciting issues in rapid-reaction chemistry to be 
examined. One of these concerns coherence effects, where 
extremely fast reverse ET might, in fact, provide a fundamen- 
tal dephasing mechanism for rapidly created states. On a 
slightly longer timescale (i.e. 0.5 to 5 picoseconds) the avail- 
ability of tunable transient infrared spectrometers may permit 
one to understand how ET product states dispose of excess 
vibrational energy and reconfigure themselves for further 
reaction. Indeed, preliminary studies along these lines (by 
Doorn, Dyer, Woodruff and Stoutland at Los Alamos”) appear 
promising. Transient IR methods should also permit chemists 
to examine how photochemically generated, vibrationally hot 
systems behave in redox transformations. This problem could 
be of particular importance in molecule-based studies of solar 
energy conversion. Emerging work from Spears and co-work- 
ers™ suggests that highly unusual reactivity patterns will be 
seen. Other solution-phase work of interest concerns reactions 
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where the formation of one or more bonds directly accompa- 
nies the transfer of a single electron. A specific example, from 
the Barbara group, involves charge transfer between a halogen 
atom and a dimeric (van der Waals) benzene cation.” Here 
two bonds appear to be formed. Interestingly, adequate under- 
standing of the problem — in both a dynamic and a vibronic 
sense — would appear to lie beyond the realm of current theory. 

Returning to surfaces or interfaces, the unexpected obser- 
vation of electrode lattice activation during electron transfer 
(Eq. 15) merits further investigation. Is the phenomenon re- 
stricted to semiconductor materials, or might it also be signifi- 
cant at metal surfaces? For either class of material or surface, 
are there vibrational activation consequences for quantum 
confinement, i.e. for geometrically constraining the trans- 
ferred electron or hole within an electrode or particle of 
smaller size than the characteristic e or h* radius? Both 
questions could be of importance in emerging redox-based 
applications of colloidal particle dispersions. We further sug- 
gest that both may be amenable to investigation by systematic 
electrochemical, spectral and theoretical studies of cluster 
systems. Finally, one may ask whether there are likely to be 
unusual local, or collective, vibrational effects in reactions 
involving self-assembled monolayers. Preliminary studies by 
Abruna and co-workers, based on an unusual X-ray standing 
wave method, indicate that there are.’ In the X-ray study clear 
evidence is found for coupling of an “accordion” type bridge 
compression motion to ET between an electrode and a tethered 
transition-metal reactant. Also of interest is a very recent 
observation by Bedzyk and Mirkin of long range two-dimen- 
sional ordering in a redox-active monolayer (again via stand- 
ing-wave techniques).’' These new experiments should open 
up the possibility of detecting and evaluating collective coor- 
dinate displacement effects in interfacial ET processes. Inter- 
estingly, collective — or at least correlated — coordinate 
displacements appear to play a very important role in the 
kinetics of ET in solid-state environments (i.e. three-dimen- 
sional crystalline environments), as demonstrated by Hen- 
drickson and co-workers.*° 


Driving Force Effects 


Just as important (in a kinetic sense) as solvational and 
vibrational reorganization energies are the effects of thermo- 
dynamic driving forces (reaction free energies). As suggested 
by eqs. | and 2, the way in which driving forces affect ET rates 
is primarily through classical barrier height manipulation (al- 
though other effects, such as electronic tuning, are known). 
For small changes in driving force (AG*) a linear variation in 
AG" is predicted by Marcus theory; for larger changes the 
expected variation is quadratic. In other words, an exponential, 
or nearly exponential, increase in ET rate with increasingly 
favorable driving force is predicted — much as in other transi- 
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tion-state theories. A very large number of examples of this 
type of behavior have been described experimentally. 

More unusual is the prediction, already noted, that for 
very large driving forces (-AG®° > A), electron-transfer rates 
should decrease with increasing exoergicity, i.e. “inverted” 
rate versus driving force behavior should be seen. This well- 
known prediction proved quite resistant to verification, lead- 
ing many to conclude that the underlying theory was not fully 
correct. In 1979, however, inverted rate behavior — consistent 
with Marcus theory — was finally observed by Beitz and Miller 
for a series of bimolecular ET reactions in low-temperature 
glasses.”* Nevertheless, because of the exotic medium, many 
researchers were skeptical. Ultimately, in 1984 — nearly 30 
years after the initial theoretical proposal — Miller, Calcaterra, 
and Closs provided an even more convincing demonstration 
of rate inversion.”> These experiments involved intramolecu- 
lar ET and were based on a clever application of pulsed 
radiolysis methods to solution-phase kinetics. In retrospect, 
the functional equivalent of inverted ET was, in fact, seen 
some time before in various studies of nonradiative decay with 
organic” (and later inorganic’®) chromophores as both Mar- 
cus’°and Meyer” have noted. It is now generally appreciated 
that the so-called “energy-gap law”, often used to characterize 
these experiments, can be obtained from a low-temperature 
quantum version of Marcus’ classical inverted-region analy- 
sis. 

Since 1984, more than two dozen additional examples 
have appeared, encompassing a variety of chemical contexts. 
For illustrative purposes we will focus on two of the more 
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recent. The first is represented schematically by the photo-re- 
dox sequence shown below (where Ir, is [Ir(1,5-cyclooctadi- 
ene)(\-pyrazolyl)}, and py” is any of several alkyl pyridinium 
ions): 


tr, Sir (16a) 


Ir + py*— Ir3*+ py” (16b) 


Ir, ** + py*—> Irs" + py* (16c) 
Note that ET is bimolecular and that it can occur in two 
directions. Furthermore, Gray and co-workers have shown 
that the driving force can be made to exceed A in both.” 
Surprisingly, however, rate inversion occurs only in the re- 
verse direction (i.e. diiridium cation reduction, eq. 16c). In the 
forward direction, rates increase to a diffusion-limited value 
with increasing -AG®, but are then unresponsive to further 
increases in driving force. Several explanations have been 
devised. But, as yet, none have been confirmed. The second 
example is shown in Figure 6 and comes from Northwestern.” 
The reaction involved is Eq. 15, where the driving force has 
been varied by replacing one of the six available CN” ligands 
of the parent compound with any one of several pyridyl 
ligands. The absolute driving forces are unknown at present, 
but clearly they will scale as the formai potentials (E,) of the 
molecular reactants in solution. In this case it also has been 
possible to calculate the slope of the rate vs. driving-force plot 
by using Raman-derived displacement and force-constant data 
(Table IT) and low-temperature theory. The calculated slope 
(roughly a factor of ten decrease in rate for a 300 mV increase 
in driving force) is in good agreement with experiment. Per- 
haps surprisingly, the experiment represents the first in which 
rate inversion has been observed at an interface. 

Among other important current issues and problems re- 
lating to driving-force phenomena are: 1) “gating” effects (see 
article by Ratner) which may decouple ET from very slow 
nuclear coordinates and make certain ET processes driving- 
force independent,’** 2) exploitation of rate inversion to 
inhibit reverse ET — and enhance efficiency — in artificial, 
molecule-based solar energy conversion devices,® and 3) 
adaptation of Marcus-type analyses to atom and ion transfer 


reactions.***5 


Reaction Dynamics 


We have so far focussed attention primarily on the ener- 
getic aspects of ET chemistry, that is the factors and strategies 
utilized to understand and evaluate energy barriers. A com- 
plete description of ET kinetics, however, also necessarily 
involves evaluation of the reaction dynamics, that is, the 
rapidity by which the system is able to move along the reaction 


coordinate so to yield products from reactants. As already 
noted, several factors may influence the reaction dynamics, 
and hence the preexponential factor for activated ET proc- 
esses, not the least of which is the electron-tunneling prob- 
ability for nonadiabatic processes. Much experimental 
information on the kinetic consequences of electronic cou- 
pling in nonadiabatic ET has been gained from intramolecular 
D-A systems where the donor and acceptor sites are separated 
by extended and/or saturated organic linkages. These include 
deliberately synthesized systems,’> as well as modified bio- 
logical systems such as metalloproteins,*’ both of which also 
have been featured prominently in tests of driving-force ef- 
fects (see above). In the context of electronic coupling effects, 
primary tactics involve discerning the kinetic response to 
alterations in the spatial and/or structural D-A separation. 
Besides homogeneous-phase intramolecular ET, substantial 
advances have been made along these lines for structurally 
well-defined films at electrochemical surfaces.**” 

We confine our comments here primarily to the reaction 
dynamics anticipated for adiabatic processes, where nuclear 
(especially solvational) rather than electronic factors should 
prevail. Up to the late 1970's, theoretical treatments of solvent 
effects, including Marcus theory, focussed on activation ener- 
getic aspects. Beginning around 1980, and continuing to the 
present, substantial research effort has been invested in under- 
standing the dynamical roles of reactant solvation.2’” Pro- 
gress has been driven, in particular, by the recognition of 
“solvent-friction” effects in chemical kinetics (which, as men- 
tioned above, constitute an interesting deviation from the usual 
TST approach). While the precise physical meaning and inter- 
pretation of solvent-friction effects depend somewhat on the 
type of reaction being considered, the usual notion is that an 
impediment exists to net progress along the reaction coordi- 
nate, because of irreversible energy dissipation from the “re- 
active” (possibly solvent molecular rotation) modes to the 
surrounding solvent “bath”. The reaction coordinate motion is 
then no longer smooth and unidirectional (as in TST) but is 
characterized instead by “fits and starts”, commonly termed 
“overdamped motion”, perhaps even involving several re- 
crossings within the barrier-top region before the reaction is 
finally consummated. 

While ET was historically not the first class of reactions 
for which the notion of solvent friction was considered in 
dynamical descriptions, the central role of solvation in the ET 
reaction coordinate attracted attention in this context by theo- 
reticians by the early 1980’s.°'*? Since then, the literature has 
expanded enormously. Most of this work emphasizes the idea 
that “dielectric friction”, associated physically with collective 
rotations (and related motions) of solvent dipoles that neces- 
sary attend solution ET processes, should exert significant or 
even substantial influences upon the net barrier-crossing dy- 
namics, and hence the ET reaction rate. 
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The emergence of this theoretical framework has inspired 
a large number of experimental studies aimed variously at 
probing the nature of, and the extent to which, collective 
solvent polarization can influence the dynamics of ET and related 
processes. Broadly speaking, the experimental inquiries fall into 
four categories. The first (type I) involves the evaluation of 
time-dependent fluorescence Stokes shifts (TDFS) for solute 
chromophores forming suitable charge-transfer excited 
states.'®'°' In optimal cases, these ultrafast (> 50 ps) laser-in- 
duced measurements can probe directly the real-time dynamics 
of polar solvent relaxation around a newly formed solute dipole, 
and have yielded some fascinating information of relevance to 
ET kinetics. The second type concerns a related use of ultrafast 
pulse lasers to follow intramolecular ET processes emanating 
from phototexcited states.'®' So far, however, almost all of the 





Figure 7. 


Logarithmic plots of “barrier-corrected” rate constants (M 
s', extracted from rate and optical barrier data) versus inverse 
longitudinal relaxation time (s') for five metallocene self-ex- 
change reactions in 11 solvents, taken from ref. 29. See Table | 
of ref. 107 for key to solvents. Key to redox couples: filled circles, 
Cpz Co*/0 (Cp = pentamethyicyclopentadieny/); filled squares, 
Cp2°Co*” [Cp® = (carboxymethyl) cyclopentadieny\]; filled trian- 
gles, CpeCo*”: open triangles, CpeFe*”: open squares, (hy- 
droxylmethyl) ferroceniuny-ferrocene 
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ET reactions studied in this fashion have evolved only small 
free-energy barriers, AG” > k,T (k, is Boltzmann's con- 
stant).''-' (The additional role of inertial solvent motions, 
therefore, can probably not be overlooked.) 

The remaining two experimental categories constitute 
attempts to extract information conceming dynamical solvent 
effects for activated ET processes (say, for AG” > 5k,T). The 
most common approach (type IIT) involves the acquisition of 
solvent-dependent kinetic data for bimolecular outer-sphere 
electron-exchange reactions, especially symmetrical self-ex- 
change, and for related reactions at electrode surfaces. '* The 
kinetics of the former are most commonly evaluated by means 
of magnetic resonance (nmr, epr) line-broadening methods. 
These reaction systems feature significant (commonly 5-10 
k,T) solvent reorganization barriers (i.¢., they truly are acti- 
vated processes). In addition, these systems provide opportu- 
nities to examine the dynamical consequences of 
systematically altering solute-solvent interactions, reactant 
electronic structure, vibrational barriers, etc. Despite the clear 
fundamental importance of activated ET processes, however, 
the experimental extraction of information on barrier-crossing 
dynamics is somewhat problematic. The core of the problem 
concerns the need to separate the dynamical (preexponential 
factor) and energetic (barrier height, AG") contributions to the 
measured solvent-dependent rate constants. In most cases, the 
extraction of the desired dynamical information from rate data 
requires one to use theoretical estimates for AG”. An additional 
problem in these conventional analyses is that the geometry as 
well as stability of the precursor complexes for bimolecular 
(or electrochemical) outer-sphere processes is often unknown. 

The conclusions regarding the role of solvent dynamics 
reached in this fashion are therefore often critically dependent 
on the validity of the theoretical models utilized to separate the 
dynamical and energetic factors. These type III tactics have 
consequently tended to be regarded as a “poor cousin” in com- 
parison with the more direct (and sophisticated-looking) insight 
into real-time dynamics obtained from approaches I and II. Most 
of the uncertainties faced with such dynamical analyses for 
activated ET reactions could be circumvented by evaluating 
solvent-dependent rate data for intramolecular (preferably sym- 
metric) ET reactions (labelled here as type I'V). The evaluation of 
both the unimolecular ET rates and the barrier heights (the latter 
from optical ET energies) would provide the preferred route to 
the separation of dynamical and energetic factors. Unfortunately, 
however, relatively few type IV systems have been scrutinized in 
the context of solvent dynamics.'™ 

Despite the complications, sufficient information on bar- 
rier height and related factors can be obtained for certain 
bimolecular ET reactions to allow reliable, and even quantita- 
tive, deductions to be made regarding solvent-dependent bar- 
rier-crossing dynamics. The specifics have been reviewed in 
some detail recently.'?-" Broadly speaking, the extant experi- 
mental information lends support to the earlier theoretical 





predictions that overdamped solvent motion (“solvent fric- 
tion”) can often influence and even control the ET barrier- 
crossing dynamics, as prescribed most simply by the 
longitudinal solvent relaxation time, t, . While a large number 
of type III reactions (both homogeneous and electrochemical) 
have now been examined, one series of systems studied by the 
Purdue group a few years ago provides an interesting example 
of what may be discerned experimentally in optimal cases. The 
series consists of metallocenium-metallocene self-exchanges: 


Cp,M* + Cp)M ~ Cp)M + Cp)M” (17) 


where Cp is a cyclopentadienyl ligand (or derivative), and M 
is iron or cobalt.'*'*!® The reactions were selected for 
several reasons, including the dominance of outer-shell (sol- 
vent) reorganization in the activation energetics, and the op- 
portunity to modulate the relevant electronic properties via 
both metal and ligand alterations. 

Figure 7 shows a logarithmic plot of the self-exchange 
rate constant kK (M ~'s™') for five metallocene couples ver- 
sus the inverse longitudinal relaxation time for a range of 
solvents (see caption and ref. 107 for details). The rate con- 
stants k’ were corrected in an approximate way for solvent- 
induced variations in barrier height via comparisons to optical 
ET data for a related mixed-valence ferrocenium-ferrocene 
system; solvent-induced variations in K_ for a pair of given 
self-exchange reactions, therefore, should reflect only differ- 
ences in the reaction dynamics (i.e. Kv, in Eq. 5). The +" 
parameter, extracted from dielectric loss spectra, provides an 
approximate assay of the overdamped solvent dynamics in 
so-called “Debye” media (those characterized by only a single 
relaxation time). While the solvents selected for Fig. 7 gener- 
ally exhibit only approximate Debye behavior, the solvent 
friction dynamics (i.e. the .- values) are seen to span a 
substantial numerical range (ca. 1 x 10''s' t0 4 x 10's"). 

Interestingly, the sensitivity of the ET rates to the solva- 
tion dynamics varies systematically with the nature of the 
metallocene: The slowest reactions (ferrocenium/ferrocene) 
are insensitive to the solvent dynamics, while the fastest 
(cobaltocenium/cobaltocene) vary almost proportionally with 
1" — with intermediate behavior seen for others. This spec- 
trum of reactivity can be rationalized on the basis of a transi- 
tion from nonadiabatic to adiabatic behavior (see Background 
section above). The slowest reactions follow largely nonadi- 
abatic ET pathways, so that the barrier-crossing probability is 
essentially proportional to (H,,)* and independent of the solva- 
tion dynamics. For reactions featuring larger D-A electronic 
coupling (larger H,,), reaction adiabaticity is approached or 
reached, and the reaction dynamics become highly sensitive 
to the nuclear dynamics (in this case, solvent friction). 

The observation ultimately of almost pure solvent dy- 
namical control (cobaltocene reactions; see figure) also allows 
approximate estimates of H,,to be obtained for the entire series 
of metallocenes.'”’ The resulting H,, estimates, varying from 


about 0.075 to 1 kcal mol", are in reasonable agreement with 
theoretical predictions by Newton and coworkers,'!'"!? and 
can be understood qualitatively in terms of the primarily 
metal- and ligand-centered nature of the redox orbitals in the 
ferrocene and cobaltocene couples, respectively. The observa- 
tion of clearly solvent-controlled adiabatic ET rates also en- 
ables the detailed nature of the solvation dynamics for 
activated ET processes to be explored. One aspect of this latter 
endeavor concems the extent to which faster solvent dynami- 
cal components can accelerate the reaction rate beyond that 
prescribed by more overdamped motions. Comparisons be- 
tween such kinetic data and subpicosecond TDFS measure- 
ments in non-Debye media (such as methanol) have 
demonstrated the importance of such effects.'*!'” 

Prompted by theoretical work by Marcus and cowork- 
ers,'*5 an important issue also receiving attention is the extent 
to which solvent dynamics can limit ET reaction rates even in 
the face of faster (and inherently underdamped) reactant vi- 
brational components of the activation barrier. Experimental 
data has so far provided at least qualitative support to the 
theoretical predictions.'*"'® 


Concluding Remarks 


Although the foregoing is only an incomplete (and sub- 
jectively selective) survey of some experimental aspects of 
fundamental research in solution-phase ET processes, it is 
evident that much has been learned concerning the physical 
and chemical factors that control activation energetics and 
dynamics, especially for small molecule reactions. In particu- 
lar, the interplay between theory and experiment, initiated 
almost forty years ago by Marcus, is now stronger and more 
diverse than ever. This situation reflects not only the present 
buoyancy and health of the subject, but bodes well for the rapid 
development of our understanding of ET processes in much more 
complex systems. The future promises to be an exciting one! 
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Electron Transfer Rate 
Theory, and Control of 
Electron Transfer 


Processes 


Mark A. Ratner 
Northwestern University 


Introduction 


The process of electron transfer is one of the most impor- 
tant and fundamental in all of chemistry’’. It occurs over a 
very broad range of time scales (from subpicoseconds to many 
seconds), and in a bewildering variety of natural and synthetic 
systems ranging from bacteria to biosensors, from photosyn- 
thesis to photomultipliers. The fundamental understanding of 
rates of the electron transfer process, and what microscopic 
parameters and behaviors control the rate of electron transfer, 
has been the focus of Rudy Marcus’ work for nearly forty 
years, and is the basis for his having been awarded the 1992 
Nobel Prize in Chemistry. Despite Rudy’s enormous contribu- 
tions,'°'? and important early work on electron transfer rates 
by Noel Hush,” and by Levich and his school,2"™ and a 
number of other workers, many important issues of under- 
standing electron transfer reactions still remain to be clarified. 

One of the great strengths of chemistry is that, based on 
increased understanding, chemists can create and design new 
molecules, new structures, and new processes. In electron 
transfer, we are only beginning to see designed molecular 
systems whose electron transfer rates can be controlled. These 
designed electron transfer systems, one component of an area 
often referred to as molecular electronics,”>*° may lead to 
important advances in energy transduction, in information 


storage and retrieval, in chemical sensing and in integrated 
molecular electronic devices. 

In this article, I will discuss the understandings of electron 
transfer rates that have been developed in the past four dec- 
ades, and the beginnings of molecular system design to control 
electron transfer rates and processes. 


The Beginnings: Activated 
Complex Theory and Solvent 
Polarization 


Modern experimental studies of electron transfer began 
after World War II, when the availability of separated isotopes 
made it possible to measure experimentally the rates of elec- 
tron transfer exchange reactions between metal complexes.”' As 
such measurements became available, chemists searched for a 
simple rate theory to describe electron transfer. Generically, one 
can write an electron transfer reaction as in equation 1, 


A+D>~A+D* (1) 
where A and D denote acceptor and donor, respectively. In fact, 
as Marcus realized, for separated species in solution the rate 
of electron transfer will depend upon the distance and orien- 
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tations of the donor and acceptor species. For this reason, it is 
often simpler to study intra-molecular electron transfer reac- 
tions, whose experimental investigation was pioneered by 
Taube and his students. A characteristic intra-molecular trans- 
fer system*? is that indicated by equation 2, 


Ru* e Co* 
1 (2) 


Ru’* - Co** 


in which electrons move from the stable Ru”* to form the labile 
Co**. In such intra-molecular transfer reactions, there is no 
work term to assemble the reactants into a configuration 
appropriate for electron transfer. 

Marcus’ original formulation of the electron transfer rate 
problem was based on activated complex theory. That is, he 
wrote the electron transfer rate as: 


k., = Zexp [-E, / RT] (3) 


a typical Arrhenius expression. Here the prefactor Z is the 
molecular collision rate, generally taken as 10'7/sec, while E, 
is the activation energy for electron transfer. Marcus’ critical 
realization was that when electrons move, the changing mul- 
tipole charge distribution is coupled to polarization fluctua- 
tions in the solvent. Indeed, the activated complex state is that 
in which the polarization for the reactant and the product are 
the same. This is shown schematically in figures 1 and 2, which 
show the reactant well, the product well, the activated com- 
plex, and the energy quantities AE° = A (reaction exoergicity) 





Figure 1. 


Ground state potential energy surface, in one schematic dimen- 
sion Q, for electron transfer. The activation barrier and overall 
exoergicity, AE, are defined. The minima corresponding to reac- 
tant and product are denoted Q; and Q, respectively. 
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and A (reorganization energy). Marcus realized, in 1956, that 
the polarization energy in both reactant and product states is a 
quadratic function of the polarization coordinate. Given this, 
it is a simple matter of high school algebra (equating reactant 
and product parabolic energy surfaces) to derive the important 
result for the activation energy, 


E, =(A+A4A (4) 


Thus the three quantities activation energy, exoergicity 
and reorganization energy are clearly related, and knowledge 
of any two will yield the third. Two important early contribu- 
tions led to increased understanding of the reorganization 
energies. Marcus used continuum dielectric theory to derive 
the result, for geometric models of precursor complexes for 
electron transfer, that 


o= (SZ) (Am)7F (RY (5) 


Here, €_ and €, are the optical and static dielectric constants 
of the solvent, An is the amount of charge transferred, and 
F(R) is the geometric factor depending on the relative sizes of 
reactants and products. This form gives the solvent contribu- 
tion (often referred to as the outer sphere contribution, hence 
i.) to the reorganization energy. Noel Hush pointed out’? that 
optical excitation from the reactant well yields an excitation 


whose frequency is /4 Ti for exchange reactions; this leads to 
a direct experimental estimate of the reorganization energy. 

These early understandings form the basis for all of mod- 
ern electron transfer theory. Nevertheless, there clearly were 
some difficulties and inadequacies of these early formulations; 
since the electron transfer reactions occur over such broad 
times and in such widely differing systems, it is not surprising 
that different control mechanisms are relevant. 


Tunneling Behavior: 
Non-adiabatic Transfer, 
Electron Tunneling, and 
Nuclear Tunneling 


The conceptual underpinning of equation 3 is that once 
the solvent has rearranged itself to attain the top of the activa- 
tion barrier in figure 1, the reaction will proceed. The number 
of times per second that the system attempts to cross the barrier 
is given by the prefactor Z, and the value 10'°/sec, generally 
assumed for Z, is a characteristic nuclear vibrational fre- 
quency. In a number of electron transfer systems, however, the 
system does not simply cross from the reactant side to the 
product side once the barrier top is attained. It was realized in 
the 1960s that electron transfers resemble other sorts of curve 
crossing phenomena - that is, there is an excited state curve 








Figure 2. 


Same as figure 1, but showing the excited state potential 
surface alsuv. The diabatic surfaces are defined by the light 
lines, the adiabatic ones by the heavy lines. The exoergicity 
AE” is roughly equal to the free-energy change AG? they differ 
by small entropic terms. The distance between the ground and 
excited state surfaces at the crossing point is 2Hrp. 











that “goes with” the ground state curve figure 1. Figure 2 
sketches the more general situation; note here that half the 
splitting between the two curves is given by H_, an electronic 
tunneling matrix element. Landau and Zener, in the 1930’s, 
had analyzed this sort of curve crossing process.** Under 
reasonable conditions, the rate of crossing from reactant to 
product on the lower potential surface will depend on H.» for 
small enough H_. For very large H_ (as is generally assumed 
in activated complex theory), the rates are independent of this 
electronic tunneling. On the other hand, for small H_ (such as 
are found in long distance electron transfers), the rate becomes 
proportional to the square of this electronic mixing matrix 
element. This sort of result is characteristic of quantum me- 
chanical perturbation theory, in which the rates are given by 
how strongly mixed the reactant and product are, times the 
density of product states, averaged over all of the reactant 
states. The generalization of the simple Marcus formula of 
equation 3 to the non-adiabatic regime is simply given by 


Mp | . 


vA ~ 
a Xk,T exp(-E,/RT ) © 


This rate law takes into account the electron tunneling, 
through the presence of the mixing element H.- Generally,’ 
one distinguishes electron transfer rates of eq. 6 as non-adiabatic 
transfer, as compared to adiabatic transfer, described by eq. 3. 


At very low temperatures, it was observed that many 
electron transfer rates are not activated — for example, figure 
3 shows the historical observation®> by Chance and DeVault 
of electron transfer process in chromatium; note that, though 
activated at high temperatures, the rate is essentially tempera- 
ture independent at low temperatures. To explain this, it is 
necessary to invoke quantum mechanical tunneling at low 
temperatures. Although classical chemical kinetics, as ex- 
pressed in equation 3, is based on the idea that this system will 
overcome the activation barrier in equation 1, quantum me- 
chanically the system can also tunnel through the barrier. 
While such tunneling contributions to chemical reactions have 
been known since Wigner’s work in the 1930’s, the under- 
standing of the role of nuclear tunneling in chemical reactions 
of electron transfer type is much newer. Vibronic models of 
electron transfer, that give nuclear tunneling contributions at 
low temperatures and go over to the appropriate Marcus-like 
form at high temperatures, were given in the 1970’s by a 
number of investigators, including Jortner,*® Fischer and Van 
Duyne,*”** Schmidt®? and (of course!) Marcus.'> These po- 
laron-type models***! are based on an assumed linear cou- 
pling between the electron population density and all the 
vibrations of the system, including solvent polarization (A) 
and intra-molecular vibrations (A,, or inner-shell reorganiza- 
tion energy). The formulas for non-adiabatic electron transfer 
rates that take nuclear tunneling into account properly, based 
on this polaron model, do not yield results that are simple as 
equations 6 or 3, though they go over to 6 in the appropriate 
high temperature limit. They do, however, correctly charac- 
terize the low temperature, non-activated rate seen in many 
non-adiabatic electron transfers.'® 


Solvent Dynamics, and 
Dynamic Gating 


Like activated complex theory itself, all of the preceding 
discussion has been based on relating the rate of a chemical 
reaction (in this case electron transfer) to the energetics in- 
volved in that chemical reaction. Both the original Marcus 
formulation and the polaron vibronic theory use parameters 
such as reorganization energy, exoergicity, and electronic tun- 
neling matrix elements that characterize the equilibrium prop- 
erties of the reactant or product species. A bit of reflection, 
however, suggests that if the rates at which the nuclear mo- 
tions, especially solvent rotations, explore their potential sur- 
faces are slow enough then the rates of such relaxation 
processes might dominate the overall electron transfer rate. In 
the Soviet Union, Burshtein and his students, especially Zus- 
man,” pointed out the role that solvent dynamics might play 
in reducing the rates of electron transfer: naively, one can 
simply imagine that, in the case of adiabatic electron transfer, 
when the system reaches the top of the barrier in figure 1, the 
solvent may not relax to the appropriate geometry to trap the 
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product on a rapid enough time scale. Under these conditions, 
the entire reaction dynamics can be dominated not by electron 
tunneling or nuclear attempt frequencies, but simply by how 
rapidly the solvent can relax into the geometry appropriate to 
the product state.'*'” An approximate equation for this limit 
is given by“ 


k 


k= 
ss 1+4x(H,,)? tM A 





(7) 


where k_, is the nonadiabatic rate of eq.(6), and T is a charac- 
teristic relaxation time. Thus when the relaxation time is very 
rapid, and the equilibrium state indeed always prevails, one 
retains the standard non-adiabatic rate expression as given, say, 
by vibronic theory; alternatively, when the solvent relaxation rate 
is slow compared to the non-adiabatic reaction rate, the rate 
constant is given simply by the rate of solvent relaxation.“° 
Experimentally, different fluids have different relaxation 
rates. Lovely picosecond and subpicosecond experiments by 
Fleming,” Barbara,’ Simon,** Weaver’ and others have 
demonstrated quite clearly the substantial regime in which 
solvent dynamics controls the rate of electron transfer. 


The Many Dimensions of Electron 
Transfer—Solvent Gating 


Figures 1 and 2 are highly schematic, in that the reaction 
coordinate is indicated as being unidimensional — in Marcus’ 
original work, it was a solvent polarization coordinate. The 
vibronic models, and the experimental realization that there 
were inner and outer sphere contributions to the reorganization 
energy, focused specific attention on the fact that many vibra- 
tional modes of the molecules, as well as polarization modes 
of the solvent, contributed to the reorganization energies and 
might contribute nuclear tunneling terms. Agmon and Hop- 
field suggested that differential friction, corresponding to 
diffusion rates along different nuclear distortion directions, 
might give highly complex chemical kinetics, such as is seen 
in, for example, ligand rebinding to myoglobin. In the simplest 
model for processes of this kind, one can envision two nuclear 
displacement coordinates coupled to the electron transfer re- 
action. For instance, figure 4 shows a schematic potential 
surface with the abscissa a fast, intra-molecular vibrational 
coordinate and the ordinate a slow, solvent conformational 
coordinate. If the rates of relaxation, or of diffusion, along 
these two coordinates are substantially different then the sys- 
tem point may not simply cross over the saddle point, or 
mini-max, on the potential surface - rather, it can cross the 
ridge line between reactant and product many times along the 
fast coordinate while slowly diffusing along the slow. Such 
behavior, whose theory has been extensively explored by a 
number of groups,*'**>? would among other things result in 
non-exponential rates of reaction. There are a number of other 
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Figure 3. 


Temperature dependence of the electron transfer in C. vino- 
sum cytochrome oxidase, along with theoretical fits. The 
dashed line is the experimental data from the original work of 
DeVault and Chance. Taken from reference 1, page 72. 
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situations in which reaction rates are also non-exponential, and 
some experimental evidence has indeed been found (again, 
often in biological systems) for such non-exponentiality. 

There is a second situation in which multi-dimensional 
behavior can yield surprises in the overall electron transfer 
rate. It was pointed out by groups at Northwestern™ and at 
Brookhaven” that subsidiary minima on the potential surface 
may permit the system to avoid what appears to be the saddle 
point, and, instead, to proceed from reactants to products by 
the motions that first move along one reaction coordinate and 
subsequently along another. This is indicated schematically in 
figure 5, and, indeed, such processes are well known in elec- 
trochemistry, where the kinetics that proceeds via the indicated 
subsidiary minima, rather than crossing over the global saddle 
point at the center of the diagram, is referred to as a square 
scheme”® process (perhaps because the system point motion 
looks like part of a square dance or a box step). 


Golden Rules and Correlation 
Functions - A General 
Dynamical Model 


The several stages of improved understanding that have 
been recounted so far lead to a contemporary view of electron 
transfer reactions that is both physically appropriate and com- 
putationally accurate. Still, however, several crucial aspects 
are missing. Perhaps the most important of these are that for 
dissociative electron transfer reactions, like that of equation 2, 
the potential surfaces are certainly not in any sense harmonic; 
the polaron model,***! which is the basis for the most sophis- 
ticated of the theories described thus far, assumes harmonic 





vibrations. Second, the electronic tunneling matrix element 
may well depend on the particular nuclear configuration - for 
instance, early work'*~’ on electron transfer between planar 
entities such as phthalocyanines indicates an important de- 
pendence of the electron transfer tunneling on the relative 
geometries, especially rotational geometries, of the two phtha- 
locyanine units. Third, general solvent behaviors should in- 
clude both homogeneous and inhomogeneous broadening, 
arising from solvent motions on the one hand, and from 
differing solvation energetics on the other.** All of these can 
be included,” within a general golden rule rate formula- 
tion, by replacing the energy conserving delta function that 
appears in Fermi’s golden rule for the rate constant by its time 
equivalent, and recasting the entire expression in terms of 
correlation function. This treatment was first done for spec- 
troscopy by Lax™ in the early 1950’s, and applied, very 
effectively, in a chemical context first by Gordon™ and later 
(using a particularly elegant wavepacket formulation) by Hel- 
ler.*© The generalization from spectroscopy to electron 
transfer has been recently given by cooperative efforts of 
groups from Evanston and Tel Aviv. *! 

The rate constant expression in the time-dependent correla- 
tion function context is quite simple, and is given by equation 8. 


k= if ae a ae i a 


There are several aspects of this expression that are quite 
remarkable. First, notice that the gap parameter A is now 
given as an explicitly time-dependent function; this permits 
information on correlations within the solvent, including 
inhomogeneous and homogeneous broadening, dephasing, 
and slow dynamics, to be included directly. Second, the 
tunneling matrix element appears in the correlation function 
— that is, when the averaging is performed it can be per- 
formed over whatever dynamical variables may affect the 
tunneling (that is, the Condon approximation is not in- 
voked). Third, the hamiltonians H, and H, of the initial and 
final states enter directly into the propagator — that is, one 
explicitly considers the time evolution of the initial and 
final states. If the potential surface is harmonic, this evolu- 
tion is essentially trivial; on the other hand, generalization 
to anharmonic behaviors, to multiple tunneling events, to 
superexchange, to dissociative electron transfer, and even 
to multiple minima is relatively straightforwardly handled 
within the time domain.*?*! 


Rate Control: Insights from 
Theory 
Nearly four decades of electron transfer theory have re- 


sulted in a clear picture of electron transfer rates at a number 
of different levels depending on the conditions of the actual 





Figure 4. 


A schematic sketch of the potential surface for ligand protein 
rebinding. R denotes the ridge, and S the saddlepoint. Stable 
reactant and product configurations are noted by a and b, 
respectively. Note that the shapes and orientations of the ridge 
are schematic, but also that the contour lines for diffusion are quite 
different in reactant and product states. From reference 51 














Figure 5. 


Schematic potential energy surface for electron transfer with 
conformational equilibrium. The abscissa and ordinate are, 
respectively, internal motion and configurational motion. The 
reaction goes from the reactant in the upper left to the product 
in the lower right; note the presence of the saddlepoint in the 
middle of the diagram, but also the two subsidiary minima in 
the upper right and lower left. Transfer from the precursor state 
in the upper left to the successor in the lower right can occur 
by a two step process, either across and down or down and 
across, rather than the direct path over the saddlepoint. From 
reference 54. 
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Figure 6. 


A series of synthetic intra-molecular electron transfer systems, 
containing rigid sigma bridges between donor and acceptor 
groups. From reference 68. 
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experiment. With respect to the traditional chemical view of 
rate determining steps, this theoretical picture permits us to 
characterize conditions under which particular physical ef- 
fects may control the rate of electron transfer. These control 
processes include: 


1) Franck-Condon factors, that limit the acceptance of the 
exoergicity by the vibrational modes of molecules and 
the solvents. This was the motivation for the gap law and 
the so-called Marcus inverted regime behavior, that 
reflect the counter-intuitive notion that under particular 
conditions increasing exoergicity results in decreasing 
rate. 


Solvent relaxation dynamics, that can slow down an 
electron transfer rate due to the inability of the solvent 
slow modes to follow the equilibrium evolution pre- 
dicted on the basis of the potential surface. These may 
be particularly important in polymeric, glassy, or protein 
systems. 


In non-adiabatic electron transfers, the electron transfer 
tunneling matrix element H» (whether direct or assisted 
by bridging groups) will limit how effectively the sys- 
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tem passes over the barrier top. Synthetic engineering 
based on the size of H,, permits fairly precise control of 
electron transfer rates. 


Gating effects are a particular manifestation of solvent 
dynamics, in which motion along one coordinate is so 
slow that it essentially determines the overall rate be- 
havior. Gating conditions are, again, expected in glass- 
forming materials. In particular, one strongly 
over-damped, slowly diffusing coordinate may com- 
pletely dominate the overall reaction, resulting in non- 
exponential behavior and very slow electron transfer. 
The presence of subsidiary minima on the potential 
surface can result in a completely different set of reac- 
tion pathways, leading to so-called “square scheme” 
dynamics. This can lead to the system completely avoid- 
ing the global saddle point, and, as a function of vari- 
ables such as solvation or energy, evolving either 
through the subsidiary minima or over the saddle point. 
Thus by emplacing subsidiary minima, the entire reac- 
tion mechanism can be changed. 


Energetic control. The overall exoergicity, A, is a critical 
factor in electron transfer reactions. Changing exoer- 
gicity can either increase (in the normal regime) or 
decrease (in the abnormal regime) the transfer rate. Both 
chemical synthesis and such weak effects as hydrogen 
bonding, selective solvation or ion pairing can substan- 
tially change the overall exoergicity, and therefore the 
rate of the reaction. One interesting form of energetic 
control involves shifts of local energy levels either by 
selective complexation™ or by ionic Coulomb interac- 
tion™ (“ion gated electron transfer”). 


Given these elements of rate control, it becomes a fasci- 
nating and promising synthetic target to design electron trans- 
fer reactions whose rates can be varied by the experimenter, 
either statically (particular design structures), or even dynami- 
cally (particular design structures whose behavior can be 
changed by controls such as photon switching, ion switching 
or voltage switching). 





Figure 7. 


A tetrad molecule, containing a carotenoid, a porphyrin, and 
a diquinone, synthesized to provide artificial photosynthesis - 
that is, long-life separation of an excited state into isolated 
cationic and anionic structures. From reference 30. 
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Controlled Electron Transfer 
Structures — Experimental 
Molecular Electronics 


By controlling the relative magnitudes of the six possibly 
rate controlling steps just discussed, it should be possible to 
prepare systems whose electron transfer rates can be control- 
led. The solid state analogs of such molecular systems have 
been known for a long time; indeed, transitor-based technol- 
ogy essentially arises from controlling directional electron 
flow across p/n junctions, analogous to electron donor/ac- 
ceptor interfaces. The first steps in attaining molecular devices 
based on controlled electron transfer entail demonstration of 
controlled electron transfer as an elementary process, in syn- 
thetic structures. In the remainder of this article, we will 
mention both completed experimental work on controlled 
electron transfer, and current investigations on application of 
controlled electron transfer rate processes to the understanding 
and preparation of molecule based devices. 


Controlled Single Step Electron Transfer 


Both intra and inter-molecular transfers have been studied 
in which rates can be changed, often by many orders of 
magnitude, by control of the exoergicity parameter . More 
recently, an elegant series of synthetic methods has culminated 
in the preparation of a number of species (figure 6) in which 
changing the nature of the bridge between donor and acceptor 
sites changes the efficiency of intra-molecular electron trans- 
fer.©” Finally, it has been shown that dropping the environ- 
mental temperature such that the solvent becomes glassy limits 
the available outer sphere reorganization energy, and can 
essentially turn off electron transfer processes.” 


Electronic Shuttles, and Molecular Wires 


Covalent structures to provide facilitated electron transfer 
between particular donor and sites has been used in 
the preparation of biosensor devices.>°”""-" Here, modified 
electrodes are attached to immobilized electron transfer en- 
zymes, by means of a molecular wire feature, that facilitates 
electron transfer between the sensing electrode and the redox 
enzyme. More generally, Lehn and co-workers” have inves- 
tigated conjugated systems for facilitating very long-range 
electron transfer, by what they call a “molecular wire.” 


Charge Separation and Storage 


The natural photosynthetic apparatus leads to separated 
electrical charges following the absorption of a photon. A 
number of artificial systems have now been prepared,>*”> in 
which charge separation is induced by controlled electron 
transfer and trapping, subsequent to initial excitation by pho- 
tons at a particular chosen wavelength. Obvious applications 
of such structures to energy transduction and solar energy 
Capture are under intensive investigation, using such species 
such as that shown in figure 7. 


Molecular Switches, Based on Two-Step 
Electron Transfer 


Wasielewski and collaborators”® have designed an elec- 
tron donor/acceptor/donor molecule, in figure 8, which exhib- 
its optical switching based on two ultrafast electron transfer 
reactions. The molecule absorbs light strongly at two different 
frequencies, 713 nm and 546 nm, in the singly and doubly 
reduced states of the electron acceptor end. This photophysical 
behavior is the basis for a light intensity dependent optical 
switch. As the light intensity increases, the molecule switches 
from being a strong transient absorber at 713 nm to a strong 





Figure 8. 


A bridged diporphyrinic system that constitutes a light-inten- 
sity dependent optical switch. This molecule switches its tran- 
sient absorption frequency with the intensity of the applied 
light field; after absorption that photon, an ion pair located on 
the left porphyrin and the bridge becomes stable, following 
which the second electron absorbs at the bluer wavelength. 
From reference 76. 
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Figure 9. 


The pi donor-sigma bridge-pi acceptor structure first pro- 
posed as a molecular rectifier. The acceptor on the left end 
has a lower unoccupied first level than the donor on the nght 
so that current would pass in a more facile fashion from left to 
right than vice versa. From reference 77. 
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Figure 10. 


The zwitterionic species recently demonstrated to provide 
molecular rectification. Except for the nature of the bridging 
structure, it bears a remarkable similarity to the species of 
figure 8. From reference 78. 
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absorber at 546 nm. The molecule can then be used to modu- 
late two beams of light at different colors on the picosecond 
time scale and even to perform logic operations. This stands 
as an interesting prototype of integrated molecular opto-elec- 
tronic devices, that would be both very small in size and very 
fast in performance. 


Molecular Rectification 


It was proposed twenty years ago” that use of an intra- 
molecular donor/acceptor structure with a bridge that permits 
only weak mixing of donor/acceptor levels could result in 
rectification of applied current. The first proposed structure is 
shown in figure 9: the forward bias would result in oxidation 
of the donor and reduction of the acceptor, but on reverse bias 
the more energetically demanding processes of reduction of 
the donor and oxidation of the acceptor would not be observed 
until a much higher driving field is applied. Very recently, the 
structure shown in figure 10 has been prepared,”* and appar- 
ently demonstrated to act, indeed, as a molecular rectifier 
(figure 11). If verified, this would be an important advance in 
the understanding and preparation of single molecule electron- 
ics. 


Macroscopic and Microscopic 
Coupling — Scanning Tunneling 
Microscopy and Electron Transfer 


Single molecule electronics has as one goal the prepara- 
tion of molecular circuits that operate in a fashion similar to 
those of macroscopic solid state circuits. A close analogy 
between electron flow in engineered molecular structures and 
in solid state device applications is strikingly indicated by 
consideration of scanning tunneling microscopy.’**° In this 
technique, electrical currents pass between an electrode and a 
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metallic or semi-conductive tip; the intervening structures 
facilitate the electron transfer process, and study of the cur- 
rents as a function of applied voltage, or of tip to electrode 
distance, provides information (in the form of images) on the 
dimension and nature of the intervening material. In the very 
weak perturbation limit,®’ the STM current is proportional to 
the density of electronic states on the electrode/overlayer 
structure, evaluated at the position of the tip and the energy of 
the applied voltage. If, however, the tip is brought closer to the 
surface, simple perturbation analysis fails and one must con- 
sider the direct calculation of the tunneling current.*? 

A simple and effective picture for discussing the STM 
experiment, in the limit of small applied voltage, involves 
transport through an overlayer consisting of discrete local 
electronic levels (such as atomic orbitals) coupled by usual, 
extended Huckel-type tunneling matrix elements. At the elec- 
trode end and the tip, these local sites are coupled to continua, 
corresponding to the conduction bands biased by the applied 
potential. Under these conditions, the rate can be described** 
in terms of an effective scattering matrix element that acts 
between an initial state on the electrode and a final state that 
has traversed from electrode through overlayer to tip. Devel- 
opment of such a formalism, based™ on applications of formal 
scattering theory, is in fact fairly straightforward.*° 





Figure 11. 


he measured current voltage curve from a structure consisting 
of a silver surface, several monolayers of the molecule in figure 
10, and a magnesium junction. The voltage sweep was set at 
44mV/s. Note the obvious rectification behavior, with forward 
bias providing much larger current. From reference 78 
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The striking resemblance between STM and non-adi- 
abatic electron transfer, at a molecular level, is that in the 
electron transfer process one again deals with a series of 
tunneling matrix elements between local atomic orbital sites. 
In place of the electronic continuum, however, these electronic 
structures are coupled to a quasi-continuum of vibrational 
energy states, and the Franck-Condon weighted density of 
States structures replace the continua arising, in the STM 
experiment, from the metal bands. This analogy suggests not 
only that similar formalisms can be used, but that the distin- 
guishing hallmarks of the electron transfer process (inverted 
region slow down, and the other rate-controlling steps men- 
tioned above), can be observed, appropriately, in STM meas- 
urements. Conversely, it also appears that some of the 
characteristic behavior of mesoscopic structures, in particular 
formally negative resistances, coulomb blockades, and reso- 
nance assisted tunneling processes should be observed in 
electron transfer phenomena. These lead to an additional level 
of control, and may well provide a very rich area for real 
molecular electronics. 

Electron transfer or, more generally, the flow of electrical 
charge in designed structures, lies at the very heart of current 
research and applications in chemistry, materials science, bio- 
chemistry and solid state physics. Exploration of the common- 
alities of these processes, and possible control mechanisms 
and design structures, makes electron transfer research one of 
the most exciting, forefront topics in each of these fields. 
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Solvent Effects in Quantum 
Chemistry: The Self- 
Consistent Reaction Field 


Model 


Michael C. Zerner 
University of Florida 


Nearly all of quantum chemistry relates directly to the 
chemistry of isolated molecules in the gas phase at absolute 
zero temperature. In spite of the fact that quantum chemistry 
and quantum chemical concepts now permeate the chemical 
literature, most everyday chemistry takes place in solution, 
where it is influenced, sometimes considerably, by the sur- 
rounding solvent. This article presents some simple concepts 
that help to bridge the gap between quantum chemistry and the 
chemistry that most of our colleagues examine experimentally. 
We examine here the reaction field model, and show the impact 
these considerations have on structure, thermochemistry, and 
on electronic spectroscopy. Examples range from simple 
equilibria to the photochemistry involved in the initial photo- 
chemical event in photosynthesis. 


Introduction. 


Solvent effects. 


Now this seems at first glance a subject far from the heart of 
most quantum chemists, and I am a Quantum Chemist. By-and- 
large quantum chemistry refers to isolated molecules in the gas 
phase. And the history of modern physical chemistry focuses 
around isolating systems — low pressures, molecular beams, 
matrix isolation. Everything possible is being done by modern 
physical chemists to uncover the “basic” nature of a molecule or 
a pair of molecules reacting in isolation — free from complicating 
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environmental factors that make these studies complex. It is 
these “modern” physical chemists that most of my crowd want 
to talk to. These 200 or so scientists are creating the data where 
we can test our “fundamental” theories. Unfortunately, most 
of chemistry simply isn’t like this. It is messy. Molecules are 
generally surrounded by many others, in “solution”, and they 
react with others also in solution. For every low pressure gas 
phase physical chemist, there are hundreds of bench chemists 
happily creating things, and their results are not directly pre- 
dicted by modern quantum chemistry if solvent effects are 
ignored. Quite the contrary. There are now many examples of 
state-of-the-art quantum chemistry yielding incorrect results. 
Not necessarily incorrect results for the gas phase experiment, 
which usually has not been done, but incorrect when compared 
to the only experiments that have been done — in solution, or 
surrounded by protein material, or on surfaces. These environ- 
mental effects are not necessarily small. 

That these effects are not small is hardly news to biochem- 
ists Or to anyone interested in charge or electron transfer. 
Fundamental processes that occur in biology or in solution 
often do not occur in the gas phase. Certainly solvent partici- 
pation, or more precisely dielectric relaxation of a solution or 
surrounding protein materials is at the heart of the theories 
developed by Rudy Marcus for charge transfer’, our unusual 
colleague to which this article is devoted. 

That solvent effects are not necessarily small may not be 
news to the biochemist, or the statistical mechanician or Rudy 





Marcus. It should not even have been news to me, as I follow 
much of the statistical mechanical literature. But it was. It was 
a surprise to me when some of the things I was doing didn’t 
turn out as expected. I will write briefly about two of these 
below: merostabilization of dyes, and the initial photochemi- 
cal event in photosynthesis. Both depend crucially on the 
surrounding media. 

Now what are these solvent effects? Below is a useful 
classification, although it is certainly ordered in a way useful 
for a quantum chemist. 


None. 


Then, of course, this paper need not be written! There is 
a chemistry at low pressure with but small environmental 
effects, the experiments we alluded to above. Also solvent 
effects on non-polar systems MIGHT have only a slight effect 
on structure or reactivity, and even vibrational spectroscopy. 
But as more careful experiments are done, we notice that even 
systems in “matrix isolation” are not as isolated as one might 
hope. 


Weak. 


By weak we usually refer to binding energies between solute 
and solvent below thermal energies, below “kT”, the Boltsman 


constant time temperature. Since liquids generally have binding 
energies also of the order of kT (or they would be solids or gases 
at room temperature), they are influenced by the field of the 
solute, but continue to maintain to good degree of approximation 


the properties of the pure solvent. From the quantum mechani- 
cal point of view we want to be able to consider the solvent 
and solute as separable. This implies no charge transfer and no 
exchange between solvent and solute. 


Strong. 


Here strong refers to interactions greater than KT. In such 
acase the nearest neighbor solvent molecules have a relatively 
fixed geometric relationship with the solute. The solute is best 
considered as the dissolved molecule and its strongly bound 
neighbors. The binding between this super-molecule and the 
rest of the solvent is then weak. 


Reaction. 


The reacting species must all be considered as a super- 
molecule, and the rest weakly interacting. The reacting species 
may include the solvent directly, or as a catalyst. For example, 
many proton transfer reaction exchange protons with water 
when it is present *. 


Mechanical. 


Two species can be trapped in a solvent cage and forced 
to influence one-another, or even react with one another, as 
they are in proximity for a longer time than they would be in 
the gas phase. For example, they collide with one another, and 
would bounce apart if it were not for the fact that they collide 
with the solvent and bounce together again. 





Figure 1. 


Three regions to be considered in simulating solvent effects. 
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Figure 2. 


A hypothetical charge distribution of a molecule in a cavity, and 
the image charges generated from in a dielectric continuum. 
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There is an important operational difference between 
weakly and strongly interacting species. In the former the 
changing orientation between solute and solvent, and solvent 
and solvent, is important, and effects the properties of the 
solute. This is correctly modeled only through a careful statis- 
tical mechanical averaging of all possible orientations. In the 
latter, a fixed geometric relationship generally pertains, and 
this is the usual domain of quantum chemistry. 

Much can happen to a molecule as it dissolves in solution. 
It can dissociate, react, complex, form associations or dimers, 
etc. Still all species present can be treated as weakly or strongly 
binding with the solvent, if the reaction is slow, or if it is fast. 
This becomes more difficult if what happens happens on the 
same time scale as do the thermal motions of the solvent itself. 


How does one treat a “solution?” 


In principle, we should perform: 

1. A Monte-Carlo or molecular dynamics study of the 
solute and solvent using potentials derived from quantum 
mechanics. In general, we cannot perform a study of this 
nature on Avogadro’s number of molecules to mimic systems 
of macroscopic size (we would need 6x10” molecules, and 
just keeping track of so many coordinates exceeds the memory 
of any computer). Generally a solvent and, say, 200 solute 
molecules with the appropriate boundary conditions that 
mimic the remaining solution is sufficient to accurately repre- 
sent the thermodynamics *:*. But now the quantum chemistry 
becomes a problem. The thermodynamics requires that one 
study millions of spatial configurations in order to calculate 
the average macroscopic property. Since there are about 31 
million seconds in the year, if we could perform each calcula- 





Figure 3. 


The Born model of solvation for an ion, see text. 
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tion in 30 seconds this becomes feasible if we had a computer 
dedicated to this task. Hartree-Fock Molecular-orbital inde- 
pendent-electron calculations using semi-empirical quantum 
chemistry can perform in this time frame on these big systems 
on today’s largest computers. But the interactions in the system 
are weak, and it is dubious if independent-particle calculations 
are of any value as they do not correctly yield weak interac- 
tions. What will be needed are “correlated” theories, those that 
take more specific care of electrons that then yield correct 
weak interactions. Considering the problems that then arise 
with basis sets, basis set errors and corrections, correlation —* 
® this is not possible today and not in the near future. 

2. Monte-Carlo or molecular dynamic studies with clas- 
sical potentials. This is a well established branch of theoretical 
chemistry, and has been used to yield a great wealth of infor- 
mation on solvation and on biological interesting systems ***. 
The accuracy of these results rely on the accuracy of the 
potentials assumed, and these potentials are generally selected 
to reproduce a set of well known properties. These potentials 
are relatively simple, treating atoms as balls, and bonds as 
springs. This is a technique generally for weak interactions. 
Classical mechanics can not be used to correctly describe 
reactions or electronic spectroscopy, for examples, as these 
phenomena are electronic in nature. Atoms are not balls, and 
bonds are not springs: molecules are made of electrons and 
nuclei. 

There is a middle ground here, and one that is beginning 
to be successfully exploited. Here one obtains the potential for 
isolated pairs (or even triplets) of molecules, and than fits these 
potentials, without necessarily recourse to ball and spring 
models. This was shown very successfully in modeling the 
famous SN2 reaction * 


Cr+H,C-Cl >CI-CH,+Cr 


which reacts 20 orders of magnitude more rapidly in water 
than in the gas phase. The water-water potential was treated 
from a fit classical potential, while the potential for the reac- 
tion itself came from quantum mechanics. 

3. Cluster models. Here one surrounds the solute molecule 
with an ever increasing number of solute molecules and treats 
the system using quantum mechanical potentials and statistical 
averaging. One then watches for the convergence of calculated 
properties with cluster size. 

This is a very difficult approach, and is fraught with the 
very many pitfalls that have accompanied similar studies on 
solids: how many atoms or molecules make a solid? This is 
still an unanswered (unanswerable?) question, and the far 
greater thermal motion of a solution will make the answer to 
this question even more difficult. Most likely we will not like 
the answer. 

4. Peturbative models. A picture is given in Figure 1. A 
quantum mechanical supermolecule is surrounded by a region 





of solvent molecules treated by classical potentials, and the 
entire system placed in a dielectric continuum. 

Examining the utility of the dielectric continuum model 
will form the bulk of this manuscript. The energy (Hamil- 
tonian) of this system can be written as the sum of its parts plus 
the pairwise interactions between parts. 


H=H(QM)+H(Classical)+H(QM-Classical)}+H(QM-Dielectric) 
+H(Classical-Dielectric) (1) 


This is a particularly easy way to describe the problem of 
solvation and to decompose the problem into more manage- 
able pieces. In some sense, it has no content, for it clearly 
includes case 1., above, if the QM-supermolecule is made big 
enough, case 2. if the QM moiety is dropped completely, case 
3 if the QM-supermolecule is developed as an isolated cluster. 

In this work we will focus mostly on the quantum me- 
chanical supermolecule within the dielectric continuum. This 
is the easiest to describe, and already yields interesting results. 
In the final section we will describe calculations for spectros- 
copy that contain all three moieties. 


The Dielectric Model for 
Ground State Properties 


When a molecule is dissolved in a dielectric material, the 
electric field it establishes sets up a reaction field, one that 
counters it, and stabilizes it. Figure 2 pictorially represents 
this. Figure 3 presents the Born model ° for the solvation of an 
ion. The ion and atom are assumed to be spherical. The net 
energy for the Born-Habor cycle is then given by 


AH(sol)=AH(gas)+AH(cavitation)+AH(solvent) (2) 


~ -Q’/2a+0.0+Q’/2ea 


where a is the radius of the ionic sphere, Q is the charge of the 
ion, and the cavitation energy, a few kcal/mol., is forgotten for 
the moment when examined in the context of these other larger 
numbers. AH(gas) and AH(sol) are the energies required to 
charge a sphere of radius a reversibly, i.e. 


Q 
aH=-J adq (3) 


where € is the dielectric constant of the medium, unity for the 
gas phase, 80 for water at room temperature, 2.2 for a non-po- 
lar solvent like cyclohexane. 

Assuming a simple ion, the radius of the largest orbital 
can be estimated from 


<r> = (2n+1)/2 (4) 


where n is the principle quantum number (n=1 for Li*) and 
is the orbital exponent for this orbital, which can be, for 
illustrative purposes, taken from the Zener-Slater screening 
rules (€ =Z-0.3 for 1s orbital, 2.7 for Li*)'® "". This yields for 
Li* a solvation energy of 13.1 eV(leV=23kcal/mol.). An ex- 
perimental value of 15.1 is reported for Li. This is remarkably 
close for such a simple theory. On the other hand, 
2eV.=46kcal/mol., a large error for some purposes! Such a 
simple and appealing theory has, as might be expected, led to 
tables of effective radii for atoms, ions and molecular groups 
fit to reproduce experimental heats of solution’. 

This theory can, and has been generalized by Kirkwood". 
The interaction between a molecule and its environment is 
given by". 


E=-1/2 Jpmupe) dV (5) 


where p(r) is the electronic density of the molecule (or super- 
molecule) and U, is the reaction potential that this charge 
distribution generates within the medium. Assuming the solute 
is embedded in a sphere of radius a, and expanding the ficld 
generated by the molecular density in terms of electric mo- 
ments yields the expression’? 


E=-122e TMM, 


l m=-] 





Figure 4. 


Amore realistic model for the shape of a molecule in a solvent 
is generated from surrounding each atom by a sphere, rather 
than the entire molecule by a sphere or an ellipse. The external 
surface of each sphere is then divided into patches in which 
the surface charge density is estimated, and the boundary 
equation 9 is solved. 
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In the above, M, _, are the moments and g, are the reaction 
“tensor” given by 


8, = 2 [0+1)(e-1) } /[ (+ 1e+! J (6b) 


The first term is simply the Born charge term. The second 
term, the dipolar term, was first explored by Onsager". 


E = 1/2Q’/a - 1/2 g,(e) pep /a® (7) 


where 1 is the molecular dipole with components M, |= px, 
M, = Hy and M, |= iz, etc., and, from the above equation, 


g(e) = 2(e-1) / (2e + 1) (8) 


Analytical expressions analogous to those above can also 
be derived for an ellipse’* '* '”, More general shapes can be 
generated also, as depicted in Figure 4, in which each molecu- 
lar group is surrounded by a sphere of a given size (either a 
Van Der Waal radius or a radius fit to reproduce experimental 
values'* '*). In this case, the electrostatic potential must be 
generated on each surface point, and the reaction charge 
“facing” it generated by the continuity equation. 


e(internal)(5V/Sv)internal = e(external) (85V/5n)external(9) 


In equation 9, the derivatives represent the electrostatic 
potentials normal to the surface, internal to the sphere, with 
usually equal to unity, and external to the sphere, with e(exter- 
nal) usually equal to the bulk dielectric constant. 

Generally this is done using finite element procedures'*. 
The generation of the electrostatic potential at the surface is 
difficult, and has been approximated often just from atomic 
charges’®, and the number of surfaces patches often suffi- 





Figure 5. 


A pictorial representation of mero-stability. Molecules 1 and 2, 
symmetrically substituted, have small dipole moments along 
the symmetry axis, here “z". The asymmetric molecule 3 has 
a much larger dipole moment. Since the stabilization in a 
dielectric medium is roughly proportional to the dipole moment 
squared, molecule 3 is stabilized more than 1 or 2. 





My and wu? < 3 
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Table I: 


A Study of Merostability in the Gas Phase (e=1) and in “Water 








e=1 =80 


Molecule kcal/mol kcal/mol 





CH(CN)OH 
INDO 3.1 
4-31G 3.3 
CH(CN)OCH, 
INDO 1.2 
CH(CN)NH, 
INDO 5.7 13 
4-31G 58 -7.0 


-24.4 
-22.2 


-28.8 





ciently large as to suggest further tricks to cause more rapid 
convergence. 

Another criticism of this approach, in addition to attempts 
to describe the shape of wiggling solute molecules in a solvent, 
is that the dielectric strength passes directly from that of the 
gas phase to that of bulk dielectric. One correction is to create 
two spheres (or ellipsoids) and more gradually pass from 
dielectric constant 1 to €’ to € (bulk)”®. This is somewhat more 
successful, but involves two additional parameters, €’ and the 
corresponding a’. Amore realistic approach would be to model 
the dielectric as a value that smoothly went from 1 to bulk as 
a function or r. This has also been developed successfully”! but 
has not yet been systematically examined. 

One can establish in the context of quantum mechanics 
the equivalent Hamiltonian (Schrodinger equation), and solve 
the equation 


Hy = Ew (10a) 


H=Ho- 1/2 g(e)/ <p>. wa® (10b) 


in which Ho is the gas phase super-molecular Hamiltonian 
operator, and the remaining terms mimic the effects of solva- 
tion. y is the many electron wave function. The charge term 
is classical, and can be added to the energy at the end of the 
calculation (i.e., the net charge on the supermolecule does not 
depend on the quantum mechanics). If the arbitrary shape of 
Figure 4 is used, then 


H=Ho- 1/2 >, Dalia) /r(ia) 


a ie€a 


(10c) 


is used, where q(i,a) are the surface charges established 
through solving the continuum equation 9, and the sum if over 
all point “i” generated on the patches “a” on the surface. 
Equations 10 are iterative in nature, as we do not know 
the charge distribution of the molecule that generates the 





moment expansion of equation 10b when we begin, nor do we 
know the electrostatic potential needed for equation 10c. But 
this is not a real concem, as the solution of Ho, itself, is also 
iterative in nature. This includes (along with the sum of the 
momentum of the electrons, the electron nuclear electrostatic 
attraction, and the additive nuclear-nuclear repulsion) elec- 
tron-electron repulsion, and this also can not be established 
without the electron distribution. 

This is a brief introduction to the major ideas behind the 
quantum mechanical description of the “self consistent reac- 
tion field model, SCRF. 


Some Examples of the Simple 
SCRF Theory for Ground 
states: 


Table Il: 





Reaction Field Effects on LogK from AM1 Calculations. Experi- 
mental values are in parentheses*. 





e=1 e€=78 
(gas) (“water”) 
0.4 -2.8 

(0.4) (-3.1) 


Equilibrium 





1.7 47 
(4.7) 


-3.0 
(-3.0) 


6.5 -73 
(-4.6) 
a. From reference 26. In this reference 12 systems are ex- 
amined, of which the above are typical. The experimental 
values are from A.R. Katritzky, Handbook of Heterocyclic 
Chemistry, Pergamon Press, New York (1985) pp. 47-50. 








Table til. 


An Estimate of Boiling Temperatures from Cavitation Energy 
(Scaled Particle Theory) and Calculated Dispersion Ener- 
gies* 





calc. (K) obs. (K) 
CH, 123 113 
CH, 177 188 
Benzene 327 353 
Naphthalene 520 491 
*) From M. Karelson and M. C. Zerner, work in progress. 





There are many examples obtained from SCRF theory”, 
but I first got “hooked” while trying to examine a phenomena 
called mero-stability”. This has to do with the stability of a 
series of dyes with general formula X(CH)Y. Merostability is 
then defined as 


. 2 om 
E(mero) = E(XCY) - [E(XCX) + E(YCY)}/2 (11) 
H H H 
This is a quantum mechanician dream system! All system- 
atic errors should cancel (at least hopefully) in the subtraction. 
This is borne out. Ab-initio calculations of nearly any quality 
yield nearly the same results as do semi-empirical calcula- 
tions-and ALL of these calculations showed that E(mero) is 
clearly positive. There is less of the mixed species present than 
the symmetrically substituted ones. By now my poor colleague 
Alan Katrisky™, was having his experiments attacked “left 
and right” by knowledgeable theorists. But now we are pre- 
pared to worry about the solvent. Table 1 shows a few exam- 
ples from our study”. Figure 5 shows our explanation for these 
results. The asymmetric substitution in these systems gener- 
ates a large charge asymmetry, and this receives considerable 
stabilization in water. The symmetrically substituted com- 
pounds have only relatively small dipole moments along their 





Figure 6. 


A comparison of the free energy of solvation with AEag calcu- 
lated from AHe-80) - AHe=1) using the AM1 model and the 
mult-cavity self consistent reaction field model; from reference 
26. Molecules range from CH4, point 2 and CH3CONHe, point 11 





AE. 














AG(g--w) 
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Figure 7. 


The influence of a dielectric medium on absorption spectros- 
copy. The ground state of the chromophore is assumed to be in 
equilibrium with the solvent. The absorption process is assumed 
so fast that only the electronic polarization of the solvent can 
adjust. The image dipole only partally rotates in this short time. 
Complete rotation requires the nuclei of the solvent to translate 
and rotate, a process much slower than the absorption process. 


Y 





YY 


o 


A simulated spectrum of pyrimidine in water. The spectrum is 
calculated using the INDO/S model (J. E. Ridley and M. C. 
Zerner, Theoretica Chemica Acta (1973) 32, 111) and includes 
pyrimidine and two water molecules all embedded in a reac- 
tion field simulating bulk water. The geometric conformations 
were generated using the water-shell option of QUIPU (from 
George Purvis Ill). 


Figure 8. 
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35,300 35,400 


frequency in crv 1 


35,500 


axis of symmetry and receive a relatively small amount of 
stabilization. 

Table 2 shows a study of the relative stabilities of the 
tautometers of hydroxy-pyridine-pyridone systems”®. These 
systems are of importance in biology and many studies have 
been made of them, some with the best of quantum mechanical 
technology. These studies yielded incorrect results. Again, as 
shown in this table, solvation favors the largest electrical 
asymmetry and yields quantitative results for the Log(K), with 
K the equilibrium constant, the ratio of product to reactant. A 
negative Log(K) means that the molecule on the right (reac- 
tant) side of the equilibrium dominates, a positive value that 
the molecule on the left(product) side dominates. 

A final example for ground state properties examines the 
use of a multi-cavity SCRF, MCa-SCRF, scheme for solva- 
tion””. In this, 


AG(sol)=AG(es)+AG(pol)+AG(cav)}+AG(disp}+AG(S) (12) 


AG(es) from the pure electrostatics will come from our self- 
consistent reaction field calculations. AG(pol) refers to polari- 
zation effects-changes in the electronic distribution due to the 
presence of the solvent, and this is also included (this is the 





Table IV: 


Some Calculated Spectroscopic Shifts in x - x* Transitions, in 
1000cm.'. The quantum chemistry is from ZINDO, reference 
42. See also Reference 32. 








Molecule 


OO) 


CH 


UNM 


Calc 
1. Gas 20.1 
2. CHCI, 21.8 
3. water 24.6 
«(3)-@(2) 2.8 





1. Gas 32.0 
2.C,Hi, 32.8 
3. water 29.1 


w(3)-@(2) -3.7 


1. Gas 
LCS 


29.7 

ra 263 274° 
3. water 223 233 
w3)-(2) 4.04.1 


H5C2 
\ 
H5C2 


a. Daehne, S.; Schob, F.; Nolte, K.-D. Z. Chem. 1973, 
13,471 

b. Kosower, E.M.,; Ramsey, B.G. J. Am. Chem. Soc. 1959, 
81,856 

c. Brooker, L.G.S.; Craig, A.C.; Haseltine, D.W.; Jemkins, 
P.W.; Lincoln, L.L. J. Am. Chem. Soc. 1965, 87,2443 
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Table V: 


INDO/S SCRF CI Calculated and Experimental n-x Transition Energies of Some Heterocycles in Different Media (from Reference 34). 
Molecule Solvent € n,, vem") 
1. Pyrimidine Gas Phase 1.000 1.0000 32966 

Isooctane 1.940 1.3915 33559 

Diethyl! Ether 4.355 1.3527 34127 

Acetonitrile 37.50 1.3416 34697 

Water 80.10 1.3330 34743 

2H,0 1.000 1.0000 30982 

Water + 2H,0' 80.10 1.3330 36572 








Gas Phase 1.000 1.0000 28329 
Isooctane 1.940 1.3915 29460 
Diethyl Ether 4.355 1.3527 30382 
Acetonitrile 37.50 1.3416 31296 
Water 80.10 1.3330 31368 
2H,O0 1.000 1.0000 26490 
Water + 2H,O 80.10 1.3330 33927 


Gas Phase 1.000 1.0000 30387 
Isooctane 1.940 1.3915 30387 
Diethyl Ether 4.355 1.3527 30387 
Acetonitrile 37.50 1.3416 30387 
Water 80.10 1.3330 30387 
2H,,0* 1.000 1.0000 32900 
Water + 2H,0° 80.10 1.3330 33301 


'The calculation gives two close transitions at 35465 and 37678 cm" of equal oscillator strength. In the table the mean value 
of transition energy is given. 


*The calculation gives two close transitions at 30505 and 33562 cm" of inequal oscillator strength. In the table the weighted 
mean value of transition energy is given. 








Table Vi: 


Calculated Absorption Spectra of the Reaction Center of R. Viridis, See Figure 9 





Isolated - Molecule, Solvent - Simulated 
Avicm') _Av(Debye) Avicm') _Av(Debye) 
11,600 46 11,400 46 

19,700 789 12,600 75.0 
19,800 49.9 16,600 46.1 
20,600 769 13,500 748 
20,200 55.5 16,700 484 
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“self consistent” part. A non-iterative solution of equations 10 
would not include this term). When a cavity is made in the 
solvent, this requires energy, represented by AG(cav)”. Fi- 
nally, AG(S) is a term representing changes in the solvent 
structure due to the presence of the solute. 

Experimental AG(sol) are plotted against the results from 
the MCa-SCRF approach calculated as AE(sol)= AH,, (€=80)- 
AH,, (e=1) in Figure 6. The cavities surround each rotationally 
independent moiety of the molecule: for example, CH,F is one 
cavity, CH,-CH, two cavities and the CH,-CO-CH, three 
cavities. The total volume of each molecule is determined from 
mass density and the volume of a CH, group determined as 
half that of CH,-CH,, etc. Such a plot assumes that the cavi- 
tation energy and the dispersion contributions are roughly the 
same and of opposite sign, which they are, and cancel, and that 
the entropy term is similar for all these similar molecules. The 
linearity of this NO FREE PARAMETER plot is remarkable. 

The single cavity SCRF model gives about the same result 
for all of these systems (at least up to the dipole term), 
demonstrating that it is important to identify rotationally free 
sections of a molecular system that have orientation times of 
the order of the orientation times of the solvent itself. They 
establish independent “reaction fields”. 

Recently we have established a theory to calculate the 
dispersion energies for a molecule in a reaction field”® *°. 
Using these values, and cavitation energies from scaled parti- 


cle theory”*, we can estimate heats of vaporization for a pure 
solvent. If there is little polarity, then equation 12 becomes 


AG(sol) = AG(cav) = AG(disp) (13) 


as the AG(S) term is zero. We are examining this equation now, 
but an interesting sideline is, that at the boiling temperature, 


AG(cav) = -AG(disp) 


again for a non-polar solvent. Since AG(cav) at most tempera- 
tures is a linear function of T from scaled particle theory, we 


can estimate the boiling points of some pure solvents. Boiling 
point must be considered one of the most colligative of all 
properties, and this can be estimated from a single QM calcu- 
lation! Preliminary results appear in Table 3. Mati Karelson 
and I stress preliminary, for we have not had time for a careful 
examination of these results, nor what is behind the choice of 
a number density that is required in scaled particle theory *'. 


The Dielectric Model for 
Excited State Properties: 


It is well known that the spectrum (color) of a molecule 
is different in the gas phase than in solution. Solution spectros- 
copy is everyday business in any experimental lab: vapor 
phase spectroscopy is already a specialized topic. Can we 
predict the spectrum of a molecule, or, at least, what size 
“solvent shift” to expect? 

A molecule can absorb light, and then re-admit light. 
Absorption spectra is generally estimated using the Franck- 
Condon principle that states that the nuclei do not move during 
the absorption process. Absorption of light is general of the 
order of 10'* seconds, much more rapid than the molecular 
vibrations required to move nuclei relative to one-another, 
generally of the order of 10°? seconds, or 100,000 times 
slower. Once excited, however, the excited state can often live 
from seconds to, perhaps, 10°* seconds, and this is enough time 
for the molecule to assume the geometry of the excited state. 
For this reason, emission spectroscopy is calculated from the 
geometry of the excited state. This idea of relative time scales 
becomes important when considering the effects on solvation 
on absorption or emission. The model we will assume for 
absorption is that the solute molecule (the chromophore) is in 
equilibrium with the solvent before it absorbs radiation. The 
absorption process is then assumed so rapid that the solvent 
molecules cannot rotate or translate during the absorption 
process. The orientation and translational degrees of freedom 
of the solvent that make up the dielectric response of the 





Table VII: 


Simulated Spectra of CH2O in Water, Shift in cri'. 





Supermolecule, 
CH,0...H,0*) 
Supermolecule in SCRF 
CH,0O...H,O*) 
Simulation 
Scaled Changes”) 
Unscaled Changes 


385 


1130 


1200 
1900 


*) An average of the ten lowest lying conformations, see text. Only one H,0 is kept from the simulation 
» The scaling is simply to set e(R)=R. The charges are chosen to reproduce the dipole moment of water in water. 
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solvent are frozen. However, the electronic response of the 
solvent is considered instantaneous, and we assume that the 
electrons do follow the chromophore excitation: i.e., part of 
the solvent responds to the excitation during the absorption 
process-the electronic polarization, another part is too slow- 
the orientation and translational polarization. 

The dielectric constant consists of components of differ- 
ent frequencies. Here we need consider only two-that slow 
relative to the excitation lifetime and that rapid relative to it. 
This can be expressed as '*:*2 


g(e) = g(D) + g(n(D)’) (14) 


where D is the component of the dielectric response that is 
slow, rotations, vibrations and translations, and g(n(D)”) is the 
part that is rapid. Here n(D) is the index of refraction, and 
n(D)’ approximates the instantaneous dielectric response of 
the solvent'*: *2, 

Proceeding with this model, we solve the many electron 
problem as before in the presence of the SCRF. We then must 
correct the energy of the excited state (to first order) by 
subtracting the electronic polarization contribution that is 
proportional to the charge distribution of the ground state, and 
that is automatically incorporated in solving a Schrodinger 
equation having a Hamiltonian of the form of Equation 10, and 
then adding back that part of the solvent electronic polariza- 
tion that senses the excited state charge distribution. 


AE(absorption)=AE(SCRF)-1/2g1(n(D)*)<*>e(<p*>-<y10>){ 15a) 


where <j1*>is the dipole moment of the excited state, etc. This 
of the ground process is shown schematically in Figure 7. 

This is not the only way in which we can view the 
absorption process. For example, we could also view this 
process as if the electrons of the solute are always in equilib- 
rium with those of the solvent. This leads to a somewhat 
different results than that above”? 


AE(absorption)= AE(SCRF) - 1/4gi(n(D)’(<p*>-<po>)> (15) 


and suggests that the electron of the solvent sense both the 
ground state and excited state charge distribution in an average 
way. As the terms that we add in either 15a or 15b are similar 
in most cases, and small relative to the shifts already included 
in AE(SCRF) it is difficult to decide between these two mod- 


Using the model suggested by equation 15a., we reported 
the spectra of some pi->pi* excitations that appear in Table 4 
that show unusually large solvent shifts’. The results are 
remarkably good for such a simple straight-forward model, and 
we must assume that many of the errors inherent in such an 
approach must cancel when we subtract excited state energies 
from those of the ground state. This would include errors in 
molecule shape (these molecules are certainly not spherical), 
errors in the simple time development assumed, and errors from 


omission of the cavitation and entropic terms of equation 12. 
The dispersion energy is not constant, and so the shifts we 
compare are between two different solvents, and not between 
solvent and the gas phase. In general dispersion effects “red” 
shift (a shift to lower energies) higher excited states more than 
lower one when compared to gas phase spectra, but are rea- 
sonably constant for different solvents”: ** *, 

In Table 5 we report the shifts in some n->pi* states of 
diazabenzenes**. The calculated shifts nicely reflect the ex- 
perimental observations. But we note that the shifts calculated 
for these systems in acetonitrile and water are very similar, for 
example, 31,296cm-1 and 31,368cm-1 for pyridazene respec- 
tively. Experimentally the difference is nearly 2500cm-1. 

The line “2H,O” designates a supermolecule of two wa- 
ters hydrogen bonded to each of the lone pairs of the basic 
nitrogens of the azabenzenes, as would be expected on chemi- 
cal grounds. These hydrogen bonds are of the “strong” nature 
discussed previously and cannot be expected to be accurately 
represented by a dielectric continuum model which assumes 
thermal averaging of orientations of the solvent in the field of 
the solute. Interestingly enough, this would lead to the predic- 
tion of a red shift (to lower energy) of these systems in a 
supermolecule. But note that when the entire super-molecule 
is embedded in the reaction field, nearly quantitative results 
are obtained. For example, in the pyridazene case discussed 
above a shift of 2630cm-1 is observed between water and 
acetonitrile, compared to the experimental value of 2490cm-1. 


Figure 9. 


The structure of the reaction center of R. Viridis, from reference 
35. 








R. viridis reaction center 


ies 
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A question naturally arises from these results. How sen- 
sitive is the calculated spectra of these compounds in water to 
the relative positions of the bound waters? This is the type of 
question that can be answered in a combined classical QM 
approach. Considering a Monte-Carlo calculation on the su- 
per-molecule in the presence of a dielectric continuum yield 
for pyrimidine a shift of 2600cm-1 relative to iso-octane as a 
solvent (experiment is about 2700 cm-1) and a width of the 
band of about 1300cm-1. This width comes from sampling the 
different possible conformations of the two waters bound to 
pyrimidine. No doubt the presence of more water molecules 
treated explicitly, either quantum mechanically or through 
classical force field, would yield a greater width. The results 
we obtain for this initial study are shown as a histogram in 
Figure 8. 


Photosynthesis: 


I mentioned two problems that propelled me from a naive 
quantum chemist, to one that was forced to acknowledge that 
many interesting problems simply do not occur in the gas 
phase. The fist of these was the mero-stability problem de- 
scribed above. The second had to do with our studies on the 
initial photo-chemical event in photosynthesis. 

To say that photosynthesis is an interesting reaction (or 
series of reactions) is an understatement. Nearly all of the 


biomass on earth is the result of photosynthesis. This reaction 
is the major one-almost unique-for fixing carbon atoms in 
arrangements more suitable for life as we know it. To accom- 
plish this, photosynthesis must destroy the strong bonds found 





Figure 10. 


The structure of an idealized dimer of R. Viridis, see reference 
37. The large central balls represent magnesium atoms. The 
light smaller balls represent nitrogen atoms, and the darker 
ones carbon atoms. The hydrogen atoms are not shown 
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in the carbon dioxide that is found in the atmosphere. The 
energy to do this comes from the sun. Light is absorbed, and 
produces sugars and starches from water and carbon dioxide. 
The overall stoichiometry for plant photosynthesis is given 
below: 


hv 
nCO, + nH,0 — (CH,0)n + nO, 


A lot of wonderful problems are associated with this 
reaction, and much is known about the chemistry that follows 
the initial photochemical event. What is uncertain is what 
happens during the first second after the absorption of a 
photon. This is an unlikely problem to be tackled by quantum 
chemistry until recently, when the structure of a photosynthetic 
bacteria, R viridis, was obtained by Huber, Deisssenhofer and 
Michel**. The essential part of the “reaction center” RC, is 
given in Figure 9. Experiments indicate that within 1 to 3 
picoseconds after the energy is received in the RC, in a dimer 
of chlorophyll, Figure 10, an electron is transferred from this 
“special pair” to bacterio-pheophytin, a magnesium free bac- 
terio-chlorophyll. 

So we*’:** as well as others **, embarked on applying 
quantum chemistry to this structure of over 500 atoms and 
1500 electrons, to examine the nature of the ground states, and 
the various excited states. These calculations clearly indicated 
that the lowest excited state of the system was one localized 
to the special pair, P*, and that the formation of the dimer itself 
was sufficient to insure this**. The charge transfer states, 
however, were calculated to lie some 8000cm-1 (leV or about 
23kcal/mol.) higher in energy than the absorbing state of the 
dimer. It thus is difficult to understand how to transfer energy 
from P* to P*H’, the state in which the pair has surrendered an 
electron (P-->P*) and the bacterio-pheophytin (H-->H’) has 
accepted it. 

But there is some 17 Angstroms from the center of the 
pair to the distant pheophytin. The transfer of an electron this 
distance generates a dipole moment of nearly 82 Debye (the 
dipole moment of water, for example, is 1.8 Debye). It is naive 
to expect that the protein material that surrounds the RC does 
not notice this large electric moment and respond to it! Using 
the SCRF theory for excited states as described above suggests 
that the largest term for the stabilization of this charge sepa- 
rated state is proportional to the dipole moment squared. 

SCRF calculations yield the results given in Table 6. The 
P* state shifts to slightly lower energies, while the P*H state 
“red” shifts some 7000cm-1, and now lies but 1200 cm-1 
above P*. The crossing from P* into the important P*H state, 
the precursor to all this chemistry, could easily occur through 
relaxation of the geometries of the 6 chromophores that make 
up the RC, or through further relaxation of the dielectric media 
through positional changers (rotations, translations and vibra- 
tions, see above). An important observation in this studies is 
that excited charge-transfer states from the pair to the auxiliary 








Figure 11. 


The simulated spectrum of CH2O in water. The quantum 
mechanics is from ZINDO and includes only the atoms of 
CH20 in the presence of 99 water molecules, treated as scaled 
point charges, see text. The conformations were generated 
by the Mobi program, reference *". 
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bacterio-chlorophyll’s, B, and from B to H that are near in 
energy to the P*H state before application of the reaction field, 
are considerably higher in energy after, and are too high in 
energy to take part in the electron transfer process itself, a 
matter of some controversy in the past. This observation is 
easy to understand: the dipole moments generated by the states 
P*B and B*H are much smaller than that of P*H™ (about half 
as large) and the stabilization is roughly proportional to the 
dipole moment squared. the actual calculations show P*H” 
stabilized by 7100cm-1, P*B” by only 3200cm-1. 

Now the RC is not surrounded by a homogeneous dielec- 
tric material. It is surrounded by real atoms in reasonably fixed 
positions. We have modeled all these atoms and their structure 
with cyclohexane treated as a dielectric continuum! We have 
tried to examine this crude model by adding more and more 
of the molecules present until we have run out of reasonable 
computing resources: i.e., extended the supermolecule of Fig- 
ure 1. The results we obtained from the simple cyclohexane 
model have not changed. 

More recently, Marchi, Gehlen, and Newton“! have per- 
formed statistical mechanical calculations on the RC and 
much of its surrounding proteins using our gas phase com- 
puted values to align the excited states correctly. The “solvent 


effect” is then generated through examining the motions of the 
various atoms included in the simulation. Their results are very 
similar-although more detailed-to those we obtained for the 
energetics. Their results, however, also include the dynamics 
necessary to show that our results can actually lead to the 
observed electron transfer in under 3 picoseconds, and in 
particular that the P*B” and B*H' states are not needed for the 
observed dynamics*”*!, 


Preliminary Results from 
Combined Models: 


In this section I would like to compare results obtained 
from the SCRF model with those obtained from a model in 
which a solute is embedded in a solvent and then the spectros- 
copy of this model examined using a molecular Monte-Carlo 
method. This model should, as discussed above, approach the 
correct answer if we could treat the ensemble quantum me- 
chanically. We cannot do this. The system examined is formal- 
dehyde in water. The potentials and the calculated relative 
orientations of the 99 water molecules and formaldehyde 
come from the Mobi program“ utilizing box boundary condi- 
tions. At every tenth geometry that we obtain from the simu- 
lation consistent with room temperature (298K) we calculate 
the spectroscopy of formaldehyde in the presence of 99 waters. 
The spectroscopy is calculated treating formaldehyde quan- 
tum mechanically and treating the waters as point charges that 
reproduce the dipole of water in water using the ZINDO 
program“. The screening of the waters is checked by includ- 
ing more and more of the waters that surround the acetone in 
the quantum mechanics for several test conformations. A 
histogram of results occur in Figure 11. The average shift from 
the gas phase spectra of acetone is 1200cm". This might be 
compared with a 1900cm” shift calculated by Blair, Levy, and 
Krogh-Jespersen™ “°. We believe the differences have to do 
with different charge scaling schemes. 

In Table 7 we present the results we obtain for a formal- 
dehyde-water super-molecule in the reaction field of water. 
Two observations are very interesting! The first is that these 
far simpler calculations yield the same average shift, 1200cm” 
' The second is that the bound water molecule is responsible 
for about 1/3 of the average shift. This is also the conclusion 
that we reach in the full Monte-Carlo calculation if we include 
only one (or even two) of the nearest neighbor water(s)*“. 


Conclusions: 


This paper is long enough! Hopefully what is portrayed 
is the wealth of information that can be obtained from a 
relatively simple model-the self consistent reaction field 
model. There is now little question that this model bridges 
most of the gap between gas phase quantum chemistry, and 
what is observed in solution. By its very nature, however, we 
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model reality. This SCRF model is unlikely ever to be truly 
quantitative. More accurate models are available, such as 
Monte-Carlo or molecular dynamics models using classical 
potentials that can be used together with quantum chemistry 
as discussed briefly in this manuscript. This model should be 
quite accurate for examining situations in which the coupling 
between solute and solvent is weak, and the “solute” is defined 
as containing all the parts that must be treated quantum me- 
chanically-parts reacting, parts absorbing light, parts contain- 
ing strong interactions as described previously. A fully 
quantum mechanical Monte-Carlo or molecular dynamics 
model is not available today, nor likely in the near future. But 
I am hopeful that many of the interesting chemistry and 
physics questions that arise will find their solutions in these 
simpler models. 
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